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1.  Executive  Summary 


This  report  presents  a  comprehensive  summary  of  the  research  work  performed  under  ONR 
funding,  grant  No.  N000 14-04- 1-0 167  over  the  period  of  June  1,  2004  to  June  1,  2005. 

The  contributors  to  the  research  over  the  span  of  the  funding  period  are  Prof.  Yimin  Zhang 
(PI),  Prof.  Moeness  Amin  (Co-PI),  Dr.  Genyuan  Wang  (Post-Doctoral  Fellow),  and  Mr.  Gjergji 
Kurti  (Graduate  Student)  from  Villanova  University.  Ms.  Morgan  Watson  from  NAVSEA 
Philadelphia  has  worked  closely  with  the  research  team  at  Villanova  University  and  has  provided 
valuable  insights  into  several  issues  vital  to  the  progress  and  advances  in  research. 

The  results  obtained  in  this  project  have  been  disseminated  in  three  submitted  journal  articles 
and  seven  peer-reviewed  conference  papers.  Section  2  lists,  with  complete  citations,  all 
publications  acknowledging  the  ONR  support.  This  list  is  also  used  to  point  the  reader  to  the 
technical  material  and  the  pertinent  papers  in  the  report  underlying  each  subsection.  Section  3  of 
this  final  report  includes  selected  key  papers,  representing  our  efforts  in  different  areas  covering 
cooperative  diversity,  space-time  coding,  and  channel  modeling,  and  spatial  correlation  study  in 
the  sea  surface  environment.  Each  of  these  papers  is  cast  as  a  subsection  with  its  own  abstract, 
introduction,  figure  numbers,  equation  numbers,  conclusions,  and  references. 

1.1  Main  Contributions 

The  fundamental  objective  of  this  research  project  is  to  develop  novel  cooperative  wireless 
network  schemes  and  practical  coding  designs  that  provide  unmanned  surface  vehicles  (USV) 
swarm  reliable  wireless  links  with  the  overall  data  throughput  and  the  communication  quality 
significantly  exceeding  those  of  existing  non-cooperative  systems.  Towards  this  end,  we  have 
made  contributions  in  three  areas:  (1)  We  have  developed  a  cooperation  protocol  that  effectively 
constructs  distributed  space-time  codes  for  cooperative  diversity  systems.  Depending  on 
application  environments,  the  proposed  protocol  can  be  used  to  construct  different  distributed 
space-time  codes.  For  example,  high  data  rate  distributed  codes  can  be  used  to  improve  the 
communication  efficiency  when  the  channel  state  information  (CSI)  is  available  at  the  receivers, 
whereas  differential  distributed  space-time  codes  can  be  used  to  avoid  channel  estimations  in  the 
absence  of  CSI.  (2)  We  have  developed  full-rate,  large  diversity  product  space-time  codes  with 
non-vanishing  minimum  determinant.  The  proposed  space-time  codes  are  an  important 
component  in  high-efficiency  distributed  space-time  code  construction.  (3)  We  have  studied  the 
channel  modeling  and  spatial  correlation  between  different  antennas  in  sea  surface  environments. 
The  results  are  important  in  assessing  the  performance  of  cooperative  diversity  systems  in  sea 
surface  environments  as  well  as  in  designing  array  antenna  configurations. 

1.2  Cooperative  Diversity  Techniques 

We  have  developed  a  protocol  for  cooperative  systems  that  permits  effective  use  of  various 
multiple-input-multiple-output  (MIMO)  space-time  coding  schemes.  Based  on  high-rate  MIMO 
space-time  codes  and  differential  space-time  codes,  the  proposed  protocol  allows  the 
development  of  high-rate  distributed  space-time  codes  and  differential  distributed  space-time 
codes,  respectively,  depending  on  the  channel  state  information  (CSI)  availability,  for  the  use  of 
USV  network  systems.  A  multiuser  scenario  is  considered  where  all  users  have  their  own 
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information  to  send,  and  desire  to  cooperate  with  each  other  in  order  to  maintain  reliable 
communication  links  and  to  send  the  information  to  the  receiver  at  a  high  data  rate.  When  the 
CSI  is  available  at  the  receivers,  this  protocol  enables  the  use  of  high-rate  full-diversity  space- 
time  codes  for  the  construction  of  high  efficiency  distributed  space-time  codes.  On  the  other 
hand,  when  CSI  is  unavailable  at  the  receivers,  differential  distributed  space-time  modulations 
can  be  constructed  based  on  the  proposed  protocol. 

Publications  [1,  2,  3,  7,  8,  9]. 

1.3  Space-Time  Code  Development 

Designs  of  space-time  code  with  full  rate,  large  diversity  product,  and  non-vanishing  minimum 
determinant  of  codewords  continue  to  attract  great  attentions.  However,  in  most  available 
determinant  non-vanishing  space-time  codes  for  three,  four  and  six  transmit  antennas,  the 
average  powers  at  each  layer  are  different,  resulting  high  peak  to  average  power  ratio.  In  this 
paper,  a  new  cyclic  algebraic  space-time  design  scheme  is  proposed  and  the  optimal  codes  in  this 
class  are  provided  by  using  some  specific  cyclic  field  extensions.  The  class  of  cyclic  algebraic 
space-time  codes  not  only  include  the  available  non-vanishing  determinant  cyclotomic  space- 
time  codes  for  three,  four  and  six  transmit  antennas,  but  also  have  the  desirable  property  that  the 
optimal  codes  can  be  achieved  with  the  same  average  power  at  each  layer. 

Publications  [4,  10]. 

1.4  Channel  Modeling  and  Spatial  Correlation  Study  in  Sea  Surface  Environment 

We  have  studied  the  propagation  environment  and  technology  applicable  to  multi-sensor 
communications  in  USV  networks  located  in  a  rough  sea  environment.  We  first  investigated  the 
statistical  model  of  the  propagation  environment,  where  we  provided  the  analytical  expressions 
of  the  wave  height,  the  probability  of  having  line-of-sight  (LOS)  between  two  USVs  in  terms  of 
the  weather  condition  and  antenna  height,  and  the  spatial  correlation  of  the  propagation  channels 
related  to  different  antennas.  We  have  then  investigated  the  fading  reduction  capability  of 
multiple-input-multiple-output  (MIMO)  systems  and  the  performance  in  general  Ricean  channels 
assuming  different  spatial  channel  correlation  values. 

Publications  [5,  6], 
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Cooperative  Wireless  Network  for  Unmanned  Surface  Vehicles 

Yimin  Zhang',  Genyuan  Wang' ,  Morgan  Watsori,  and  Moeness  G.  Amiri 

t  Center  for  Advanced  Communications 
Villanova  University,  Villanova,  PA  19085 
e-mail:  {yimin. zhang, genyuan. wang, moeness .amin}@villanova. edu 

http:  /  /  www.villanova.edu/ 

i  NAVSEA 

4700  S.  Broad  Street,  Philadelphia,  PA  19112-1403 
e-mail:  WatsonMAOnswccd.navy.mil 

Abstract 

As  part  of  Future  Naval  Capabilities  (FNC),  the  Navy  is  envisioned  to  be  a  lean,  distributed  force 
of  reconfigurable  and  multi-mission  capable  platforms.  At  the  heart  of  this  vision,  cascading  unmanned 
surface  vehicles  (USVs)  into  a  “SWARM”  of  connectivity  nodes  offers  great  flexibility  and  modularity. 
The  USVs  must  be  able  to  provide  a  network  of  communication  data  relays,  a  common  tactical  picture 
and  mapping  for  situational  awareness.  The  vehicles  must  also  possess  robust  communications  resistant  to 
jamming,  fading  and  interference  while  acquiring  a  high  speed  data  exchange  for  processing  and  fusion  of 
information.  Reliable  wireless  network  with  high  data  throughput  is  essential  in  various  military  wireless 
communication  systems.  While  the  multiple-input  multiple-output  (MIMO)  communication  schemes  have 
the  potential  to  exceed  the  single-input  single-output  (SISO)  Shannon  capacity  bound,  equipping  a  large 
number  of  antennas  in  a  USV  is  practically  impossible.  To  increase  the  diversity  gain  beyond  the  number 
of  physical  antennas,  cooperative  wireless  networks  that  exploit  cooperation  among  multiple  terminals 
become  essential.  This  paper  proposes  novel  cooperative  wireless  networking  schemes  and  practical  coding 
designs  that  provide  the  USV  Swarm  reliable  wireless  links  with  the  overall  data  throughput  and  the 
communication  quality  significantly  exceeding  those  of  the  existing  non-cooperative  systems. 

I.  Introduction 

As  part  of  Future  Naval  Capabilities  (FNC),  the  Navy  is  envisioned  to  be  a  lean,  dis¬ 
tributed  force  of  reconfigurable  and  multi-mission  capable  platforms.  At  the  heart  of 
this  vision,  cascading  unmanned  surface  vehicles  (USVs)  into  a  “SWARM”  of  connectiv¬ 
ity  nodes  offers  great  flexibility  and  modularity.  USVs  are  on  the  threshold  of  playing 
key  roles  in  the  digital  battlespace.  Within  the  context  of  a  larger  group  of  unmanned 

This  work  is  supported  in  part  by  the  ONR  under  Grant  No.  N00014-04-1-0617. 
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systems,  US  Vs  have  been  in  use  for  years  or  are  in  the  late  stages  of  development.  Rela¬ 
tively,  Unmanned  Air  Vehicles  (UAVs)  are  commonplace  in  many  military  operations  as 
both  weapons  and  reconnaissance  platforms.  Unmanned  Ground  Vehicles  (UGVs)  are  be¬ 
ing  developed  for  high-risk  operations  such  as  surveillance,  while  Unmanned  Underwater 
Vehicles  (UUVs)  are  new  and  in  development. 

US  Vs  are  envisioned  to  be  a  multi-mission,  rapidly  configurable,  long  endurance  un¬ 
manned  platform  providing  offboard  sensor  and  weapon  capability  to  perform  critical 
missions  like  Intelligence,  surveillance,  reconnaissance,  Anti-Submarine  Warfare,  mine 
countermeasures,  communications  and  navigation.  The  technical  challenges  during  these 
USV  missions  are  to  provide  optimal  performance  while  subjected  to  adverse  motions, 
multi-vehicle  control,  data  fusion,  obstacle  avoidance  as  well  as  robust  and  reliable  com¬ 
munications.  More  importantly,  an  autonomous  environment  where  electromagnetic  ra¬ 
diation  is  increased,  brings  with  it  the  possibility  of  fading,  jamming  and  interference. 
This  paper  addresses  methods  for  maximizing  communication  and  reducing  fading  and 
interference.  In  particular,  this  paper  focuses  on  novel  high-rate  full-diversity  cooperative 
diversity  schemes  for  wireless  communication  systems  through  the  use  of  space-time  codes. 
The  proposed  systems  enable  simultaneous  transmission  of  data  from  all  cooperative  ter¬ 
minals,  providing  high  diversity  gain  and  potentially  high  data  rate  beyond  the  limit  of 
non-cooperative  systems. 

In  combating  channel  impairment  problem  caused  from  multipath  fading,  receive  diver¬ 
sity  techniques  using  different  space  branches  at  the  receiver  have  been  shown  effective 
for  many  years.  When  the  channel  is  fast  time-varying,  time  diversity  is  also  an  option. 
Recently,  transmit  diversity  schemes  add  the  capability  of  diversity  to  the  transmitter  to 
spread  the  information  across  the  transmit  antennas  [1],  [2].  As  a  result,  the  information 
capacity  of  wireless  communication  systems  increases  dramatically  by  employing  multiple 
transmit  and  receive  antennas  [3].  A  wireless  system  consisting  of  multiple  transmit  and 
receive  antennas  is  referred  to  as  a  multi-input  multi-output  (MIMO)  system,  whereas  a 
system  with  multiple  transmit  antennas  and  a  single  received  antenna  is  referred  to  as  a 
multi-input  single-output  (MISO)  system.  When  the  correlation  among  different  channel 
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path  gains  is  low,  a  MIMO/MISO  system  can  achieve  higher  diversity  gain  and  higher 
data  rates,  compared  to  a  single  transmit  antenna  system,  without  increasing  premium 
bandwidth  or  power  [3],  [4], 

When  the  channel  state  information  (CSI)  is  available  at  the  receiver,  several  space- time 
codes  have  been  developed  to  achieve  high  data  rate,  while  maintaining  a  high  diversity 
gain  [14],  [15],  [16],  [17],  [18].  In  particular,  a  space-time  code  is  said  to  achieve  full-rate 
when  the  data  rate  (i.e.,  the  number  of  transmitted  symbols  per  time  slot)  reaches  the 
number  of  transmit  antennas.  Use  of  the  full-rate  full-diversity  space-time  codes  allows 
simultaneous  signal  transmission  without  reducing  date  rates. 

In  certain  situations,  for  example,  when  the  fading  channel  environment  is  rapidly  time- 
varying,  the  CSI  may  not  available  or  cannot  be  reliably  estimated.  Differential  space- 
time  codes  are  useful  which  imply  differential  coding  between  two  adjacent  space-time 
code  matrices.  In  this  case,  the  decoding  at  the  receiver  is  independent  of  the  underlying 
channels  [22],  [23],  [24],  [25],  [26],  [27],  [28]. 

It  may  often  be  impractical,  however,  for  a  transmitter  in  the  underlying  USV  networks 
to  be  equipped  with  multiple  antennas.  For  example,  the  physical  limitation  of  a  USV 
makes  it  only  ready  to  host  a  single  antenna,  particularly  when  the  frequency  is  relatively 
low.  In  addition,  depending  on  the  fading  characteristics,  for  example,  when  shadowing 
effect  occurs,  the  use  of  multiple  antennas  may  not  be  effective  to  combat  the  attenuation. 
To  increase  the  reliability  and  throughput  of  a  system  in  such  situations,  cooperative 
diversity  exploiting  cooperation  among  multiple  terminals  is  a  promising  techniques.  It 
was  first  introduced  in  [5],  [6]  and  further  developed  by  others  (e.g.,  [7],  [8],  [9],  [10]).  In 
a  multi-user  environment,  cooperative  diversity  allows  mobile  terminals  to  achieve  high 
transmit  diversity  gain  beyond  the  limitation  of  the  number  of  transmitter  antennas  at  each 
terminal.  It  has  been  shown  that  cooperative  diversity  can  improve  the  communication 
capacity  and  enhance  the  robustness  of  .  wireless  link  when  a  single  channel  is  not  reliable. 

In  this  paper,  we  develop  a  protocol  for  cooperative  systems  that  permits  effective  use 
of  various  MIMO  space-time  coding  schemes  [11],  [12].  Based  on  the  aforementioned 
high-rate  MIMO  space-time  codes  and  differential  space-time  codes,  the  proposed  pro- 
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tocol  allows  the  development  of  high-rate  distributed  space-time  codes  and  differential 
distributed  space-time  codes,  respectively,  depending  on  the  CSI  availability,  for  the  use 
of  USV  network  systems.  A  multiuser  scenario  is  considered  where  all  users  have  their 
own  information  to  send,  and  desire  to  cooperate  with  each  other  in  order  to  maintain 
reliable  communication  links  and  to  send  the  information  to  the  receiver  at  a  high  data 
rate. 

II.  Cooperative  Diversity 

A.  Concept 

To  illustrate  the  concept  of  cooperative  diversity  in  a  wireless  network,  consider  a  simple 
model  as  depicted  in  Fig.  1(a).  The  user  cooperates  with  other  users  and  serves  as  a  relay 
terminal  for  them.  Therefore,  each  transmit  user  receives  an  attenuated  and  noisy  version 
of  the  partner’s  transmitted  signal  and  relays  it  to  the  destination  or  other  relays.  The 
destination  receives  a  noisy  version  of  the  sum  of  the  attenuated  signals  from  all  users. 
Denote  the  number  of  users  as  M.  The  received  signal  at  the  destination  during  one 
symbol  period  t  is  expressed  in  baseband  model  as 

M 

Vo(t)  =  hn&iit)  +  no(£),  (1) 

i 

whereas  the  signal  received  at  relay  k  is 

Vk{t)  =  Y2  hikXi{t )  +  nk{t).  (2) 

i^k 

We  used  subscript  0  for  the  destination  terminal  for  notational  convenience.  In  (1)  and 
(2),  Xi(t)  is  the  signal  transmitted  by  user  i,  for  i  —  1, 2,  •  •  • ,  M,  and  Uk(t)  are  the  additive 
channel  noise  terms  at  the  destination  ( k  =  0)  and  the  relay  receivers  (A;  =  1,2,-*-,  M), 
respectively.  It  is  assumed  that  a  terminal  does  not  receive  signal  from  other  parties  when 
it  is  transmitting.  The  fading  channels,  h^,  remain  constant  over  a  period  of  several  time 
slots  (the  length  of  one  or  two  codewords  depending  on  the  protocols),  and  when  observed 
over  time,  they  form  independent  stationary  ergodic  stochastic  processes,  resulting  in 
frequency  non-selective  fading. 


8 


(a)  System  model.  (b)  Data  transmission  through  relays. 

Fig.  1.  System  model  illustration. 

The  cooperation  process  can  be  divided  into  two  phases  (refer  to  Fig.  1(b)).  In  the  first 
phase  (broadcast  phase),  the  information  is  transmitted  from  a  source  user  to  the  relay 
terminals,  and  the  destination  may  also  receive  a  copy  of  the  same  information.  In  the 
second  phase  (relay  phase),  the  relay  terminals  transmit  the  signal  to  the  destination. 

B.  Amplify-and-Forward  and  Decode- and-Forward  Algorithms 

Depending  on  how  the  relay  terminals  relay  the  signals  from  other  terminals,  there 
are  two  major  algorithms,  namely,  amplify-and-forward  and  decode-and-forward.  In  the 
amplify-and-forward  algorithm,  a  relay  terminal  amplifies  the  attenuated  and  noisy  signals 
it  receives  and  retransmits  them  to  the  destination  and  other  possible  relay  terminals.  The 
operation  at  the  relay  terminal  is  limited  to  amplification  and,  in  some  cases,  some  simple 
computations  such  as  complex  conjugation  which  is  required  in  certain  orthogonal  codes 
including  Alamouti’s.  On  the  other  hand,  when  the  decode-and-forward  algorithm  is  used, 
the  information  is  first  decoded  at  a  relay  terminal,  and  then  retransmitted  from  the  delay 
terminal  after  proper  coding.  At  the  expense  of  introducing  additional  complexity  at 
the  relay  terminals,  the  decode-and-forward  algorithm  allows  the  removal  of  relay  noise 
(i.e.,  the  noise  the  relay  terminals  receive),  and  provides  the  flexibility  of  encoding  the 
information  at  the  relay  phase  so  that  higher  spectral  efficiency  can  be  achieved  [8],  [21]. 

It  is  clear  that  the  decode-and-forward  algorithm  requires  higher  complexity  and,  in  some 
applications,  may  cause  security  concerns  because  a  relay  terminal  decodes  the  information 
it  relays  and,  therefore,  it  is  technically  possible  for  a  relay  terminal  to  monitor  other  users’ 
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information.  It  is  contrast  to  the  amplify-and-forward  algorithm,  where  a  relay  terminal 
does  not  decode  the  information  it  relays.  In  addition,  the  decode-and-forward  algorithms 
requires  the  CSI  for  effective  decoding,  unless  the  modulation  schemes  used  do  not  require 
CSI  for  decoding  (such  as  the  differential  modulation  schemes  to  be  discussed  later). 

C.  Protocols 

Several  cooperative  diversity  protocols  have  been  investigated.  The  early  protocols  used 
repetition-based  approaches,  which  is  illustrated  in  Fig.  2(a).  After  the  source  broadcasts 
its  information  to  the  destination  and  the  relays,  all  relays  repeats  this  information  in 
a  sequential  order.  Therefore,  a  diversity  gain  of  up  to  M  is  achieved  at  the  expense 
of  reducing  the  degrees  of  freedom  by  a  factor  of  M.  Arbitrary  permutations  of  these 
allocations  in  time  and  frequency  does  not  alter  the  performance,  as  long  as  causality  is 
preserved  and  each  of  the  subchannels  contains  a  fraction  of  1/M2  of  the  total  degrees  of 
freedom  (DOFs)  in  the  channel.  Such  protocols  does  not  require  decoding  and  encoding  at 
the  relay  terminals  and,  therefore,  is  suitable  for  the  amplify-and-forward  algorithm.  The 
most  significant  disadvantage  of  such  protocol  is  their  low  throughput  rates.  As  is  clear 
from  Fig.  2(a),  transmission  between  different  nodes  at  each  time  slot  utilizes  a  fraction 
of  1/M  of  the  total  DOFs  in  the  channel. 

Recently,  more  effective  protocols  have  been  developed  to  take  advantages  of  the  ad¬ 
vances  of  MIMO  space-time  codes.  For  example,  in  [8],  [21],  proposed  space-time  coop¬ 
eration  protocols  illustrated  in  Figs.  2(b)  and  (c).  These  protocols  provide  more  effective 
bandwidth  use,  particularly  when  the  number  of  cooperative  users  is  large.  These  protocols 
allow  the  use  of  a  fraction  of  up  to  1/2  of  the  total  DOFs  in  the  channel.  The  protocols  re¬ 
quire  decoding  of  other  users’  information  and,  therefore,  should  be  incorporated  with  the 
decode-and-forward  algorithm.  As  addressed  earlier,  the  use  of  the  decode-and-forward 
algorithm  not  only  significantly  increases  the  complexity  of  the  terminals  that  is  highly 
undesirable,  it  may  also  raise  security  concern  in  some  applications. 


10 


Ch.  1 
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1  repeats  2 

• 

*  .  4 

; 
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(a)  repetition-based  diversity  scheme 


(b)  space-time  cooperation  scheme  I 


Frequency 

Ch.  1 
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2,3,.. .,M  repeat  1 

Ch.2 

2  transmits 

1,3,...,M  repeat  2 

: 

j 

i 

Ch.  M 

M  transmits 

1,2,...,M-1  repeat  M 

Time 

(c)  space-time  cooperation  scheme  II 


Fig.  2.  Cooperative  diversity  schemes. 


III.  Proposed  High-Rate  Full-Diversity  Cooperation  Protocol 

Figure  3  shows  the  proposed  cooperation  protocol.  There  are  M  users  present,  each 
is  equipped  with  a  single  antenna  that  can  be  used  for  transmit  and  receive.  When 
considering  user  i  as  the  source  user,  the  other  M  —  1  ones  act  as  the  relay  users  for  user 
i.  At  the  source  user,  a  set  of  M  x  M  multi-layer  cyclotomic  space-time  codewords,  X*, 
are  formed  from  the  source  information  sequence.  Before  the  codeword  is  modified  for  the 
cooperative  diversity  use,  the  data  rate  is  full  (i.e.,  R  =  M ).  However,  because  the  M  -  1 
relays  do  not  have  the  information  to  be  relayed,  the  source  user  must  first  broadcast  this 
information  to  the  relays. 

In  this  protocol,  through  the  first  M  —  1  time  blocks,  the  Ah  row  of  the  M  x  M 
codeword  is  transmitted  during  the  zth  time  block  from  each  source  transmitter,  where 
i  =  1,  -  •  • ,  M  —  1.  In  the  Mth  time  block,  the  full  codeword  is  transmitted  from  the  M 
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virtual  antennas  (i.e.,  one  source  transmit  antenna  and  M  —  1  relay  antennas). 

The  proposed  method  is  amendable  for  the  use  of  both  amplify-and-forward  and  decode- 
and-forward  algorithms.  When  the  former  one  is  considered,  this  protocol  does  not  require 
decoding  and  encoding  at  the  relay  terminals  while  it  takes  the  advantage  of  high  data 
rate  and  high  diversity  gain.  In  this  case,  each  of  the  M  —  1  relay  terminals  receives  one 
different  row  of  the  codeword  transmitted  from  the  source  user  through  the  broadcast 
phase  (i.e.,  during  the  first  M  —  1  time  blocks),  and  the  row  is  retransmitted  in  the  relay 
phase  which  spans  the  Mth  time  block.  The  use  of  amplify-and-forward  algorithm  does 
not  require  the  relays  to  have  the  CSI.  It  is  evident  that,  in  the  proposed  scheme,  unit 
degree  of  freedom  is  achieved,  compared  to  1  /M  (repetition-based  diversity  scheme  in  Fig. 
2(a))  and  1/2  (space-time  cooperation  schemes  I  and  II  in  Figs.  2(b)  and  (c))  decode-based 
diversity  scheme  in  Fig.  2(a))  in  the  aforementioned  space-time  cooperation  schemes. 

When  the  decode-and-forward  algorithm  is  considered,  the  M  —  1  relay  users  can  decode 
the  signal  from  the  information  transmitted  through  the  first  M  —  1  time  slots.  The  signal 
will  be  retransmitted  after  encoding  during  the  Mth  time  slot  either  from  all  the  M  users 
(space-time  cooperation  scheme  I)  or  from  the  M— 1  relays  (space-time  cooperation  scheme 
II).  Unless  differential  or  other  non-coherent  coding  schemes  (to  be  introduced  later)  are 
exploited,  the  use  of  decode-and-forward  algorithm  requires  each  relay  to  have  the  CSI 
of  the  channel  between  the  source  and  itself.  Such  decoding  and  encoding  can  effectively 
eliminate  the  relay  noise  when  the  channel  between  the  source  user  and  the  relay  terminals 
are  reliable  to  the  extent  that  the  information  can  be  correctly  detected  at  each  relay. 

It  is  emphasized  that,  because  the  codeword  X,  is  full  rank,  the  relay  users  can  decode 
the  signal  with  a  full  or  partial  set  of  the  M  rows  of  the  M  x  M  codeword.  This  permits 
flexible  protocol  design  to  further  increase  the  data  rates  by  truncating  the  broadcast  phase 
when  the  channels  between  the  source  user  and  the  relay  terminals  are  highly  reliable.  This 
comes  at  the  expense  of  requiring  additional  complexity  to  both  the  physical  (PHY)  and 
media  access  control  (MAC)  layers  in  order  to  monitor  the  channel  quality  and  decoding 
performance,  and  to  adaptively  control  the  resource  allocations  in  the  broadcast  phase. 

From  a  2  x  2  MIMO  space-time  codeword  Xi,  the  protocols  equivalently  transmit  the 
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1st  row  ofX2 

2nd  transmits 
2nd  row  of  X2 

transmit 
codeword  X2 

* ' 

. 

• 

Ch.  M 

M  transmits 

1st  row  ofXM 

M  transmits 

2nd  row  of  X  ^ 

transmit 
codeword  XM 

Fig.  3.  The  proposed  space-time  cooperation  scheme. 


following  2x4  codeword  X',  i.e., 


*i(l)  Xi( 2) 
rci(3)  Xi(4) 


v/ _  /Zi\x^)  Xi( 4)  Xi{  1)  Xi( 2) 

*  V9  [  0  0  Xi( 3)  *i(4)J  ’ 


(3) 


where  Xi(t)  is  a  replica  of  Xiit)  and  g'  is  a  power  normalization  factor.  When  the  amplify- 
and-forward  algorithm  is  used,  Xi(t)  is  a  complex  scaled  and  noisy  version  of  Xi(t).  The 
scale  comes  from  the  channel  attenuation  and  the  amplification  in  the  relay  terminal, 
whereas  the  channel  noise  is  added  before  the  amplification.  When  the  decode-and-forward 


algorithm  is  used,  the  scale  and  noise  can  be  removed,  given  the  channels  between  the 
source  and  the  relay  terminals  are  reliable. 

When  we  assume  Xi(t )  is  the  exact  replica  of  Xi(t),  we  can  equivalently  consider  the 


resulting  codeword  as, 


Xi( 1)  Xi( 2) 
Xi(  3)  Xi(  4) 


v/  _  [Z,  I  Si(3)  Xi{ 4)  Xi(  1)  Xi( 2) 

1  *  9  q  q  Xi(3)  Xi(4) 


(4) 


A  code  like  X'  in  (4)  is  referred  to  as  the  distributed  space-time  code  corresponding  to 
the  space-time  code  X*. 


IV.  High-Efficiency  Space-Time  Cooperation  Schemes 


In  this  section,  we  consider  effective  space-time  cooperation  schemes  when  the  CSI  is 
available  at  the  receivers.  We  first  review  the  general  criteria  of  MIMO  space-time  code 
design  and  some  class  of  high-rate  full-diversity  space-time  codes.  We  then  propose  a  new 
space-time  cooperation  protocol  that  can  effectively  incorporate  the  high-rate  full-diversity 
space-time  codes  and  is  applicable  to  both  amplify-and-forward  and  decode-and-forward 
algorithms. 


13 


A.  Space- Time  Codes  in  MIMO  Systems 

Assume  a  transmitter  that  has  Nt  transmit  antennas,  and  consider  that  the  following 
Nt  x  Nt  square  codeword 


is  transmitted.  In  a  codeword,  the  rows  of  the  coding  matrix  stand  for  “space”  (i.e., 
antennas),  and  the  columns  stand  for  “time”.  When  codeword  E  is  transmitted,  the 
probability  that  a  maximum-likelihood  receiver  decides  erroneously  in  favor  of  another 
signal  codeword  E  in  independent  fading  channels  is  bounded  by  [4] 

P(C  — ►  E)  <  (j[. Aij  (Es/ANo)~rNr ,  (5) 

where  r  is  the  rank  of  DCe  =  C  —  E,  A *,  i  =  1,  •  •  • ,  r  are  the  r  nonzero  eigenvalues  of 
DceDce,  Nr  is  the  number  of  receive  antennas,  Es  is  the  average  transmitted  symbol 
energy,  and  N0  is  the  noise  variance  of  N0/2  in  each  of  its  two  dimensions. 

There  are  several  criteria  for  MIMO  space-time  code  design.  The  most  commonly  used 
ones  are  the  rank  criterion,  determinant  criterion,  and  symbol  rate  criterion.  In  (5),  the 
rank  r  determines  the  order  with  which  the  codeword  error  rate  (CER)  decreases  with  the 
input  signal-to-noise  ratio  (SNR)  when  the  SNR  goes  to  infinite.  Maximizing  the  diversity 
product  fj£=1  A i  also  reduces  the  CER.  The  symbol  rate  is  defined  as  R  =  M/T  symbols 
per  channel  use  (PCU)  if  an  information  sequence  of  M  information  symbols  is  mapped 
to  a  space-time  codeword  X  with  T  time  slots  (i.e.,  each  codeword  is  a  Nt  x  M  matrix). 
It  is  said  that  a  space-time  code  achieves  full  rate  if  its  symbol  rate  is  R  =  Nt  symbols 
PCU  for  a  transmit  array  of  Nt  antennas. 

B.  High- Rate  Full- Diversity  Space- Time  Codes 

With  these  design  criteria  in  mind,  several  space-time  codes  have  been  developed  that 
achieve  full  data  rate  and  full  diversity  gain.  For  example,  multi-layer  cyclotomic  space- 
time  codes  have  been  proposed  which  bear  high  or  even  full  data  rate  information  [13], 
[14],  [15],  [16],  [17],  [19], 
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Consider  a  simple  situation  where  two  transmit  antennas  are  present.  The  following 
two-layer  code  provides  data  rate  of  two  (i.e.,  full  rate  R  =  2), 


V  = 


r  (i)  (2) 

p\v\  p2v\ 

(2)  (1) 

Ypm  pivy 


(6) 


where  vf\  i  =  1, 2,  l  =  1,2,  are  the  ith  point  in  the  cyclotomic  lattice  of  the  Zth  layer,  and 
Pi  are  some  complex  numbers  (see  [19]  for  the  determination  of  values  of  pi).  For  each  of 
the  two  users,  i  =  1,2,  the  codeword  matrix  in  (6)  is  mapped  to  the  following  data  matrix 
to  be  transmitted, 


Xi(l)  Xi{ 2) 

Zi(3)  Xi(  4) 


Vgv  =  ^\plVh)  p2Vh 

L  P2^2  Pm 


where  g  is  a  factor  to  normalize  the  averaging  codeword  energy.  The  multi-layer  cyclotomic 
space-time  codes  are  full-rate,  and  are  shown  in  [19]  to  be  of  full  diversity  gain  and  optimal 
diversity  product.  These  codes  have  full-rate  and  full-diversity  features  and,  when  used 
instead  of  the  original  data  stream  in  an  independent  MIMO  fading  channel  environment, 
improvements  in  system  performance  and  power  efficiency  can  be  expected. 


V.  Differential  Cooperation  Schemes 

The  previous  section  considered  effective  space-time  cooperation  schemes  when  the  CSI 
is  available  at  the  receivers.  Such  information  is  obtained  usually  through  channel  esti¬ 
mation,  either  using  training  (pilot)  signals  or  utilizing  blind  methods.  However,  channel 
estimation  is  often  complicated  and  may  reduce  the  transmission  efficiency  if  pilot  sig¬ 
nals  are  used.  In  addition,  there  is  an  issue  of  reliability  and  practicality  if  the  channels 
experience  fast  fading.  This  section  develops  differential  cooperation  schemes  based  on 
the  differential  space-time  codes.  By  eliminating  the  needs  of  channel  estimation,  such 
schemes  are  useful  in  such  situations  in  which  the  CSI  are  either  unknown  or  cannot  be 
reliably  estimated. 

We  first  review  the  differential  space-time  codes  developed  for  MIMO  systems,  and  then 
show  the  feasibility  of  extending  them  to  the  space-time  cooperation  applications. 
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A.  Differential  Space-Time  Codes 

Differential  space-time  codes  are  useful  in  dealing  with  rapid  MIMO  fading  environ¬ 
ments.  Differential  coding  is  applied  between  two  adjacent  space-time  code  matrices  so 
that  the  decoding  at  the  receiver  is  independent  of  the  underlying  channels  [22],  [23],  [24], 
[25],  [26],  [27],  [28], 

For  example,  denote  Q  as  a  set  of  K  different  unitary  matrices  used  for  differential 
coding,  k  =  1,  •  •  • ,  K,  and 

Q  =  {Ufc,  k=  1,  •  •  • ,  K}  .  (8) 

The  information  steam  is  mapped  to  one  of  the  K  space-time  codewords  in  Q.  To  send 
the  message  G(t )  G  Q  at  time  f,  the  transmitter  sends  C (t)  where 

C(t)  =  C(t  —  l)G(t),  C{0)  =  I,  (9) 

with  C  (t  —  1)  denoting  the  matrix  transmitted  over  the  time  block  t  —  1  and  I  being 
the  identity  matrix.  The  differential  space-time  coding  scheme  works  well  in  time-varying 
channels,  provided  that  the  channels  do  not  change  significantly  over  two  adjacent  code 
matrices. 

B.  Differential  Space-Time  Cooperation  Schemes 

We  now  consider  the  feasibility  of  extending  the  differential  space-time  codes  to  coop¬ 
eration  networks.  For  simplicity,  we  consider  a  two-user  case  in  the  absence  of  relay  noise, 
and  user  2  is  a  terminal  to  relay  user  l’s  information.  The  following  2x4  distributed 
differential  space-time  code  C'(t)  is  generated  from  a  2  x  2  differential  space-time  code 
C(t): 

C'(f)  _  Cl(0  c2{t)  .  —  C3(^)  c4^)  cl(^)  02^)  fim 

C(t,-[c3(f)  c(t)J  ^  0  0  c,(()  C(t) J-  (10) 

Note  here  that  we  use  t  to  denote  the  time  index  for  a  codeword.  Corresponding  to  (9), 
the  relationship  between  two  adjacent  distributed  codewords,  C r(t)  and  C'{t  —  1),  and  the 
information  code  matrix  G'(t)  is 

C'(«)  =  C'(t -  l)G'(t),  G'(t)=[G(jt)  G°}].  (11) 
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Denote  h (t)  —  [hi(t)  whose  components  are  assumed  to  remain  almost  constant  for 

any  two  adjacent  codewords  but  are  assumed  to  be  independent  and  stationary  ergodic 
stochastic  processes  (i.e.,  h (t)  =  h (t  —  1)).  Then,  the  1x4  received  signal  row  vector  at 
the  destination  is  expressed  as 

y (t)  =  h(t)C'(t)  =  h(t)C(t  -  1  )G'(t)  =  h {t  -  l)C'(t  -  l)G'(t) 

=  [h(t  -  1)C'(«  -  1)]  G'(t)  =  y(f  -  l)G'(t).  •  (12) 

Therefore,  it  becomes  evident  that  the  differential  distributed  space-time  coding  schemes 
can  be  applied  to  a  cooperative  network  because  G'(f)  and,  equivalently,  G (t),  can  be 
recovered  from  y(t)  and  y (t  —  1),  without  knowing  the  CSI. 


VI.  Numerical  Results 

In  this  section,  we  show  some  numerical  results  to  support  our  findings.  The  first 
set  of  simulations  considers  a  set  of  2  x  2  cyclotomic  space-time  codes  as  the  prototype 
for  the  development  of  the  distributed  space-time  codes.  The  prototype  codeword  has 
256  constellations  over  two  time  slots.  For  comparison,  Alamouti’s  codes  with  the  same 
information  rate  are  also  considered.  All  the  channels  are  assumed  to  be  independent 
and  identically  distributed  (i.i.d.)  random  Gaussian  processes  with  zero  mean  and  unit 
variance.  The  relay  noise  is  not  considered  in  the  simulations. 

Figure  4(a)  compares  the  CER  performance  of  the  cooperative  diversity  schemes  using 
the  Alamouti’s  and  the  cyclotomic  codes.  The  SNR  is  defined  as  the  ratio  of  the  averaged 
codeword  energy  to  the  noise  energy  over  the  same  period.  It  is  seen  that,  the  cyclotomic 
codes  outperforms  the  Alamouti’s  by  about  1.5  dB. 

The  second  set  compares  different  detection  methods  for  the  Alamouti’s  codes,  where 
the  QPSK  modulation  is  used.  As  evidenced  in  Fig.  4(b),  the  expense  of  using  differential 
detection  is  the  3  dB  noise  amplification. 

VII.  Conclusion  and  Remarks 

We  have  proposed  a  novel  cooperative  protocol  and  several  distributed  space-time  codes 
for  application  to  wireless  networking  in  unmanned  surface  vehicles  (US Vs).  The  proposed 
protocol  is  applicable  to  both  amplify-and-forward  and  decode-and-forward  algorithms. 


i 
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(a)  Comparison  of  prototype  codes  (b)  Comparison  of  detection  methods 

Fig.  4.  CER  performance  of  different  distributed  space-time  schemes  and  detection  methods. 


The  two  proposed  distributed  space-time  coding  schemes  are  useful,  respectively,  when 
the  CSI  is  or  is  not  available.  The  advantage  of  the  proposed  methods  were  confirmed 
through  numerical  simulations. 

The  swarm  of  USV  architecture  has  similarities  to  the  sensor  network  and  ad  hoc  net¬ 
works.  While  our  contribution  in  this  paper  focused  on  the  effective  data  transmission, 
other  developments  in  several  fronts  of  these  areas  can  be  applied.  For  example,  routing 
algorithm  is  another  important  issue  to  be  considered  [29].  The  information  sensed  at  dif¬ 
ferent  nodes  are  usually  correlated  and  distributed  source  coding,  originated  from  Slepian 
and  Wolf’s  work  [30],  can  be  used  to  reduce  the  total  amount  of  data  without  compro¬ 
mising  the  information  [31],  [32].  The  importance  of  such  data  compression  is  particularly 
significant  when  large  amount  of  data,  e.g.,  acoustic  imaging  information,  is  to  be  trans¬ 
mitted.  However,  unlike  typical  sensor  network  application  where  energy  consumption  is 
of  great  concern,  it  is  of  secondary  priority  in  the  USV  networking  applications. 
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Abstract 

In  this  paper,  we  develop  a  protocol  for  the  construction  of  cooperative  networks  when  the  channel 
state  information  is  not  available  at  the  transmitters  and  the  receivers.  In  the  proposed  protocol,  differ¬ 
ential  space-time  codewords  are  generated  at  the  source  terminal.  In  the  broadcast  phase,  each  row  of  the 
differential  space-time  codeword  is  transmitted  to  a  different  relay,  whereas  in  the  relay  phase,  the  relaying 
terminals  retransmit  the  codeword  through  simple  amplify-and-forward  algorithm.  The  performance  of 
the  cooperative  diversity  system  is  analyzed  for  a  two-user  case  for  different  channel  environments  in  terms 
of  the  diversity  gain  and  the  diversity  product.  The  optimization  of  the  power  allocation  between  source 
and  relay  terminals  is  considered  for  the  maximization  of  the  diversity  product.  When  the  same  modu¬ 
lation  scheme  is  used,  the  performance  of  differential  detection  is  degraded  by  3  dB  noise  enhancement 
compared  with  coherent  detection. 
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I.  Introduction 


The  multi-input-multi-output  (MIMO)  technology,  incorporating  appropriate  space- 
time  coding  schemes,  allows  considerable  increase  of  the  information  capacity  in  fading 
channels  [1]— [4].  Without  increasing  premium  bandwidth  or  transmit  power,  a  MIMO 
system  can  achieve  higher  diversity  gains  and  data  rates,  compared  to  single-input  or/and 
single-output  system  counterparts. 

It  is  often  impractical,  however,  for  a  transmitter  to  host  a  large  number  of  antennas 
for  a  desirable  diversity  gain,  especially  for  small  size  transceivers  operating  at  relatively 
low  frequencies.  In  addition,  delivery  of  information  from  a  transmitter  to  a  receiver 
can  be  compromised  when  they  are  distant  apart,  or  when  the  wireless  channel  links  are 
highly  impairing.  This  difficulty  arises  in  a  wide  class  of  wireless  networks,  such  as  local 
area  networks  (LANs),  sensor  networks,  and  ad  hoc  networks.  For  these  applications, 
cooperative  diversity  exploiting  cooperation  among  multiple  terminals  proves  desirable. 

The  cooperative  diversity  techniques  have  recently  attracted  considerable  attentions.  In 
responding  to  the  increasing  needs  of  effective  and  reliable  wireless  networks  in  various 
applications,  the  development  of  cooperative  diversity  techniques  has  benefitted  from  the 
recent  advances  of  space-time  codes,  transmit  diversity,  and  MIMO  technologies.  Recent 
research  work  has  shown  system  feasibility  and  provided  capacity  analyses  of  different 
cooperative  diversity  techniques  (see  for  example,  [5]— [21])-  These  methods  assume  that 
channel  state  information  (CSI)  is  available  at  the  receivers,  although  few  assume  CSI 
knowledge  at  the  transmitters. 

CSI  knowledge  at  the  receivers  is  usually  obtained  through  channel  estimation,  either 
using  training  (pilot)  signals  or  utilizing  blind  methods.  However,  channel  estimation  is 
often  complicated  and  may  reduce  the  transmission  efficiency,  if  pilot  signals  are  used.  In 
addition,  the  estimation  becomes  unreliable  and  impractical  in  fast  fading  environments 
[22].  In  such  situations,  MIMO  systems  often  adopt  differential  space-time  coding  schemes 
which  use  differential  codes  between  two  adjacent  space-time  code  matrices  so  that  the 
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decoding  at  the  receiver  is  independent  of  the  fading  channels  [22]— [30]- 

As  the  cooperative  diversity  schemes  involve  both  broadcast  and  relay  phases,  multi¬ 
ple  cooperative  terminals  should  be  considered.  The  propagation  channels  are,  therefore, 
more  complicated  than  those  encountered  in  MIMO  scenarios  [31].  Accordingly,  the  con¬ 
sideration  of  system  design  without  assuming  knowledge  of  CSI  at  a  receiver  becomes  a 
practical  requirement. 

So  far,  several  differential  cooperative  diversity  schemes  have  been  developed  for  appli¬ 
cations  in  the  absence  of  CSI  knowledge  at  the  receivers.  A  repetition-based  differential 
amplify-and-forward  scheme  is  proposed  in  [32]  for  a  single-relay  scenario,  where  differen¬ 
tially  encoded  BPSK  signals  are  used.  In  [33],  a  differential  modulation  scheme  is  proposed 
for  two-user  cooperative  diversity  systems,  where  the  relay  terminals  uses  QPSK  modu¬ 
lation  to  transmit  two  BPSK  streams.  The  information  of  one  user  stays  in  the  in-phase 
axis  (I-axis)  while  that  of  the  other  user  stays  in  the  quadrature-phase  axis  (Q-axis).  In 
[34],  a  non-coherent  decode-and- forward  scheme  is  developed  for  binary  frequency  shift 
keying  (BFSK)  modulations. 

In  this  paper,  we  develop  a  new  distributed  space-time  modulation  scheme  for  coopera¬ 
tive  systems  that  require  no  knowledge  of  the  CSI  at  both  transmitters  and  receivers.  In 
the  proposed  protocol,  differential  space-time  codewords  are  generated  at  the  source  ter¬ 
minal,  and  each  row  of  a  codeword  is  relayed  by  a  different  relay  terminal.  Fundamentally, 
any  differential  space-time  codes  can  be  employed  in  the  proposed  scheme,  whereas  the 
relay  terminals  implement  a  simple  amplify-and-forward  relaying  algorithm.  An  entry  of 
the  user’s  differential  space-time  codeword  is  not  based  on  any  specific  modulation  such  as 
PSK  or  FSK.  The  performance  of  the  cooperative  diversity  system  is  analyzed  for  a  two- 
user  case  for  different  channel  environments  in  terms  of  the  diversity  gain  and  the  diversity 
product.  The  optimization  of  the  power  allocation  between  source  and  relay  terminals  are 
considered  for  the  maximization  of  the  diversity  product.  Roughly,  the  performance  of 
differential  detection  is  degraded  by  3  dB  noise  enhancement  compared  with  the  coherent 


23 


detection,  when  the  same  modulation  scheme  is  used. 

In  this  paper,  we  consider  memoryless  channels  and  assume  that  the  terminals  transmit¬ 
ting  the  same  distributed  space-time  codeword  are  synchronized  at  the  destination  receiver. 
When  the  channels  are  not  perfectly  synchronized,  the  synchronization  error  between  the 
terminals  can  be  considered  as  channel  dispersion  effect  [35].  The  effect  of  channel  dis¬ 
persion  in  cooperative  network  as  well  as  the  effective  approaches  using  equalization  and 
orthogonal  frequency  division  multiplexing  (OFDM)  techniques  have  been  considered  in 
[35],  [36],  [37], 

This  paper  is  organized  as  follows.  In  Section  II,  the  system  model  is  introduced.  In 
Section  III,  we  review  the  existing  cooperative  protocols,  and  then  propose  a  new  protocol 
which  is  effective  to  transmit  distributed  information  when  the  CSI  is  unavailable  at  the 
receivers.  The  differential  space-time  cooperative  structure  is  shown  in  Section  IV  and  the 
performance  analysis  is  provided  in  Section  V.  Numerical  results  are  provided  in  Section 
VI. 

II.  System  Model 

To  illustrate  the  concept  of  cooperative  diversity  in  a  wireless  network,  consider  a  simple 
narrowband  communication  model  as  depicted  in  Fig.  1(a).  Each  user  cooperates  with 
the  other  users  and  serves  as  a  relay  terminal  for  them.  Therefore,  each  user  receives 
an  attenuated  and  noisy  version  of  the  partners’  transmitted  signal  and  relays  it  to  the 
destination  or  other  relays.  The  destination  terminal  receives  a  noisy  version  of  the  sum 
of  the  attenuated  signals  from  all  users  (see  Fig.  1(b)).  Denote  the  number  of  users  as  M. 
The  received  signal  at  the  destination  during  the  Ith.  symbol  period  is  expressed  in  the 
baseband  model  as 

M 

y(l)  =  J2rYiohio(l)xi(l)  +n(l),  (1) 

i=l 

where  we  used  subscript  0  to  denote  the  destination  terminal  for  notational  convenience. 
In  (1),  Xi(l )  is  the  source  or  relaying  signal  transmitted  by  the  ith  user,  for  i  =  1, 2,  •  •  • ,  M, 
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7io  is  the  large-scale  attenuation  factor  of  the  channel  between  the  ith  terminal  and  the 
destination,  whereas  hi0(l)  represents  the  respective  unit- variance  small-scale  time- varying 
statistics  of  the  channel.  In  addition,  no(l)  is  the  additive  channel  noise  term  at  the 
destination. 

It  is  assumed  that  a  relay  terminal  does  not  receive  signals  from  other  users  when  it  is 
transmitting.  Therefore,  the  signal  received  at  relay  terminal  k  is  expressed  as 

r(k\l)  =  £  TStM0*i(0  +  n?°(0-  (2) 

ij^k 

where  7 &  and  hik{l),  respectively,  represent  the  large-scale  attenuation  factor  and  the 
small-scale  time- varying  statistics  of  the  channel  between  user  terminals  i  and  fc,  and 
(/)  is  the  additive  channel  noise  term  at  the  relay  receiver  k,  where  i,  k  =  1, 2,  •  •  • ,  M. 

Depending  on  how  the  relay  terminals  relay  the  signals  from  other  users,  two  major 
schemes  can  be  devised,  namely,  amplify-and-forward  and  decode-and-forward.  In  this 
paper,  we  focus  on  the  simpler  amplify-and-forward  scheme  which  requires  the  least  com¬ 
plexity  at  the  relay  terminals. 

III.  Cooperation  Protocols 

A.  Existing  Protocols 

Several  cooperative  diversity  protocols  have  been  investigated.  The  early  protocols  used 
repetition-based  approaches,  which  is  illustrated  in  Fig.  2(a)  [8],  [19].  After  source  terminal 
i  broadcasts  a  block  of  information  symbols  (e.g.,  Sj(£)  =  [si(l),  Si(l  +  1),  •  •  • ,  Sj(/  +  L  —  1)], 
where  L  is  the  length  of  the  data  block)  to  the  destination  and  the  relay  terminals,  all 
relay  terminals  repeat  the  same  data  (i.e.,  Si(l )  =  [si(l),Si(l  +  1),- •• ,  Si(l  +  L  —  1)])  in  a 
sequential  order. 

Recently,  more  effective  protocols  have  been  developed  to  take  advantages  of  the  ad¬ 
vances  of  MIMO  space-time  codes.  For  example,  [8]  and  [19]  proposed  space-time  coop¬ 
eration  protocols  which  are  illustrated  in  Fig.  2(b)  and  (c).  In  the  broadcast  phase,  the 
source  terminal  i  broadcasts  a  block  of  information  symbols  (e.g.,  s i(l)  =  [sj(/),  Si(l  + 
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1),  •  •  • ,  Si(l  +  L  —  1)])  to  the  destination  and  the  relay  terminals.  In  the  relay  phase,  for 
scheme  I  depicted  in  Fig.  2(b),  the  source  and  the  relay  terminals  transmit  different  rows 
of  a  space-time  codeword  Xi( u),  which  is  constructed  from  s *(/),  where  v  denote  the  index 
of  a  coded  data  block.  Space-time  cooperative  scheme  II  is  very  similar  to  scheme  I.  The 
only  difference  between  them  lies  in  the  fact  that  in  scheme  I,  the  source  terminal  takes 
part  in  the  relay  phase,  whereas  in  scheme  II  it  does  not.  Such  protocols  provide  more 
effective  bandwidth  use,  particularly  when  the  number  of  cooperative  users  is  large. 

B.  Proposed  Protocol 

Most  existing  cooperative  protocols  assume  that  the  CSI  is  known  at  the  receivers.  In 
this  section,  we  develop  a  protocol  that  is  applicable  to  space-time  cooperation,  where  the 
CSI  is  unavailable  at  the  receivers.  The  proposed  protocol  performs  differential  distrib¬ 
uted  space-time  coding  and  information  detection  using  the  relatively  simple  amplify-and- 
forward  algorithm.  Fig.  3  depicts  the  proposed  protocol.  When  considering  user  i  as  the 
source  user,  the  other  M  —  1  ones  act  as  its  relays.  At  the  source  terminal,  an  M  x  M 
space-time  codewords,  Xj(^),  is  formed  from  the  source  information  symbols.  Each  of  the 
M  —  1  relay  terminals  receives  a  different  row  of  the  codeword  Xj(i/)  that  is  transmitted 
from  the  source  user  through  the  broadcast  phase  (i.e.,  during  the  first  M  —  1  time  blocks). 
All  the  M  rows  of  Xj(i/)  are  transmitted  from  the  M  terminals  in  the  relay  phase  which, 
spans  the  Mth  time  block.  It  is  noted  that  the  destination  listens  to  the  source  during  the 
broadcast  phase,  whereas  it  listens  to  the  source  and  the  relay  terminals  during  the  relay 
phase. 

Consider  the  simple  case  of  two  users,  i.e.,  M  =  2.  The  source  terminal  forms  a 
2x2  MIMO  space-time  codeword  X,(i^)  from  the  source  information  symbols  s i(l)  = 
[si(l),Si(l  +  1),  •  •  • , Si(l  +  L  —  1)].  Then,  using  the  cooperative  protocol,  the  source  and 
relay  terminals  equivalently  transmit  the  following  2x4  codeword  X((j/)  in  the  cooperative 
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diversity  system,  i.e., 

\Xi(utl)  Xi(i/,  2) 

XAu)  = 

.Xi(u,3)  Xi(v,  4) 


3)  Xi{y,  4)  Xi(is,  1)  2) 

0  0  Xi(u,  3)  Xi(v,  4) 


where  r)  is  the  complex  scaled  and  noisy  version  of  Xi(u,  r).  In  the  above  space-time 
codewords,  the  columns  specify  the  time  dimension,  whereas  the  rows  specify  the  space 


dimension  (i.e.,  antennas  or  relay  terminals). 


IV.  Differential  Cooperation  Schemes 

In  this  section,  we  first  briefly  review  the  differential  space-time  codes  developed  for 
MIMO  systems,  and  then  show  the  feasibility  of  employing  them  in  the  space-time  coop¬ 
eration  applications  through  the  proposed  protocol. 


A.  Differential  Space-Time  Codes 

In  order  to  deal  with  rapid  MIMO  fading  environments  where  the  CSI  is  unavailable 
at  the  receiver,  differential  space-time  codes  can  be  used.  For  essence,  differential  coding 
between  two  adjacent  space-time  code  matrices  is  used  so  that  the  decoding  at  the  receiver 
can  be  performed  without  the  knowledge  of  channels  [22]— [30].  Differential  space-time 
codes  are  most  commonly  used  in  the  form  of  unitary  matrices  [22],  [27]. 

Let  Q  —  {Gj,  i  =  1,  *  •  • ,  N}  be  a  set  of  M  x  M  unitary  matrices  (i.e.,  G,Gf  =  GfG  *  = 
IM,  where  I m  is  the  M  x  M  identity  matrix),  and  C(0)  is  an  M  x  M  matrix.  log2  N 
information  bits  are  mapped  to  one  of  the  N  unitary  matrices  Gj  in  Q.  To  send  the 
message  Gfy)  e  Q  over  time  block  the  transmitter  sends  C (u)  where 

CM  =  C(i/  -  1)GM  =  C(0)  n  G(t).  (4) 

T— 1 

Matrix  Cfy)  takes  the  following  general  format 

cu(y)  cnW)  -  ciAf(i/) 

C2i(v)  c22(o)  ...  c2M(v) 

cmiH  cM2(u)  ...  cMM(v ) 


27 


The  spectral  efficiency  of  this  code  is  (1/M)  log 2N  bits  per  channel  use.  When  the  channels 
corresponding  to  different  transmit  antennas  are  independent  and  identically  distributed 
(i.i.d.),  C(z/)  is  often  designed  to  take  the  form  of  unitary  matrices  for  improved  system 
capacity  [23],  [24],  [26].  That  is, 

CMCV)  =  I„  (6) 

In  this  case,  C(0)  is  chosen  as  unitary. 

The  differential  space-time  coding  scheme  works  well  in  time- varying  channels,  provided 
that  the  channels  do  not  change  significantly  over  a  period  of  two  adjacent  codewords,  i.e., 
h(i/— 1)  =  h(z/),  where  h(u)  =  [hi o(r'),  •  •  • ,  h,Mo(v)]  is  a  vector  channel  with  the  ith  element 
representing  the  channel  coefficient  between  the  ith  transmit  antenna  and  the  destination 
receiver.  Denote  the  signals  received  at  the  receiver  corresponding  to  the  (//  —  l)th  and 
uth  codewords  as 

y(u  -  1)  =  h {v  -  l)C(v  -  1)  +  n [y  -  1) 

and 

y(u)  =  h(i/)C(i/)  +  n(z/), 

then  we  have 

y{u)  =  h(z/  -  l)C(i/  -  l)G(u)  +  n(u)  =  [y(v  -  1)  -  n(u  -  l)]G(i/)  +  n{v).  (7) 

Because  G(z^)  is  unitary,  —n(u  —  l)G(i/)  can  be  considered  as  additional  noise  term  with 
the  same  variance  as  n(i/).  Therefore,  the  information  G(u)  can  be  recovered  from  y(u—  1) 
and  y(z')  through  the  above  relationship  at  the  expense  that  the  equivalent  noise  power  is 
doubled. 

B.  Differential  Space-Time  Cooperative  Schemes 

We  now  consider  the  feasibility  of  extending  the  differential  space-time  codes  to  cooper¬ 
ative  networks.  For  the  clarity  of  the  presentation  of  the  proposed  method,  we  ignore  the 
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noise  effect  in  this  section.  The  effect  of  noise,  and  the  optimization  of  power  allocation 
in  different  antennas  as  well  as  different  broadcast  and  relay  phases,  will  be  considered  in 
Section  V  for  a  two-user  scenario. 

Assume  that  there  are  M  users  cooperating  with  each  other.  Without  loss  of  generality, 
we  only  focus  on  the  first  user’s  information  and  the  user  index  is  omitted  without  confu¬ 
sion.  Following  the  protocol  depicted  in  Section  III,  C'  £  CMxM2  js  constructed  from  the 
corresponding  differential  space-time  codeword  C  £  Q,  expressed  as 

Ci\{v)  •••  c2M(y)  •••  CmiW)  •••  Cmm{v )  Cn(u)  •••  clM(v) 

0  •••  0  C2i(l/)  •••  C2M(v) 

0  .  0  cm\(v)  •••  cmm(^) 

(8) 

It  is  assumed  that  signals  transmitted  from  different  terminals  are  synchronized  at  the 
destination  receiver. 

Because  different  rows  of  the  DSTC  codeword  undergo  different  propagation  channels 
with  different  attenuations,  it  is  desirable  to  allocate  energy  adaptively  to  different  rows. 

For  this  purpose,  (8)  is  generalized  to 

ctiVKonty)  •••  aM\/P[cMMiy)  a\>/Kc\i(v)  ■■■  aiVPiCiM(v) 

0  •••  0  VKc2l(v)  •••  \fP2C2Miy) 

0  •••  0  \J  P V  P MCMA'l(y) 

where  >  0  is  the  transmit  power  of  different  users,  i  —  1,  -  •  • ,  M,  and  at  >  0  is  used 
to  adjust  the  transmit  power  at  the  source  terminal  corresponding  to  different  codeword 

M  2 

rows.  Because  C {u)  is  unitary,  J2  \cik\  =  1  for  i  =  1,-  •  •  ,M.  Therefore,  the  total 

fc=i 

M  M 

codeword  power  of  the  DSTC  in  (9)  is  P  =  Pi  Y,  a?  +  Y,  Pi-  The  optimization  of  the 

i= 1  i— 2 

power  allocation  is  considered  in  Section  V. 

Next,  we  consider  the  detection  of  differential  DSTC  codeword  C'(^)  without  using  the 
CSI.  From  (4)  and  (9),  the  relationship  between  two  adjacent  differential  DSTC  codewords, 
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C'(u)  and  C'(u  —  1),  and  the  information  codeword  G(v)  is  expressed  as 

C'H  =  CV  -  l)G'(iz)  =  C'(0)  n  G'(r),  (10) 

r=  1 

where 

G'(u)  =  diag  [  G(v),  -  •  • ,  G(i/)  ]  =  Iu  ®  G (*/),  (11) 

with  0  denoting  the  Kronecker  product  operator.  Because  G(^)  is  unitary,  it  is  clear  that 
the  M2  x  M2  matrix  G'(z')  is  also  unitary  (i.e.,  G'(u)(G'(i'))H  =  1m2)- 
Note  that  the  encoding  process  of  (10)  is  performed  at  the  source  terminal  and,  therefore, 
it  does  not  require  the  relay  terminals  to  perform  encoding  or  to  possess  information 
required  in  the  differential  coding  process.  In  addition,  it  is  noted  that  the  distributed 
space-time  codeword  C'(z')  is,  in  general,  not  unitary.  However,  because  G'{y)  is  unitary, 
we  have 

C»[C'(0]H  =  C'(i/-l)G'(t')[G»]"[C>-l)]''  =  C'(i/-l)[C'(i/-l)]i'  =  C'(0)[C'(0)]H 

(12) 

Therefore,  the  transmit  power  of  G(v)  is  unchanged  for  different  v  and,  subsequently,  the 
robustness  of  the  differential  code  is  guaranteed. 

Denote  h(i')  =  [ hio{v ),  hi2(^)/i2o(i/),  •  •  • ,  hi,M(v)hMfi(v)\  as  the  lx  M  channel  vector 
whose  components  are  assumed  constant  for  any  two  adjacent  codewords  (i.e.,  h(z/)  = 
h(i/  —  1))  and  are  independent  and  stochastic  processes  over  time.  Then,  the  1  x  M 2 
received  signal  vector  at  the  destination  terminal  is  given  by 

y(i /)  =  h(v)G(v)  =  h(i/) C'(u  -  l)G'(u) 

=  h(i/-l)C;(z/-l)G'(i/)  (13) 

=  y{v-l)G'(v). 

From  (13),  it  becomes  evident  that  the  differential  detection  is  feasible  in  cooperative 
networks  because  G'(u)  and,  equivalently,  G(z/),  can  be  recovered  from  y (u)  and  y (u—  1), 
without  the  knowledge  of  the  channel  characteristics. 
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V.  Performance  Analysis 


A.  Signal  Model 

In  this  section,  we  consider  the  performance  of  the  proposed  system.  For  simplicity,  we 
only  consider  a  two-user  scenario  (ref.  Fig.  4),  while  the  extension  to  multi-user  scenarios 
is  straightforward.  Because  of  the  symmetry,  we  focus  only  on  the  transmission  of  user  l’s 
information,  whereas  the  second  user  is  regarded  as  the  relay. 

Consider  the  generalized  codeword  C'(^)  which  considers  the  relay  noise  and  takes  the 


following  form 


C»  = 


VKc»i(v)  \TPiC22{v)  OL\ft\cn{v)  cc/PiAaO')' 

0  0  \/p>C2 i(u)  sJP2C22{v)  J  ’ 


where  denotes  the  attenuated  and  noisy  version  of  c^,  i,k  =  1,2.  The  expression 
of  Cik  will  be  implicitly  given  later  in  (17).  The  total  energy  used  to  transmit  C'(z')  is 
P  =  (l+»2)P[-f  Note  that,  while  we  only  use  a  single  a  in  the  first  row  (compared  to  a\ 
and  a.2  in  (9)),  it  is  obvious  that  a,  Pi,  and  P2  together  provide  enough  degrees-of-freedom 
to  adjust  the  power  over  different  terminals  and  different  time  slots. 

We  first  consider  the  signal  flow  during  time  block  v  which  spans  four  time  slots.  In  the 
first  two  time  slots,  c2i(^)  and  022(1')  are  transmitted  from  user  l’s  antenna.  As  the  result, 
the  signals  received  at  relay  terminal  (user  2)  and  the  destination  receiver  are  expressed, 


respectively,  as 


[i’i(i'),r2(i/)]  =  sJlPiiuhwiy)  [ c2i(z/) , c22(i')  ]  +  [nri(u),nr2(^)]  (15) 

ll/i(i/),3/2(i/)]  =  yfPiliohio(v)  [C2i{iy),c22^)}  +  [ni{v) ,n2{v)} ,  (16) 

where  Ui{u)  and  nri(u )  denote,  respectively,  the  receiver  noise  at  the  destination  and  the 
relay  at  time  slot  i  within  time  block  v. 

At  the  second  time  block,  Cn( v)  and  012(1')  are  transmitted  from  user  l’s  antenna,  and 
ri(v)  and  r2(i')  are  relayed  from  user  2’s  antenna  after  the  amplification.  The  received 
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signal  at  the  destination  receiver  becomes 

r  _  _  ,  [ acn(u)  aci2(^)l 

[lfc(*,)>!fe(*')]  =  VPiliohioiv),  V^72oh20(^)  „  A  ,,  +[n3(v),n4(v)] 

L  c2\{y)  c22{u)  J 


—  [\/A7lO  hio{v) ,  V^A7207l2^2o(l/)^12(f/)fil2 


ctcu(v)  acn(u ) 


If  J-  i  \  /  >  V  ^  V  ^  "  1  /V  yrx 

L  C2l(u)  C22M 

+\fKl2oh2o{v)g2[nTl{v),nr2(v)]  +  [n3(i/),n4(u)] 
r—  ,  ,  r—  ,  s  ,  .  i  \®cn(y)  acn{u)] 


'/Pi'Yio  /iio(^),  Vp272oh2o(v)hn{v)g2 


C21M  c22(^) 


I P2  720 , 


+VA^20^[nrl(l/),nr2(l/)1  +  > 


where 


<?f  =  (£  [MAI2  +  MAI2]}-1  =  {E  [A7LIMAI2  +  n*rl(u)  +  nj2(u))y 


[PiTu  +  2crrn  —  A  A122  1  +  p"^ 


is  a  power  normalization  factor  such  that  the  total  energy  transmitted  from  the  relay 
terminal  over  the  two  symbols  is  P2,  and 

(19) 

is  the  power  ratio  at  which  the  relay  terminal  transmits  desired  signal  in  the  presence  of 
relay  noise.  In  (18)  and  (19),  cr2n  —  E  [|nrl(A[2]  =  E  [|nr2(A|2]  is  the  variance  of  the  relay 
noise. 

Equations  (15)  -  (18)  can  be  combined  into  the  following  format 

y(A  =  h(t/)C'(A  +  n(A  +  \h^-—h20g2nr(i'),  (20) 

V  t'- 1  712 


where 


c'(A  = 


y(A  =  MA^A^A^AI 
h(A  =  [7loAo(A;72oAo(AA2(A<?2] 

n(u)  =[ni(v),n2(u),n3{i'),n4{i')] 
nr(u)  =  [0, 0,  nri(i/),  nr2{u )], 

\[P\C2  i(A  VPiC22(v)  aJPlcn{v)  asfP[cX2{v) 
0  0  sfP2c2X{v)  \fP2C22tv) 
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It  is  noted  that,  because  of  the  relay  noise,  the  average  power  that  the  relay  terminal  uses 
to  transmit  the  useful  information  is  P^gl/Z  for  each  symbol,  whereas  the  average  power 


transmitted  from  the  relay  terminal  is  P2/2  per  symbol. 

Similar  to  (20),  for  time  block  v  —  1,  we  have 

I p2  720 
Pl7l2 
I P 2  720 


y(y  1)  =  h(u  -  1)0 (v  -  1)  +  n(i/  -  1)  + 

=  h(u)C'(u  —  1)  +  n(u  -  1)  + 


Pi  712 


hohnr{^  ~  1) 
h2ohnr{v  -  1), 


(22) 


where  again  h(i/)  —  h (v  —  1)  is  assumed.  With  the  noise  components  considered,  y (u) 


and  y(v  —  1)  are  now  related  through  the  following  relationship 


y(u)  =  y{v  -  l)G’(u)  +  n (1/), 


(23) 


where 


n{u)  = 


n(i')  +  \PE-~h2ohnr(v) 
V  a  712 


n(i/  -  1)  +  4 /5— /i2o52nr(^  -  1) 
V  r 1  712 


G'W.  (24) 


Because  G’{v)  is  unitary,  it  is  evident  from  the  above  equation  that  the  differential  de¬ 
coding  doubles  the  noise  power  or,  equivalently,  reduces  the  SNR  by  3  dB. 

B.  Maximum  Likelihood  Detection 

The  covariance  matrix  of  the  combined  noise  vector  h(u)  is  obtained  as 


R.n  =  E[  nH(v)n(v)]  =  2  al 


T  .  P‘2  720  arn  -2r> 

14  T\iD<92 


= 


(25) 


where  <rj  =  E  [|m(i/)|2]  =  E  [|n2(  i/)|2j  is  the  variance  of  the  destination  receiver  noise, 
D2  =  diag[0,  0,  1,  1],  and 

Pillar 


A2  =  2 


I4  + 


Pi  7l2  o* 


#2D2 


=  2 


I4  + 


2  „2 

rn  Do 


(Pah  +  2 cr?n) 


2  [I4  +  W  D2]  (26) 


is  a  diagonal  matrix  implicitly  defined  in  the  above  equation,  where 


W 


,2  x-r2 

rn 


vl(Pah  +  2VrnY 


(27) 


33 


The  maximum  likelihood  detection  results  in  the  estimation  of  G(u)  as 

G(u)  =  argmintr  j[y(^)A-1  —  y(u  —  1)(I2  ®  G)A_1]  [y(i/)A_1  -y (v  -  1)(I2  0  G)A_1]H} 

=  argmintr  { [y(u)  -  y{v  -  1)(I2  ®  G)]  A-2  [y(^)  -  y{v  -  1)(I2  ®  G)]H}  , 

G  L  J  (28) 

where  “tr”  denotes  the  trace  of  a  matrix. 

C.  Performance  Analysis 

While  the  accurate  performance  of  a  differentially  coded  system  requires  the  considera¬ 
tion  of  the  quadratic  receiving  structure  [22],  [23],  [26],  it  can  be  well  approximated  in  high 
SNR  situations  by  using  an  equivalent  coherent  receiver  model  (23)  with  known  channel 
vector  y(u  —  1)  and  enhanced  noise  power  of  n (is)  [25]. 

From  (25),  it  is  clear  that  the  variance  of  the  combined  noise  term  during  the  last  two 
symbols  is  larger  than  that  during  the  first  two  symbols.  Assume  that  the  covariance 
matrix  R„  =  cr2  A2  can  be  estimated  at  the  destination  receiver.  Then,  right  multiplying 
A-1  to  (23)  results  in 

y(i/)A-1  =  y(u  -  l)G/(i/)A-1  +  n(i/)A_1  (29) 

or 

y'{y)  =  y(v  -  i)G(* ')  +  n'{u),  (30) 

where  y'(^)  =  y(i/)A"1,  G(i/)  =  G'(^)A-1,  and  n'(z/)  —  ii(i/)A_1  has  a  covariance  matrix 
cr2I4.  Multiplying  y(u)  by  A-1  amounts  to  adjusting  the  signal  strength  at  the  receiver 
so  that  the  noise  power  of  the  four  symbols  becomes  equal.  It  does  not  affect  the  total 
power  transmitted  at  the  source  and  relay  terminals. 

From  (29)  and  (30),  it  is  seen  that  the  detection  problem  becomes  equivalent  to  finding 
unknown  unitary  codeword  G'(^)  in  known  equivalent  channel  vector  y(u—  1),  known  noise 
covariance  matrix  <r2A2,  and  independent  and  temporally  white  noise  components.  From 
the  above  two  equations,  we  can  derive  the  pairwise  codeword  error  probability  (CER), 
i.e.,  the  probability  of  transmitting  G'  —  diag[G,  G]  and  deciding  in  favor  of  another 
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E'  =  diag[E,  E]  at  the  detector,  conditioned  by  the  equivalent  channel  vector  y (u  —  1). 
The  CER  is  given  by 

P(G'  -  E'|y(i/  -  1))  =  Q  (\AP( G,  E)/(2ffJ))  ,  (31) 

where  d2( G,  E)  is  the  distance  between  the  received  signals  corresponding  to  respective 
codewords  G(i/)  and  E(z/)  and  is  expressed  as 


<i2(G,E)=y(I/-l)Ay"(J/-l), 


(32) 


A  =  (G(P)  -  Efi-j)  (G(i/)  -  E(i/j)  "  =  (G'lv)  -  E'(i/))  A“2  Ki'(r)  -  E »)" .  (33) 

In  calculating  the  distance  of  (32),  we  make  the  following  approximation  y(z/)  fa  h(i/)C'(v) 
for  moderate  or  high  SNR  scenarios.  Then,  (32)  becomes 

<Z2(  G,E)  =h(,-l)CV-l)A(CV-l))‘h>-l) 

=  -  1)A  (C’(v  -  1 ))"  £„h"(i/), 

where 

Sh  =  diag[(Ti,  a2]  =  diagpyio,  72052],  (35) 


and 

h(z/)  =  h(i/)Eh1/2  =  [/iio(i'),  hwiy)hto(v)]  (36) 

is  the  normalized  channel  vector.  For  simplicity  of  the  performance  analysis,  we  assume  the 
relay  channel  hniy)  is  dominated  by  time-invariant  component  (e.g.,  the  source  terminal 
and  the  relay  terminal  have  a  clear  line-of-sight  and  the  scattering  is  relatively  weak)  and, 
therefore,  can  be  approximated  by  a  unit  value.  In  addition,  hi0(v)  and  h2o(E)  are  assumed 
to  be  uncorrelated  complex  Gaussian  processes.  In  this  case,  h(z')  =  [hio(u),  /i2o(^)|  ~ 


eV(o.i2). 

In  the  underlying  situation,  the  average  power  per  symbol  is  Es  —  P/ 4  (i.e.,  the  total 
energy  of  a  DSTC  codeword  is  P).  To  clearly  show  the  relationship  between  the  pairwise 
CER  and  the  input  SNR,  we  rewrite  (31)  as 


P(G'  -  Yj\y[y  -  1))  =  Q 


/ 

<P(G,E)  2 Es) 

r-oP(G,E  )ES] 

u 

P  a2 

<  exp 

L  p  <\ 

(37) 
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where  the  last  inequality  is  the  Chernoff  bound  [2].  Because  y(u)  is  approximated  as  a 
linear  combination  of  i.i.d.  channels  h(n)  and,  therefore,  constitute  a  set  of  dependent 
channel  coefficients,  averaging  the  above  bound  with  respect  to  y(u  —  1)  results  in  [2] 

P(G'  E')  <  ft  (l  +  ,  (38) 

where  d\  and  d2  are  the  eigenvalues  of  the  following  2x2  matrix 


K  =  £hC'(v  -  1)A  (C\u  -  1  )f  Sh. 


In  the  following,  we  consider  two  different  cases.  In  the  first  case,  matrix  K  is  full  rank, 
whereas  in  the  second  case,  matrix  K  is  rank  one. 

Case  I:  Full  Rank  Channel  Environment 
When  Esdi/(P<Jn)  >  1  for  i  =  1, 2,  the  CER  can  be  approximated  as 


p(G'->Eo<n(i+§^ 

i=  1  V  °nr 


-1  ( E  \-2  2  / d  \~l 

<e  I  s  I  TT  I  Ut  \ 

-  \  ^2  )  11  {  ~p  ) 
\°nj  i=l  \r  / 


es\~2  *  (4diyl 

4 MU; 


The  term  (j^k)  determines  the  diversity  gain  which  indicates  the  slope  at  which  the 
BER  changes  with  the  input  SNR,  whereas  the  term  flLi  Up-)  ’s  a  constant  gain  in 
controlling  the  BER  and  is  referred  to  as  the  diversity  product.  In  this  case,  the  diversity 
gain  of  the  system  is  two.  Using  the  above  results,  it  is  seen  that  the  diversity  product  is 
determined  by  the  minimum  product  of  d\d2l  i.e.,  the  minimum  determinant  of  K.  It  can 
be  shown  that  (refer  to  the  Appendix) 

mindet(K)  =  —  |det(Sh)det(C2)|2  (min  |det(G  -  E)|)2 

+  P4ut°l (minCl(G  ”  E)(G  “  E)"c?)2 

(41) 

+~i(iTwy  (minc>(G  -  E)(G  -  E)H°?)2  - 


where  W  is  defined  in  (27),  and 


Cl  =  [c2i  i  C22]  1 
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c2  = 


aty/Pi  0 

Cn 

C12' 

~cxy/P[cn 

us/Picn 

.  0  y[P2. 

.C21 

C22. 

.  y/PzC2i 

V^C22  . 

(43) 


The  value  of  cL(G  —  E)(G  —  E)Hcf  is  lower  bounded  by  Amin[(G  —  E)(G  —  E)H],  which 
is  the  minimum  value  of  the  smallest  eigenvalues  with  respect  to  all  the  combinations  of 
G  and  E.  That  is, 

ci(G  -  E)(G  -  E)Hcf  >  Amin[(G  -  E)(G  -  E)H].  (44) 

Define 


P  = 


Amin[(G-E)(G-E) 


hi 


(45) 


min  |det(G  —  E)| 

as  a  factor  determined  by  the  used  differential  code.  Then,  the  diversity  product  is  lower 
bounded  by 

mindet(K)  (  «  )’  >  (min  |det  [G  -  E]|)2 .  (46) 

Now  we  consider  the  maximization  of  (46)  in  terms  of  the  power  allocations,  leading  to 
the  optimal  selection  of  Px,  P2,  and  a.  Because  (min  |det  [G  -  E] | ) 2  is  independent  of  the 
power  allocation,  it  suffices  to,  consider  only  the  following  part  of  (46)  for  the  purpose  of 
optimization  of  the  power  allocation, 

4Oa7?0720»!  [a2  +  (1  +  HO/?2] 


•  [(1  +  a2)  P,  +  PJ2  (1  +  W)2 

By  letting  dQi/da  =  0,  we  get 


(47) 


P2J20P 


2  rr-2 

rn 


P2+ IT 


+ 


(48) 


a2pt  [l  2(1  +  W)02  +  P2/Pl\  |l  2(l  +  +  ^  r  ■  A 

where  [x]+  =  max(0,  x)  (This  operator  is  used  to  ensure  that  a2  takes  non-negative  values). 
While  the  optimum  power  allocation  between  Px  and  P2  can  be  obtained  from  the  above 
expression,  the  analytic  result  is  prohibitively  complicated.  However,  the  above  results 
are  helpful  in  determining  the  optimum  power  allocation.  Specifically,  when  the  interuser 
channel  between  the  source  and  relay  terminals  has  a  good  quality,  i.e.,  7^  is  large  and 
o2n  is  small,  then  g2  &  1  and  W  «  0,  and  Qi  becomes 

4PiP27i2o722o 


Qi 


(1  +  a2)Pi  +  P2 


(a2  +  P2). 


(49) 
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(50) 


(51) 


Letting  dQ[/dPi  —  0  yields 

(1  +a2)Pl  =  P2. 

Substituting  this  result  back  to  (49) ,  we  obtain 

n,  ...  7io720  (°2  +  P2) 

In  particular,  when  01  =  1,  Q'l  opt  is  independent  of  a,  providing  that  (50)  is  satisfied. 
On  the  other  hand,  a2  assumes  a  large  value  when  j32  <  1,  resulting  in  Pi  P2  and 
a2P\  «  P2. 

Therefore,  the  optimum  value  of  a  depends  on  /3.  That  is,  given  the  same  environment, 
different  differential  codes  do  not  necessarily  allocate  the  power  in  the  same  manner.  When 
the  relay  noise  is  not  negligible,  the  optimal  value  of  Q[  is  compromised  by  the  relay  noise 
and  in  general  /3-dependent.  Therefore,  a  differential  code  with  higher  min  |det(G  —  E)| 
does  not  necessarily  imply  a  high  diversity  product  in  the  underlying  cooperative  diversity 
applications. 

It  is  clear  that,  to  optimize  the  power  allocation,  the  channel  attenuation  characteristics 
have  to  be  known.  While  in  this  paper  we  assume  that  the  exact  CSI  is  unavailable  at  the 
transmitters  and  receivers,  we  maintain  that  it  is  often  practical  to  assume  the  channel 
quality  at  both  transmitters  and  receivers  because  the  attenuation  characteristics  vary  with 
time  in  a  much  slower  manner.  It  is  emphasized  that  such  large-scale  channel  attenuation 
characteristics  are  not  required  in  employing  the  proposed  differential  distributed  space- 
time  coding  scheme,  but  a  good  estimation  of  such  information  can  be  used  to  improve 
the  system  performance. 


Case  II:  Rank  1  Channel  Environment 

When  Esdi/(Pa2)  1  whereas  P3d2/(Po-2)  <C  1,  (40)  is  no  longer  applicable.  It 
happens,  for  example,  when  either  the  direct  link  or  the  relay  link  is  highly  impairing.  In 
this  case,  (38)  becomes 


When  the  relay  channel  is  not  reliable,  that  is,  either  720  or  712  is  very  small,  or  the 
relay  noise  a^n  is  very  large,  it  is  obvious  that  the  relay  should  not  be  used  and  the  system 
degenerates  to  the  non-cooperative  system. 

On  the  other  hand,  when  the  direct  channel  between  the  source  and  the  destination  is 
unreliable,  i.e.,  710  «  0,  the  system  becomes  a  pure  relay  structure.  In  this  case,  a  should 
be  chosen  as  0,  and  K  defined  in  (39)  becomes 

P^~9l 


K  =  diag 


0, 


Ci(G  -E)(G-E)Hcf 


2(1  +  W) 


(53) 


Using  the  results  of  (44)  and  (45),  the  minimum  value  of  Adi/P  in  (52)  is  obtained  as 


mm 
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min  [Cl  (G  -  E)  (G  -  E)H  cf 
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r/3  min  jdet(G  —  E)|. 


Pi  +  P2  [al{Plii2  +  2<7*J  +  P2lWrn\ ' 

Because  j3  min  |det(G  —  E)|  is  independent  of  the  power  allocation,  we  only  consider  the 


following  term  for  the  optimization  of  the  power  allocation, 


Qi 


2PiP27?2722o 


o: 


2 

rn 


Pi  +  P2  [Vn(Ppfh  +  2<j2J  +  P2 720 cr. 
Then,  dQ2/dPi  =  0  results  in 


P> 


+  llolhPl  -  OnOrn  CTn/  (l^rn) 


(55) 


(56) 


In  particular,  when  we  assume  Pi 712  >  2<j^n,  that  is,  the  interuser  channel  has  a  high 
quality  and  the  input  SNR  at  the  relay  terminal  is  high,  the  above  result  becomes 


P2T20/W  =  Piln/vm-  (57) 

Therefore,  the  power  at  different  relay  portions  should  be  assigned  proportional  to  the 
propagation  gain  and  inversely  proportional  to  the  standard  deviation  of  the  noise  gener¬ 
ated  in  each  relay  portion. 

VI.  Numerical  Results 

In  this  section,  we  consider  a  simple  situation  where  two  users  cooperate  to  transmit  their 
respective  information  to  the  destination,  as  illustrated  in  Fig.  4.  Two  differential  DSTC 
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codes  are  considered.  In  Scheme  1,  the  differential  DSTC  code  is  constructed  based  on 
Alamouti’s  space-time  code,  where  each  symbol  uses  a  QPSK  signal  constellation,  results 
in  N  =  16  codewords  in  the  two  time-slot  space-time  code.  In  Scheme  2,  the  differential 
DSTC  is  constructed  based  on  the  differential  MIMO  space-time  code  developed  by  Liang 
and  Xia  [27].  For  the  latter,  the  constellation  size  remains  N  —  16.  Unless  otherwise 
specified,  the  power  is  optimally  allocated  by  adjusting  Pi,  P2,  and  a  to  maximize  Q\  in 
(47),  provided  that  the  average  power  per  symbol  is  unity. 

A.  Comparison  with  Co-located  MIMO  Systems 

In  the  first  set  of  simulations,  we  compare  the  CER  performance  of  the  proposed  co¬ 
operative  network  with  a  MIMO  system  consisting  of  two  co-located  antennas.  The  two 
antennas  in  the  MIMO  system  are  assumed  to  have  independent  Rayleigh  channels  to  the 
destination  with  unit  variance.  For  the  cooperative  system,  we  assume  710  =  720  =  1  and 
712  =  \/l25,  which  corresponds  to  the  situation  that  the  interuser  distance  between  the 
source  and  relay  terminals  is  one-fifth  of  the  distance  between  the  source  terminal  and 
the  destination,  and  the  power  attenuation  exponential  is  3.  In  this  case,  the  channels  to 
the  destination  terminal,  i.e.,  h[0(o)  and  h2o{u),  are  assumed  Rayleigh  faded,  whereas  the 
interuser  channel  is  considered  non-fading,  i.e.,  h\2{v)  =  1. 

Fig.  5  compares  the  performance  of  the  cooperative  diversity  system  and  the  correspond¬ 
ing  MIMO  system.  The  SNR  is  defined  as  the  reciprocal  of  the  receiver  noise  power,  as 
the  transmit  power  is  normalized  to  unity  per  symbol  period.  The  relay  receiver  and  the 
destination  terminal  receiver  have  the  same  noise  power.  It  is  seen  that,  the  cooperative 
diversity  system  achieves  the  same  diversity  gain,  and  there  is  a  performance  degradation 
due  to  the  relay  noise  in  the  cooperative  diversity  system.  The  performance  degradation 
is  very  small  when  the  channel  quality  between  the  source  and  relay  is  reasonably  good. 

We  notice  that,  however,  the  spectrum  efficiency  of  the  cooperative  network  is  half  of 
the  corresponding  MIMO  system  because  of  the  requirement  of  the  broadcast  phase. 
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B.  Effect  of  Interuser  Channel  Fading 

In  the  second  set  of  simulations,  we  consider  a  scenario  where  710  =  1  and  712  =  720  = 
y/8,  which  matches  the  situation  that  the  relay  terminal  is  in  the  mid-way  between  the 
source  terminal  and  the  destination.  The  power  attenuation  exponential  remains  to  be 
3.  Channels  hia{v)  and  /i2o(^)  are  assumed  Rayleigh  faded.  Two  different  situations 
of  interuser  channel  hnfo)  are  considered:  Rayleigh  faded  and  unfaded.  Compared  to 
the  unfaded  channel,  the  impact  of  channel  fading  in  hniv)  is  evidenced  in  Fig.  6  by  a 
reduction  in  the  diversity  gain  and  a  notable  performance  degradation. 

It  is  noted  that  the  CER  results  depicted  in  Fig.  6,  even  in  the  presence  of  the  interuser 
channel  fading,  are  generally  better  than  those  in  Fig.  5.  That  is,  by  placing  a  relay 
terminal  in  a  proper  location,  power  savings  can  be  achieved  compared  to  a  MIMO  system. 

C.  Significance  of  Power  Allocation 

The  third  set  of  simulations  uses  similar  parameters  as  in  the  second  set,  i.e.,  710  =  1, 
720  =  712  =  \/8,  and  all  the  channels  are  considered  Rayleigh  faded.  To  examine  the  effect 
of  power  allocation,  Fig.  7  compares  the  CER  based  on  optimal  power  allocation  and  a 
fix  power  allocation  with  F\  =  1  ,P2  =  2,  and  a  =  1.  Both  cases  have  the  same  total 
energy  of  four  over  the  four-symbol  interval.  It  is  seen  that,  compared  to  the  optimal 
power  allocation,  the  fixed  power  allocation  yields  an  SNR  loss  of  about  1.5  dB  when  the 
Alamouti  code  is  used,  and  up  to  about  1  dB  when  the  Liang-Xia  code  is  used. 

VII.  Conclusions 

A  novel  space-time  cooperation  protocol  has  been  developed  for  effective  cooperative 
diversity  operation  when  the  channel  state  information  (CSI)  is  unavailable  at  the  re¬ 
ceivers.  It  has  been  shown  that  the  proposed  differential  space-time  cooperation  scheme  is 
effective,  and  the  performance  loss  due  to  differential  detection  is  limited  to  3  dB  of  noise 
enhancement  compared  to  coherent  detection.  We  have  analyzed  the  pairwise  codeword 
error  probability  performance  of  the  cooperative  diversity  system.  Based  on  this  result, 
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the  optimization  of  the  power  allocation  over  the  source  and  relay  terminals  are  considered 
and  the  results  provided  useful  guidelines  in  system  designs. 


Appendix:  Derivation  of  Equation  (41) 


We  first  consider  the  minimum  determinant  of  C'(u  —  1)  A  (C'(u  —  l))H .  For  notation 
simplicity,  we  omit  ( v  —  1)  in  the  rest  of  this  appendix.  Partition  matrix  C'  as 


C'  = 


\/Plc  i 


(A.l) 


where  and  C2  are  defined  in  (43),  and  0  =  [0,0].  Also  denote  O  as  the  2  x  2  zero 
matrix.  Then  we  have 
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Z1  +  Zq  Z2 
Z3  z4  J 

and  the  last  (right-bottom)  element  of  C2(G-E)(G-E)ffCf  equals  to  P2Ci(G-E)(G 
Ei)Hcf ,  then,  we  have 

det(C'A(C')H)  -  1  u-.<^  cm2 


4(1  +  w)2  |det(Q,)det(G-E)| 
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Noting  the  fact  that 

|det(C2)|2  =  det  (C2Cf)  =  det 
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equation  (41)  follows  straightforwardly  from  (A.4). 


CC 


\/A  0 

0  VP2 


} 


=  «2PiP2, 


(A.5) 


42 


References 


[1]  G.  J.  Foschini  and  M.  J.  Gans,  “On  limits  of  wireless  communications  in  a  fading  environment  when 
using  multiple  antennas,”  Wireless  Personal  Commun.,  vol.  6,  no.  3,  pp.  311-335,  March  1998. 

[2]  V.  Tarokh,  H.  Jafarkhani,  and  A.  R.  Calderbank,  “Space-time  block  codes  from  orthogonal  designs,” 
IEEE  Trans.  Inform.  Theory ,  vol.  45,  no.  5,  July  1999. 

[3]  S.  M.  Alamouti,  “A  simple  transmitter  diversity  scheme  for  wireless  communications,”  IEEE  J.  Select. 
Areas  Commun .,  vol.  16,  pp.  1451-1458,  Oct.  1998. 

[4]  A.  F.  Naguib,  N.  Seshadri,  A.  R.  Calderbank,  “Increasing  data  rate  over  wireless  channels:  space-time 
coding  and  signal  processing  for  high  data  rate  wireless  communications,”  IEEE  Signal  Proc.  Mag., 
vol.  17,  no.  3,  pp.  76-92,  May  2000. 

[5]  A.  Sendonaris,  E.  Erkip,  and  B.  Aazhang,  “Increasing  uplink  capacity  via  user  cooperation  diversity,” 
Proc.  IEEE  Int.  Symp.  Info.  Theory,  Cambridge,  MA,  p.  156,  Aug.  1998. 

[6]  A.  Sendonaris,  E.  Erkip,  and  B.  Aazhang,  “User  cooperative  diversity  -  Part  I  and  Part  II,”  IEEE 
Trans.  Commun.,  vol.  51,  no.  11,  pp.  1927-1948,  Nov.  2003. 

[7]  J.  N.  Laneman,  D.  N.  C.  Tse,  and  G.  W.  Wornell,  “Cooperative  diversity  in  wireless  networks:  effective 
protocols  and  outrage  behavior,”  IEEE  Trans.  Inform.  Theory,  accepted  for  publication. 

[8]  J.  N.  Laneman  and  G.  W.  Wornell,  “Distributed  space-time  coded  protocols  for  exploiting  cooperative 
diversity  in  wireless  networks”,  IEEE  Trans.  Inform.  Theory ,  vol.  49,  no.  10,  pp.  2415-2425,  Oct.  2003. 

[9]  A.  F.  Dana  and  B.  Hassibi,  “On  the  power  effeciency  of  sensory  and  ad-hoc  wireless  networks,” 
submitted  to  IEEE  Trans.  Inform.  Theory. 

[10]  M.  Dohler,  E.  Lefranc,  H.  Aghvami,  “Space-time  block  codes  for  virtual  antenna  arrays,”  IEEE  Int. 
Symp.  Personal,  Indoor  and  Mobile  Radio  Communications,  pp.  414-417,  Sept.  2002. 

[11]  P.  A.  Anghel,  G.  Leus,  and  M.  Kaveh,  “Multi-user  space-time  coding  in  a  cooperative  networks,” 
Proc.  IEEE  ICASSP ,  Hong  Kong,  May  2003. 

[12]  Y.  Hua,  Y.  Mei,  and  Y.  Chang,  “Parallel  wireless  mobile  relays  with  space-time  modulations,”  Proc. 
IEEE  Workshop  on  Statistical  Signal  Processing,  St.  Louis,  MO,  Sept.  2003. 

[13]  A.  Stefanov  and  E.  Erkip,  “On  the  performance  analysis  of  cooperative  space-time  coded  systems,” 
IEEE  Wireless  and  Communications  and  Networking  Conf. ,  New  Orleans,  LO,  March  2003. 

[14]  A.  Stefanov  and  E.  Erkip,  “Cooperative  space-time  coding  for  wireless  networks,”  IEEE  Information 
Theory  Workshop,  Paris,  France,  March  -  April  2003. 

[15]  A.  Stefanov  and  E.  Erkip,  “Cooperative  coding  for  wireless  networks,”  IEEE  Trans.  Commun .,  vol. 
52,  no.  9,  pp.  1470-1476,  Sept.  2004. 

[16]  I.  Hammerstroem,  M.  Kuhn,  B.  Rankov,  A.  Wittneben,  “Space-time  processing  for  cooperative  relay 
networks,”  IEEE  Vehicular  Technology  Conf,  Orlando,  FL,  Oct.  2003. 

[17]  R.  U.  Nabar  and  H.  Bolcskei,  “Space-time  signal  design  for  fading  relay  channels”,  GLOBECOM, 
San  Francisco,  CA,  Dec.  2003. 

[18]  R.  U.  Nabar,  F.  W.  Kneubuhler,  and  H.  Bolcskei,  “Space-time  signal  design  for  fading  relay  channels” , 
GLOBECOM,  San  FVancisco,  CA,  Dec.  2003. 

[19]  M.  Janani,  A.  Hedayat,  T.  E.  Hunter,  and  A.  Nosratinia,  “Coded  cooperation  in  wireless  communi- 


43 


cations:  space-time  transmission  and  iterative  decoding,”  IEEE  Trans.  Signal  Processing ,  vol.  52,  no. 
2,  pp.  362-371,  Feb.  2004. 

[20]  Y.  Jing  and  B.  Hassibi,  “Distributed  space-time  coding  in  wireless  relay  networks-Part  I:  basic  diver¬ 
sity  results,”  Submitted  to  IEEE  Trans.  Wireless  Communications. 

[21]  Y.  Jing  and  B.  Hassibi,  “Distributed  space-time  coding  in  wireless  relay  networks-Part  II:  tighter 
upper  bounds  and  a  more  general  case,”  Submitted  to  IEEE  Trans.  Wireless  Communications. 

[22]  B.  L.  Hughes,  “Differential  space-time  modulation,”  IEEE  Trans.  Inform.  Theory ,  vol.  46,  pp.  2567- 
2578,  Nov.  2000. 

[23]  B.  M.  Hochwald  and  T.  L.  Marzetta,  “Unitary  space-time  modulation  for  multiple-antenna  commu¬ 
nications  in  Rayleigh  flat  fading,”  IEEE  Trans.  Inform.  Theory,  vol.  46,  pp.  543-564,  March  2000. 

[24]  B.  M.  Hochwald  and  W.  Sweldens,  “Differential  unitary  space-time  modulation,”  IEEE  Trans.  Com- 
mun .,  vol.  48,  pp.  2041-2052,  Dec.  2000. 

[25]  V.  Tarokh  and  H.  Jafarkhani,  “A  differential  detection  scheme  for  transmit  diversity,”  IEEE  J.  Select. 
Areas  Commun .,  vol.  18,  pp.  1169-1174,  July  2000. 

[26]  M.  Berhler,  and  M.  K.  Varanasi,  “Asymptotic  error  probability  analysis  of  quadratic  receivers  in 
Rayleigh-fading  channels  with  applications  to  unified  analysis  of  coherent  and  noncoherent  space-time 
receivers,”  IEEE  Trans,  on  Inform.  Theory,  vol.  47,  no.  6,  pp.  2383-2399,  Sept.  2001. 

[27]  X.-B.  Liang  and  X.-G.  Xia,  “Unitary  signal  constellations  for  differential  space-time  modulation 
with  two  transmit  antennas:  parameteric  codes,  optimal  designs,  and  bounds,”  IEEE  Trans.  Inform. 
Theory,  vol.  48,  no.  8,  pp.  2291-2322,  Aug.  2003. 

[28]  Z.  Liu,  G.  B.  Giannakis,  and  B.  L.  Hughes,  “Double  differential  space-time  block  coding  for  time- 
varying  fading  channels,”  IEEE  Trans.  Commun .,  vol.  49,  pp.  1529-1539,  Sept.  2001. 

[29]  A.  Shokrollahi,  B.  Hassibi,  B.  M.  Hochwald,  and  W.  Sweldens,  “Representation  theory  for  high-rate 
multiple-antenna  code  design,”  IEEE  Trans.  Inform.  Theory,  vol.  47,  pp.  2335-2367,  Sept.  2001. 

[30]  H.  Li  and  J.  Li,  “Differential  and  coherent  decorrelating  multiuser  receivers  for  space-time-coded 
CDMA  systems,”  IEEE  Trans.  Signal  Processing,  vol.  50,  no.  10,  pp.  2529-2537,  Oct.  2002. 

[31]  E.  Zimmermann,  R  Herhold,  and  G.  Fettweis,  “On  the  performance  of  cooperative  diversity  protocols 
in  practical  wireless  systems,”  IEEE  Vehi.  Tech.  Conf,  pp.  2212-2216,  Oct.  2003. 

[32]  Q.  Zhao  and  H.  Li,  “Performance  of  differential  modulation  with  wireless  relays  in  Rayleigh  fading 
channels,”  IEEE  Communications  Letters,  to  appear. 

[33]  P.  Tarasak,  H.  Minn,  V.  K.  Bhargava,  “Differential  modulation  for  two-user  cooperative  diversity 
systems”,  Globecom,  Nov.-Dee.  2004. 

[34]  D.  Chen  and  J.  N.  Laneman,  “Modulation  and  demodulation  for  cooperative  diversity  in  wireless 
systems,”  submitted  to  IEEE  Trans.  Wireless  Communications. 

[35]  X.  Li,  “Space-time  coded  multi-transmission  among  distributed  transmitters  without  perfect  syn¬ 
chronization,”  IEEE  Signal  Processing  Letters,  vol.  11,  no.  12,  pp.  948-951,  Dec.  2004. 

[36]  Y.  Hua,  A.  Swami,  and  B.  Daneshrad,  “Combating  synchronization  errors  in  cooperative  relays,” 
IEEE  ICASSP,  Philadelphia,  PA,  March  2005. 

[37]  Y.  Zhang,  G.  Wang,  and  M.  G.  Amin,  “Imperfectly  synchronized  cooperative  network  using  distrib¬ 
uted  space-frequency  coding,”  IEEE  Vehicul.  Tech.  Conf.,  Dallas,  TX,  Sept.  2005. 


44 


(a) 


relay 

(b) 


Fig.  l.  System  model. 
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Fig.  2.  Cooperative  diversity  schemes. 
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Fig.  4.  Two-user  system  diagram. 
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Fig.  5.  CER  performance  versus  SNR  for  MIMO  and  cooperative  diversity  systems. 


(a)  Alamouti  code  (b)  Liang-Xia  code 

Fig.  6.  CER  performance  versus  SNR  with  faded  and  unfaded  channels  between  the  source 
and  relay  terminals. 
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Fig.  7.  CER  performance  versus  SNR  for  different  power  allocations. 
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Abstract 

In  distributed  multiple-input-multiple-output  (MIMO)  systems,  imperfect  synchronization  causes  a 
unique  problem  in  a  coded  cooperative  diversity  system.  In  the  presence  of  a  fractional-symbol  delay 
between  the  signals  transmitted  from  different  relay  nodes,  the  channels  become  highly  dispersive  even 
at  a  flat- fading  environment.  Existing  methods  solve  such  problem  based  on  time-domain  approaches 
where  adaptive  equalization  is  required  at  the  receivers  for  combining  the  information  transmitted  from 
distributed  sources.  In  this  paper,  we  propose  the  use  of  OFDM-based  approaches  using  distributed 
space-frequency  codes.  The  proposed  schemes  are  insensitive  to  fractional-symbol  delays  and  lead  to 
higher  data  rate  transmission  and  simplified  implementation.  In  addition,  the  proposed  schemes  permit 
the  use  of  relatively  simple  amplify-and-forward  algorithm  in  multi-hop  wireless  networks  without  delay 
accumulations.  The  proposed  methods  remove  the  time  delay  in  each  relaying  hop  by  reconstructing  the 
prefix  and,  as  such,  improve  the  spectral  efficiency,  while  keeping  a  simplified  relaying  structure. 
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I.  Introduction 


In  a  wireless  network,  information  transmission  through  relaying  can  be  energy  efficient. 
The  use  of  multiple  relaying  nodes  provides  high  diversity  gain,  improving  robustness 
against  channel  impairments  [1].  The  concept  of  network  cooperation  or  cooperative  di¬ 
versity  has  benefitted  from  the  recent  advances  of  space-time  codes,  transmit  diversity,  and 
multi-input-multi-output  (MIMO)  technologies.  Recent  research  work  has  shown  the  fea¬ 
sibility  of  cooperative  protocols  and  provided  various  capacity  and  performance  analyses 
of  cooperative  diversity  systems  (see  for  example,  [2]— [5]  and  references  therein). 

To  achieve  effective  data  transmission  in  a  cooperative  network  surrounded  by  a  flat¬ 
fading  environment,  several  coded  cooperation  schemes  have  been  developed.  However, 
to  adopt  the  detection  schemes  established  for  MIMO  technologies,  a  fundamental  as¬ 
sumption  of  all  these  schemes  is  that  the  relaying  terminals  are  perfectly  synchronized. 
Unlike  the  MIMO  systems,  however,  perfect  synchronization  in  diversity  systems  is  highly 
unlikely  in  practice,  because  of  the  following  reasons  [6]:  (a)  An  accurate  synchroniza¬ 
tion  between  the  relaying  terminals  with  different  locations,  which  are  probably  subject 
to  continuous  movements,  is  impractical,  (b)  Even  for  a  network  with  fixed  terminals 
(such  as  in  a  sensor  network),  changes  in  the  parameters  of  electronic  components  may 
result  in  drifting  and  handshaking  among  transmitters  is  done  infrequently  so  as  to  save 
energy  and  bandwidth,  (c)  Although  nodes  can  synchronize  the  received  signals  as  much 
as  possible,  low  cost  implementations  may  still  make  their  timing  and  frequency  slightly 
different,  hence  cause  mismatch  in  the  long  run.  (d)  The  major  synchronization  problem, 
however,  lies  in  the  delays  of  their  signals  when  reaching  at  the  receiving  nodes  in  the 
next  hop.  Propagation  delays  may  be  unknown  to  them,  while  transmission  times  may 
be  different.  In  fact,  if  the  transmitting  nodes  try  to  synchronize  toward  one  receiver, 
they  may  increase  asynchronism  toward  other  receivers  because  of  different  transmission 
distances.  On  the  other  hand,  in  most  practical  networks,  nodes  can  (and  are  required)  to 
maintain  slot  synchronization,  which  means  that  coarse  slot  synchronization  is  available. 
The  difficulty  lies  in  the  fine  transmission  synchronization. 

Most  digital  modulation  signals  guarantee  no  inter-symbol  interference  (ISI)  only  when 
the  signals  are  accurately  sampled  at  the  symbol  center.  When  the  signal  arrivals  from 
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different  relay  terminals  are  not  perfectly  synchronized,  the  effect  due  to  pulse  tails  resulted 
from  sampling  position  error  should  be  considered.  This  is  particularly  true  when  the 
receiver  performs  a  symbol-based  processing.  To  consider  such  effect,  [6]  proposed  the 
use  of  a  equivalent  dispersive  channel  model,  whether  the  channels  are  frequency-selective 
or  frequency-nonselective.  It  also  proposed  a  useful  scheme  which  employs  distributed 
space-time  block  code  (STBC)  in  a  time-reverse  manner  to  achieve  full  diversity  while 
tolerating  imperfect  synchronization.  The  difficulty  of  using  this  method,  nevertheless, 
lies  in  the  fact  that  it  requires  a  channel  equalizer  at  each  receiver  that  makes  the  receiver 
very  complicated  and  expensive. 

In  this  paper,  we  propose  the  use  of  multi-carrier  modulation  exploiting  distributed 
space-frequency  coding  schemes.  The  proposed  schemes  allow  higher  data  rates  and  sim¬ 
plified  implementation,  as  well  as  the  distribution  of  the  system  complexity  to  transmitters 
and  receivers.  In  addition,  the  feasibility  of  distributed  space-frequency  coding  in  a  multi¬ 
hop  relaying  cooperative  network  is  also  considered.  It  is  pointed  out  that  the  use  of 
MIMO-OFDM  methods  allows  the  utilization  of  simple  amplify-and-forward  algorithm. 
While  the  time  delay  over  different  relay  hops  will  be  accumulated  in  the  amplify-and- 
forward  relaying  process  when  space-time  coding  methods  (e.g.,  time-reverse  STBC,  or 
TR-STBC  [7])  are  used,  such  delay  accumulation  problem  can  be  easily  avoided  in  the 
proposed  method  through  the  reconstruction  of  the  cyclic  prefix. 

II.  Channel  Model 

Consider  a  wireless  network  where  a  source  node  transmits  a  data  packet  to  a  destination 
node  through  a  multi-hop  wireless  network.  In  each  intermediate  hop,  the  data  packet  is 
received  by  multiple  nodes,  e.g.,  nodes  1  to  J  in  hop  1.  These  nodes  can  retransmit  the 
data  in  the  next  hop  in  a  cooperative  maimer  with  proper  space-time  (ST)  encoding. 

Without  loss  generality,  the  noise-free  signals  received  at  the  kth  node  is  expressed  as 

Xk (t)  =  f  JPihi>k(t,  r)si(t  -  r  -  TD^k)dr ,  (1) 

i=  i  J~°° 

where  Si(t)  is  the  signal  waveform  transmitted  from  node  i,  hitk{t,  r)  and  To,i,k  are  the 
channel  coefficients  and  the  time  delay  between  transmit  node  i  and  receive  node  k,  re¬ 
spectively.  It  is  assumed  that  channels  hi>k{t,  r)  are  time-invariant  during  the  transmission 


51 


of  a  packet,  but  are  randomly  time- varying  between  packets.  The  channels  are  assumed 
spatially  white,  i.e.,  r)  are  independent  for  different  indices  (i,k).  In  addition,  \/7{ 
is  the  transmit  power  of  node  i.  In  a  distributed  MIMO  system,  the  transmit  power  of 
different  nodes  may  be  different,  depending  on  the  a  priori  information  of  the  large-scale 
attenuation  characteristics  at  the  transmitters  [3],  [12],  The  total  power  is  assumed  to  be 
a  constant,  i.e.,  S/=1  Pi  —  P. 

Let  Ii(n)  as  the  discrete  information  sequence  transmitted  from  at  the  ith  node,  and 
p(t)  as  the  pulse  shaping  function.  Then,  (1)  becomes 


J  OO  -  /*oo 

Xk(t)  =  53  V  Pih{m)  J  hitk(t ,  r)p(t  -r-  TD<iik  -  mTs)dr. 


i=l  m=— oo 


In  general,  the  time  delay  consists  of  three  factors:  channel  dispersion  because  of  reflection 
and  scattering  ( rc,i,k ),  delay  due  to  different  locations  of  the  relaying  nodes  ( 'jD,i,k )>  and  the 
pulse  shaping  function  spreading  due  to  sampling  position  errors  ( Tw,i,k )•  Therefore,  the 
significance  of  IS  I  depends  not  only  on  the  channel  length,  but  also  on  the  node  locations 
and  the  pulse  shaping  functions.  Combining  all  these  factors,  we  denote  the  upper  bound 
of  the  effective  channel  length  in  terms  of  symbol  period,  as  L  =  |max(rc,-i,fc) +max(rDji|fc)  + 
ma x{Tw,i,k)]/Ts  with  T,  denoting  the  symbol  period.  We  use  Lc  =  max(ro,i,fc)/Ts,  LD  — 
max(rJD|iifc)/Ts,  and  Lw  =  max(rH/,ifc) /Ts  for  notational  simplicity,  and  L  =  Lc+Ld+Lw. 

Sampling  at  t  =  nTs,  and  taking  the  receiver  noise  into  account,  we  have 


xk(n)  =  xk(nTs )  =  5Z  12  V  /  Kk(nTs ,  r)p[(n  -  m)Ts  -  r  -  TD>i,k]dT  +  vk(n). 

i=  1  m=—oo 

(3) 

The  noise  is  assumed  to  be  of  zero  mean  and  unit  variance,  and  is  temporally  and  spatially 
white,  i.e.,  E[vk(r)vi  (u)]  =  0  for  k  ^  /  or  r  ^  u,  where  *  denotes  complex  conjugation. 


III.  Distributed  MIMO-OFDM  Schemes 
A.  MIMO-OFDM  Transmission 

We  first  review  the  MIMO-OFDM  transmission  using  space-frequency  (SF)  codes  in  a 
transmitter  with  J  co-located  antennas.  The  following  J  x  N  SF  codeword  is  formed  in 
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the  transmitter 


cil  rci(O)  ci(l)  •••  Ci(N-l) 

c2  c2(0)  c2(l)  •••  c2(A- 1) 

CsF  =  .  =  . 

.Cj.  .Cj(0)  Cj(  1)  cj(N  —  1) 

where  Ci(n)  denotes  the  channel  symbol  transmitted  over  the  nth  subcarrier  by  transmit 
antenna  i,  i  =  1,  •  •  • ,  J,  and  n  =  0,  -  •  • ,  A  —  1.  The  construction  of  space-frequency  codes 
has  been  discussed  in,  for  example,  [8]— [10].  A  detailed  review  of  broadband  MIMO-OFDM 
wireless  communications  is  provided  in  [11], 

The  OFDM  transformation  of  the  source  codeword  applies  an  A-point  IFFT  to  each  row 
of  the  matrix  Csf-  To  eliminate  the  effect  of  ISI,  cyclic  prefix  is  inserted  in  the  beginning 
of  the  transformed  matrix.  The  OFDM  symbol  corresponding  to  the  ith  row  of  Csf  is 
then  transmitted  using  the  ffh  transmit  antenna. 

At  the  receiver,  after  matched  filtering,  removing  the  cyclic  prefix,  and  applying  FFT 
to  recover  the  transmitted  SF  codeword,  information  can  be  decoded,  for  example,  using 
maximum  likelihood  detection. 

B.  Transmission  of  Distributed  SF  Codes 

Now  we  extend  the  SF  coded  MIMO-OFDM  approaches  to  the  distributed  wireless 
network.  Similar  to  the  protocol  proposed  in  [3],  [4]  for  constructing  distributed  ST  codes, 
distributed  SF  codes  can  be  constructed  at  the  source  terminal  during  the  first  hop.  This 
paper  only  considers  the  intermediate  hops  where  information  is  relayed  from  multiple 
distributed  transmit  nodes  to  multiple  distributed  receive  nodes,  forming  a  distributed 
MIMO  structure.  In  particular,  a  relay  node  transmits  a  specific  row  of  a  SF  codeword. 
The  proposed  method  can  be  used  for  both  flat  and  dispersive  channels.  In  the  latter, 
path  diversity  can  be  achieved  [9],  [10]. 

The  OFDM  transformation  at  the  source  codeword  applies  an  A-point  IFFT  to  each 
row  of  the  matrix  Csf,  resulting  in 

Xsf  =  CsfF,  (5) 

where  F  is  the  A  x  A  inverse  Fourier  transform  matrix.  After  appending  a  cyclic  prefix 
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of  length  C,  the  J  x  (N  +  C)  matrix  is  obtained  as 

XgF  —  [xSf(./V  C) ,  •  •  •  j  xSf(AT  1) ,  XSf]  ,  (6) 

where  xsp(n)  is  the  nth  column  of  Xsp-  The  frame  length  becomes  N  +  C  symbol 
durations  after  the  cyclic  prefix  is  inserted.  The  OFDM  symbol  corresponding  to  the  ith 
(i  =  1  J)  column  of  X'SF  is  transmitted  by  the  source  node  through  the  ith  time 
frame. 

C.  Transmission  Efficiency 

In  the  single-carrier  time-domain  approaches,  the  minimum  guard  interval  is  Lc  +  Lq  + 
Lw.  In  the  proposed  OFDM  approach,  the  effective  pulse  tail  Lw  does  not  play  a  role.  As 
a  result,  the  guard  interval  is  solely  determined  by  the  channel  order,  Lc,  and  the  upper 
bound  of  the  time  delay  between  different  transmit  nodes,  LD.  It  could  imply  a  significant 
reduction,  particularly  in  flat  or  not  heavily  dispersive  channels,  since  Lw  may  take  a  large 
value,  depending  on  the  synchronization  strategy  and  the  used  pulse  shaping  function.  By 
reducing  the  time  delay  effect  from  Lc  +  Lw  +  To  to  Lc  +  Lc,  the  guard  interval  can  be 
designed  to  be  smaller  than  that  in  the  single-carrier  counterparts,  resulting  in  improved 
data  rates. 

D.  System  Complexity 

The  proposed  method  requires  transmitters  and  receivers  to  perform  IFFT  and  FFT  for 
OFDM  modulations.  It  is  usually  a  great  reduction  to  the  computational  cost  of  an  adap¬ 
tive  equalizer.  It  is  noticed  that  the  proposed  approach  requires  IFFT  or  FFT  operations 
in  both  transmit  /receive  sides,  whereas  the  single-carrier  time-domain  approaches  only 
require  equalization  performed  at  the  receivers.  Therefore,  consider  a  multi-hop  wireless 
network  using  the  decode-and-forward  algorithm,  the  proposed  OFDM  method  provides 
significant  advantages  over  the  single-carrier  approaches  because  the  later  demands  high 
burden  for  equalizations  in  the  decoding  process.  On  the  other  hand,  when  we  adopt  the 
amplify-and-forward  algorithm  in  the  wireless  network,  single-carrier  relaying  approaches 
accumulate  the  channel  delays  over  multiple  hops,  whereas  the  proposed  method  can  easily 
remove  such  delay  accumulation  by  reconstructing  the  cyclic  prefix.  The  reconstruction 
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of  the  cyclic  prefix  is  discussed  in  the  following  section. 


IV.  Cyclic  Prefix  Reconstruction  in  Multi-Hop  Amplify-and-Forward 

Relaying 

At  the  kth  relaying  node,  the  received  signal  corresponding  to  the  nth  subcarrier  is 
given  by 

Vk{n)  =  Y  \fPA(n)Hitk(n)  +  vk{n ),  (7) 

i=l 

where  c*(n)  is  the  scaled  and  noisy  replica  of  c*(n)  contaminated  in  the  previous  hops, 
Hi>k{n)  is  the  Fourier  transform  of  hlik(t),  and  vk{n)  denotes  the  additive  complex  Gaussian 
noise  at  the  nth  subcarrier. 

In  matrix  format,  for  the  nth  subcarrier,  the  received  signal  at  the  J'  relaying  nodes, 
after  the  prefix  is  removed,  is  expressed  as 

y(n)  =  H(n)P1/2c(n)  +  v(n),  (8) 


where 


and 


y(n)  -  [yi(n),y2{n)i---  ,yj(n)}T , 

c(n)  =  [cl(n),c2(n),---  ,Cj{n)]T , 

v(n)  =  Mn),u2(n),  •  •  •  ,  Vj{n)]T , 
hi,i(n),  h2,i(n ),  •••,  hjA(n) 

hi,  2(^)1  h2, 2(n),  •••,  hj)2(n) 

*  5  *  ? 

hi,j'(n),  h2,j'{n ),  -••,  hj^in) 

P  =  diag[Pi,P2,---,Pj]. 


The  signal  is  retransmitted  after  the  prefix  is  reconstructed  from  y(n).  At  the  next  hop, 
the  received  signal  vector  at  the  relaying  nodes  is  expressed  as 
y'(n)  =  H'(n)  jH(n)P1//2c(n)  +  v(n)|  +  v;(n) 

=  H' (n) H (n) P 1,/2 c (n)  +  [H'(n)v(n)  +  v'(n)] ,  (9) 

where  H'(n)  is  the  channel  matrix  corresponding  to  the  next  hop.  In  general,  when  M 
hops  are  involved,  the  above  process  can  be  repeated.  Denoting  Hm(n)  as  the  channel 
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matrix  corresponding  to  the  mth  hop  and  the  nth  subcarrier,  it  is  seen  that  the  proposed 
amplify-and-forward  scheme  relays  SF  coded  information  with  equivalent  channel  matrix 
nJJ-.IMn),  whereas  the  noise  is  accumulated.  The  information  can  be  decoded  using 
standard  coherent  detection  methods  if  the  combined  channel  state  information  (CSI) 
can  be  estimated  at  the  destination,  or  using  noncoherent  detection  methods  that  do  not 
required  CSI. 

It  is  emphasized  that,  unlike  the  single-carrier  amplify-and-forward  process  where  the 
time  delay  over  different  hops  will  be  accumulated,  resulting  in  unrealistic  implementations 
in  multi-hop  systems,  the  effect  of  the  time  delay  in  each  hop  is  decoupled  by  the  insertion 
and  subtraction  of  the  prefix,  rendering  to  a  manageable  frame  design  and  high  spectral 
efficiency. 

V.  Numerical  Results 

A.  Effect  of  Pulse  Shaping  Functions  1 


To  illustrate  the  effect  of  pulse  shaping  functions  in  an  dispersive  or  imperfectly  synchro¬ 
nized  network  using  the  single-carrier  time-domain  approaches,  Fig.  1  depicts  the  impulse 
response  of. raised  cosine  filters  with  different  sampling  position  errors  of  0,  Ts/ 4,  and  Ts/ 2, 
where  the  roll-off  factor  is  0.3.  Denote  p(l )  as  the  impulse  response  of  the  raised  cosine 
filter  sampled  at  the  symbol  rate,  and  define 


e(n) 


E  |p(OI  -  E  Ml) 

all  symbols  n  dominant  symbols 

E  woi 

all  symbols 


(10) 


as  the  residual  error  of  the  impulse  response  after  n  dominant  symbols  are  counted,  the 
result  is  summarized  in  Table  I.  It  is  clear  that,  the  impulse  response  spans  roughly  five  to 
six  symbols  for  both  sampling  position  errors  of  Ts/4  and  Ts/2.  The  long  impulse  response 
results  in  high  equalizer  complicity  and,  in  a  TR-STBC  structure,  a  long  guard  interval. 

The  effect  of  sampling  position  errors  can  be  reduced  by  using  a  large  value  of  roll-off 
factor.  Fig.  2  depicts  the  impulse  response  of  raised  cosine  filters  using  roll-off  factor  of 
1.0.  While  it  implies  a  reduction  in  equalizer  complexity  and  guard  interval,  the  use  of 
large  roll-off  factor  results  in  a  significant  reduction  of  the  spectral  efficiency. 

Unlike  the  time-domain  approaches,  the  effect  of  pulse  shaping  in  an  OFDM  system  is 
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minimal  As  a  result,  the  effect  of  sampling  position  errors  to  the  proposed  MIMO-OFDM 
approach  is  negligible. 

B.  Effect  of  Cyclic  Prefix  Reconstruction 

We  consider  a  portion  of  a  multi-hop  relaying  system  as  shown  in  Fig.  3.  Two  relaying 
hops  (hop  i  and  hop  i  +  1)  are  included  in  this  diagram  to  illustrate  the  capability  of 
eliminating  delay  accumulation  in  the  proposed  scheme.  Each  relaying  hop  includes  two 
transmit  antennas  (Mt  —  2).  Only  one  receiver  is  assumed  at  the  receive  end  of  the 
( i  +  l)th  hop,  because  the  performance  of  each  node  should  be  evaluated  separately. 
While  the  relay  nodes  do  not  decode  the  signal,  the  bit  error  rate  (BER)  at  the  node  is 
computed  to  measure  the  quality  of  the  received  signal. 

The  channels  are  considered  to  be  independent  random  processes.  It  is  assumed  that 
the  channels  are  fixed  over  a  frame  period,  but  vary  randomly  and  assume  complex  Gau¬ 
ssian  distribution  over  different  frames.  The  channels  are  normalized  such  that  they  have 
a  unit  variance.  In  each  hop,  the  total  power  is  unity  and  is  equally  distributed  over  the 
two  transmit  antennas.  The  noise  power  is  assumed  to  be  the  same  over  different  relaying 
hops.  The  input  signal-to-noise  ratio  (SNR)  is  defined  as  the  reciprocal  of  the  receiver 
noise  power  at  each  hop. 

128  subcarriers  are  used  for  the  OFDM  signals,  and  the  guard  interval  consists  of  six 
symbols,  resulting  in  total  134  symbol  periods  over  an  OFDM  frame.  The  distributed 
SF  codewords  are  generated  based  on  [10]  where  F  =  4  and  P  =  16  are  used.  The 
source  stream  before  SF  coding  is  a  binary  sequence.  Random  permutation  is  applied  in 
generating  the  SF  codewords  to  achieve  path  diversity.  5000  frames  are  used  to  evaluate 
the  BER  performance. 

The  order  of  each  channel  is  assumed  to  be  five  symbols.  Therefore,  without  prefix 
reconstruction,  the  time  delay  over  the  two  hops  will  be  accumulated,  resulting  a  total 
channel  order  of  ten  symbols.  As  a  result,  the  six-symbol  prefix  does  not  provide  enough 
protection  against  the  channel  dispersion,  resulting  in  a  BER  floor  if  no  equalization  is 
applied  (the  solid  line  in  Fig.  4).  However,  when  the  prefix  is  removed  and  reconstructed 
in  each  relay,  the  time  delay  can  be  eliminated  in  each  relay.  It  is  evident  in  Fig.  4  (the 
dashed  line)  that  the  BER  in  such  case  does  not  have  a  floor  and  a  high  diversity  gain  is 
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TABLE  I 


Residual  Error  of  Impulse  Response 


n 

1 

2 

3 

4 

5 

6 

Tsf  4 

0.440 

0.261 

0.163 

0.101 

0.0617 

0.0364 

Ts/ 2 

0.658 

0.315 

0.219 

0.123 

0.0838 

0.0442 

achieved  in  the  relay  process. 


VI.  Conclusions 

In  a  cooperative  wireless  network,  channel  dispersion  exists  as  the  result  of  several 
factors,  e.g.,  channel  length,  pulse  shaping  function  spreading,  and  spatial  separations 
between  the  transmit  nodes.  We  have  considered  the  coded  data  transmission  in  a  cooper¬ 
ative  network  with  frequency-selective  fading  and  proposed  the  distributed  MIMO-OFDM 
schemes  for  effective  multi-hop  relaying.  It  is  shown  that,  by  reconstructing  the  cyclic 
prefix,  the  time  delay  effect  in  each  relaying  hop  can  be  eliminated  even  using  relatively 
simple  amplify-and-forward  algorithm.  Therefore,  the  guard  interval  can  be  minimized 
and,  thereby,  high  spectral  efficiency  is  maintained. 
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time  (symbol) 


(a)  no  sampling  position  error 


time  (symbol) 


(b)  sampling  position  error  —  Tsj 4 


time  (symbol) 


(c)  sampling  position  error  =  Ts/2 

Fig.  1.  Impulse  response  of  raised  cosine  filter  with  different  sampling  position  errors 
(roll-off  factor  is  0.3). 
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time  (symbol) 


(a)  no  sampling  position  error 


time  (symbol) 


(b)  sampling  position  error  —  Tsj 4 


time  (symbol) 


(c)  sampling  position  error  —  Tsj 2 

Fig.  2.  Impulse  response  of  raised  cosine  filter  with  different  sampling  position  errors 
(roll-off  factor  is  1.0). 
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o 
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Fig.  3.  Multi-hop  relaying  system  under  consideration. 


Fig.  4.  BER  with  and  without  prefix  reconstruction. 
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Space-Time  Code  Designs  with  Non- Vanishing  Determinants 
Based  on  Cyclic  Field  Extension  Families  * 

Genyuan  Wang,  Jian-Kang  Zhang,  Yimin  Zhang,  Moeness  Amin,  and  Kon  Max  Wong 


Abstract 

The  design  of  a  linear  space-time  code  with  full  rate,  large  diversity  product,  and  non¬ 
vanishing  minimum  determinant  of  codewords  continues  to  attract  great  attention.  However, 
in  most  available  no-vanishing  determinant  space-time  codes  for  three,  four  and  six  transmit 
antennas,  the  average  powers  at  each  layer  are  different,  which  results  in  a  high  peak  to  average 
power  ratio.  In  this  paper,  a  new  cyclic  algebraic  space-time  design  scheme  is  proposed  and 
the  optimal  codes  in  this  class  are  provided  by  using  some  specific  cyclic  field  extensions.  Our 
proposed  codes  not  only  include  the  available  non- vanishing  determinant  cyclotomic  space-time 
codes  for  three,  four  and  six  transmit  antennas,  but  also  have  the  desirable  property  that  the 
optimal  codes  can  be  achieved  with  the  same  average  power  at  each  layer. 

Keywords:  Multi-layer  space-time  block  codes,  diversity  product,  full-rate,  algebraic  num¬ 
ber  theory,  cyclotomic  number  rings,  cyclic  field  extensions,  lattices. 
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1  Introduction 


Linear  space-time  block  code  designs  based  on  algebraic  field  extensions  have  recently  attracted 
great  attention,  see  for  example  [1]-[12],  due  to  the  possibility  of  systematic  constructions  of  full 
diversity  and  high  data  rate  codes.  Diagonal  algebraic  space-time  block  codes  were  first  proposed 
in  [3],  where  an  n-dimensional  diagonal  space-time  code  diag([yi,y2>  •  -  •  , yn] )  was  generated  by 
[yu  ■  -  •  ,  yn\r  —  G[x  i,  •  •  •  ,xn\T,  with  matrix  G  and  transmitted  symbols  X\,X2>  ■■■  ,xn  being  prop¬ 
erly  chosen  based  on  algebraic  extension  theory  to  achieve  full  diversity.  The  idea  behind  the 
diagonal  algebraic  space-time  code  can  be  tracked  back  to  [1,  2],  where  the  full  diversity  multi¬ 
dimensional  signal  constellation  designs  in  both  Rayleigh  fading  and  additive  Gaussian  noise  chan¬ 
nel  were  considered.  However,  the  symbol  rate  for  the  above  diagonal  space-time  code  design  is  one 
per  channel  use.  In  [6],  a  full  diversity  space-time  code  for  two  transmit  antennas  was  proposed, 
where  the  symbol  rate  reached  two  per  channel  use.  By  employing  algebraic  number  theory  and 
threaded/multi-layer  codes  [14],  more  general  full  diversity  and  high  symbol  rate  space-time  code 
designs  were  proposed  in  [4,  6,  7,  11,  10],  Meanwhile,  another  type  of  full  diversity  and  high  rate 
space-time  code  was  presented  in  [9]  based  on  cyclic  field  extensions  and  division  algebras  . 

In  the  early  studies,  the  structure  of  code  designs  with  high  (full)  rate  and  full  diversity  received 
more  attention  than  the  high  diversity  product.  In  most  existing  codes,  the  minimum  determinant 
of  non-zero  codewords,  which  is  the  minimum  determinant  of  any  two  different  codewords,  vanishes 
as  the  symbol  constellation  size  increases. 

Other  space-time  codes  with  a  full  symbol  rate  and  high  diversity  product  have  been  recently 
generated  in  [17,  18,  19,  21],  These  codes  not  only  have  high  diversity  products,  but  also  have 
the  non- vanishing  determinant  property;  i.e.,  the  minimum  determinant  does  not  decrease  as  the 
symbol  constellation  size  increases.  Although  the  codes  in  [17,  19]  for  two  transmit  antennas  have 
the  same  average  powers  at  different  layers,  the  cyclotomic  space-time  codes  in  [21]  and  [18]  for 
three,  four  and  six  transmit  antennas  have  different  average  powers  at  different  layers,  resulting  a 
higher  peak  to  average  power  ratio.  In  this  paper,  we  propose  a  more  general  full  diversity  and 
full  rate  linear  space-time  code  design  with  the  non-vanishing  determinant  based  on  some  specific 
cyclic  field  extensions,  which  is  called  a  cyclic  algebraic  space- time  code  design.  This  class  of  cyclic 
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algebraic  space-time  codes  includes  all  the  cyclotomic  space-time  codes  in  [18]  and  [21]  for  three, 
four  and  six  transmit  antennas.  In  addition,  the  optimal  cyclic  algebraic  space-time  codes  with 
largest  diversity  products  were  constructed  with  the  same  average  power  at  each  layer.  1 

This  paper  is  organized  as  follows.  In  Section  2,  the  motivation  and  the  problem  description  of 
this  work  are  provided.  In  Section  3,  the  cyclic  algebraic  space-time  code  based  on  some  specific 
cyclic  field  extension  is  introduced.  The  codewords  with  improved  performance  and  the  same 
average  power  at  each  layer  are  then  given.  Simulation  results  are  presented  in  Section  4. 

All  lengthy  proofs  of  the  theorems  are  given  in  Appendix. 

The  following  notation  is  used  throughout  this  paper:  capital  English  letters,  such  as,  X  and 
G,  represent  space-time  codeword  or  matrix, 

Lt:  number  of  transmit  antennas 

N:  natural  numbers 

Z:  ring  of  integers 

Q:  field  of  rational  numbers 

C:  field  of  complex  numbers 

<p(n):  Euler  function  of  positive  integer  n 

<m  =  exp  (j%) 

K,F :  general  fields 

F(/?):  field  generated  by  fi  based  on  field  F 

(K/F,  /?,  a):  cyclic  field  extension  K/F  with  K  =  F(/?),  and  a  of  the  generator  of  cyclic  Galois  group 
Gal(  K/L) 

A  (K/F,  fi,  <7,  p):  space-time  code  generated  with  cyclic  field  extension  (K/F,  (3,  a),  l,p,  •••  ,pra_1 
are  the  numbers  used  for  adjusting  at  different  layers  of  the  code. 

2  Motivation  and  Problem  Description 

Before  providing  the  cyclic  algebraic  space-time  code  design,  we  first  review  cyclotomic  space-time 
designs,  based  on  which  most  of  full  rate  non-vanishing  determinant  linear  space-time  codes  are 

1  Before  the  submission  of  this  paper,  the  authors  learned  that  G.  Rekaya,  J.  C.  Belfiore  and  E.  Viterbo  recently 
also  developed  a  non-vanishing  determinant  space-time  code  design  for  three,  four  and  six  transmit  antennas  with 
same  average  power  at  different  layers  [20], 
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developed. 


2.1  Full  Rate  Cyclotomic  Space-Time  Codes  Designs 


In  this  subsection,  we  recall  cyclotomic  lattices,  cyclotomic  space-time  codes  and  some  of  their 
fundamental  properties  obtained  in  [1,  2,  12].  For  two  positive  integers  m  and  n,  let  N  =  ran  and 


r  MW 

*  tfm)* 


(1) 


The  variable  Lt  corresponds  to  the  number  of  transmit  antennas  in  a  space-time  code.  There  are 
totally  Lt  distinct  integers  k,  1  <  i  <  Lt,  with  0  =  h  <  h  <  ■  ■  •  <  <  n  —  1  such  that  1  +  km 

and  N  are  co-prime  for  any  1  <  i  <  Lt  (see  for  example  p.  75  of  [28]).  With  these  Lt  integers  and 
N  =  mn ,  we  define 


C  N 

Cl  ■■■ 

^2  (I+Z2m) 

CLt  1 

hN 

(2) 

A+lLtm 

S  N 

.2(1  +lLtm) 

SN 

Lt  x  Lt 

It  can  be  verified  that  the  above  G'm,n  is  unitary  when  n  —  Lt.  We  now  define  cyclotomic  lattices. 


Definition  1  An  Lt- dimensional  cyclotomic  lattice  T Lt{Gm,n)  is  a  set  of  Lt-dimensional  points 
[yi,---  ,yLt]T  such  that 

{y  1 ,  -  ,yLt]T  =  Gm,n[x I,-'  -  ,XLt]T,  (3) 

where  Gm,n  is  defined  in  (2)  and  xi  € 


Following  the  structure  of  threaded  space-time  codes  in  [14],  a  general  multi-layer  cyclotomic  space- 
time  code  is  defined  as  follows. 


Definition  2  Let  Lt  be  the  number  of  transmit  antennas  and  T Lt(Gmuni)  be  an  Lt-dimensional 
cyclotomic  lattice  given  in  Definition  1,  where  Gmi>ni  is  defined  in  (2)  for  l  =  1,2, •••  , Lt-  Let 
pi,  ■  •  ■  ,pit  be  Lt  fixed  complex  numbers.  Then,  a  multi-layer  cyclotomic  space-time  code  is  defined 
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as 


PlVlO)  P2|fe(l)  •••  pLt-iVLt-iO)  pLtVLti  1) 

PLtULti 2)  Piyi(2)  Pit— 22/Lt— 2(2)  PLt-lVLt-li 2) 

p2D2{Lt)  p3V3(Lt )  •••  pLtVLti^t)  Pm(Lt) 

p\X\  P2X2  pLt-lXLt-l  PLtXLt 

PLt°l(XLt )  Pl<^l(^l)  •••  PLt—2<X\{XLt—2)  PLt-l^liXLt-l) 


(4) 


L  P20*it-l(®2)  P3<TLt-l(a:3)  PLt<TLt-l(XLt)  Pl°Lt-l(Xl )  J 

where  [yfi  1),  •  -  •  ,  y;(Lt)]T  is  a  point  in  cyclotomic  lattice  r Lt(C?mi,n1)  and  oi  are  the  Lt  embeddings 
ofQiCmn)  to  C  that  is  fixed  on  Q(Cm)  /or  Z  =  1, •••  This  multi-layer  cyclotomic  space-time 
code  is  denoted  by  X(piGmi>ni,  •  •  ■  ,PLtGmLt,nLt)-  An  L-layer  cyclotomic  space-time  code  with 
1  <  L  <  Lt  is  defined  as  a  multi-layer  cyclotomic  space-time  code  X(piGmitni,  ■  ■  ■  ,  p[JtGmLl,nLt) 
when  pi  —  0  for  l  >  L  and  is  denoted  by  X(piGmi,n1,-  •  •  ,  PL,GmL,nL )■ 


2.2  Problem  in  the  Cyclotomic  Space-Time  Code  Designs 


The  key  to  designing  a  full  rate  cyclotomic  space-time  code  is  to  find  a  proper  cyclotomic  field 
extension  Q((mn)/Q((m)  for  given  Lt  and  pi, l  =  1,  •••  ,Lt,  in  (4)  such  that  the  resulting  code 
achieves  full  diversity  and  a  large  diversity  product.  Prom  the  early  studies  on  this  topic,  we  know 
that  there  are  infinite  ways  of  determining  the  values  of  m,  n,  and  pi  to  generate  full  diversity  full 
rate  space-time  codes.  However,  most  of  these  codes  do  not  have  a  large  diversity  product  and  the 
minimum  determinant  of  codewords  decrease  to  zero  very  rapidly  when  the  constellation  size  of  the 
codewords  increases. 

Recently,  the  generation  of  large  diversity  product  space-time  codes  with  a  non-vanishing  deter¬ 
minant  has  received  much  attention.  The  vital  point  to  design  this  type  of  codes  is  to  find  a  proper 
field  extension  K/F  and  pi, l  =  1,  •  •  •  ,Lt,  such  that  the  determinant  of  every  codeword  belongs  to 
the  same  lattice.  The  codewords  provided  in  [17,  18,  19,  21]  have  the  property  that  F  =  Q((4) 
or  Q(Ca)  =  QOSs)?  pi  €  ^[(4]  or  Zf^J  =  Z^],  and  the  determinant  of  every  codeword  belongs 
to  Z [(,„].  However,  in  the  code  design  of  [18,  21],  pi  =  pl,l  =  1,  •  •  •  ,Lt,  with  |p|  >  1.  Therefore,  ]p/| 
takes  different  values  for  different  Z;  i.e.,  the  average  powers  of  the  codewords  at  different  layers  are 
different.  This  results  in  a  high  peak  to  average  power  ratio.  A  close  examination  of  the  code  design 
in  [18,  21]  reveals  that  the  field  extension  Q(Cmn/Q(Cm)  based  on  minimal  polynomial  xn  —  (m, 
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m  =  3, 4,  or  6,  and  pi  may  not  be  a  good  choice.  In  this  paper,  we  develop  a  code  design  using  a 
more  general  field  extension  K/Q((m)  with  the  minimal  polynomial  xn  —  a  for  a  6  Z[(m]  and  pi- 
The  resulting  new  space-time  codes  have  non-vanishing  determinants  with  the  same  average  power 
at  each  layer;  i.e.,  \pi\  —  1,  for  three,  four  and  six  transmit  antennas. 

3  Fall  Rate  Cyclic  Algebraic  Space-time  Codes  for  Three,  Four, 
and  Six  Transmit  Antennas  with  Low  Peak  to  Average  Power 
Ratios 

The  existing  cyclotomic  code  designs  are  obtained  based  on  the  cyclic  field  extension  (Q((m™)/Q(Cn),  Cron) 
with  the  minimal  polynomial  xn  —  (m,  m  =  6  or  m  =  3.  In  this  section,  we  generalize  the  cyclic  field 
extension  to  the  (Q((mn)/Q(Cn),  /?)  with  a  minimal  polynomial  xn  —  a,  and  (3n  —  a  e  We 

find  optimal  full  rate  space-time  codes  for  three,  four,  and  six  transmit  antenna  in  this  class.  The 
results  show  that  optimal  full  rate  space-time  codes  can  be  achieved  with  \p\  =  1,  which  means  that 
the  average  powers  at  different  layers  of  the  codewords  are  identical;  i.e.,  a  lower  peak  to  average 
power  ratio  than  that  achieved  by  the  cyclotomic  space-time  codes.  Before  providing  the  optimal 
codes,  we  introduce  some  concepts  and  results. 

Definition  3  [28]  A  Galois  extension  K/F  is  called  cyclic  if  the  Galois  group  Gal(l K/F)  is  a  cyclic 
group. 

By  the  cyclic  extension  theory  we  have  the  following  proposition. 

Proposition  1  [28]  Let  ¥  .be  a  field  containing  a  primitive  nth  root  of  unity,  and  letK  =  F(  y/a) 
for  some  a  €  F.  Then  K/F  is  a  cyclic  Galois  extension.  Moreover,  m  =  [K :  FJ  is  equal  to  the  order 
of  the  coset  a¥*n  in  the  group  F*/F*n,  and  min(F,  yfa)  =  xm  —  d  for  some  d  €  F. 

From  Proposition  1,  we  know  that  for  a  field  F  containing  a  primitive  nth  root  of  unity  and  its 
extension  F[  yfa\  with  a  minimal  polynomial  xn  —  a  over  F,  K/F  is  a  cyclic  Galois  extension  of 
dimension  n.  Therefore,  there  are  a  number  of  n  embeddings  ai  =  a1  of  K  to  C  such  that  a(x)  =  x 
for  x  €  F  and  a(  tfa)  =  u  tfa.  =  f/a. 
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Definition  4  An  n-dimensional  cyclic  algebraic  space-time  code  X(K/F,(3,cr,  p)  based  on  a  cyclic 
field  extension  K/F  with  (3  =  s/a  is  a  set  of  n  x  n  matrices  with  the  form  of 


X  = 


Xi 


pn  la(xn ) 
pan~1(x2)  p2crn~l(x  3) 


px  2 
a(xi) 


Pn  2xn-i 
on~sa(xn- 2)  pn~2o(xn-i) 


Pn  lxn 


pn~1an~1(xn) 


rn—  1 


(«i) 


(5) 


where  Xi,l  —  1,  •••  ,n  are  determined  by  Xi  —  Y^k=ixk,ifik  1  with  x^i  being  algebraic  integers  in 
F.  While  p  in  (5)  can  be  any  complex  number,  in  this  paper,  it  is  chosen  from  algebraic  integers  in 
field  F. 


Notice  that  the  codes  in  [9]  and  [18]  are  generated  with  the  following  form 


xx 

X2 

Xn~  1 

xn 

pa{xn) 

a(xi) 

<?(xn-2) 

cr(xn- 1) 

<N 

1  ^ 

& 

pan~l(x  3)  ••• 

pa^ixn) 

^-'(xi) 

Although  the  codewords  in  (6)  and  (5)  differ  in  forms,  it  does  not  affect  the 
desirable  code.  The  choice  (5)  saves  the  average  transmit  energy  when  \p\ 
many  good  codes. 

Similar  to  (4),  the  cyclic  algebraic  space-time  code  defined  in  (5)  can  be 
a  complex  lattice  as  follows, 


(6) 

criterion  for  finding  a 
>  1,  as  is  the  case  for 

rewritten  as  points  of 


yi(!)  py2(i)  •••  pn  2yn- 1(1)  pn~lyn{  1) 

Pn~lVn{ 2)  2/1  (2)  •••  pn-3yn-2( 2)  pn~2yn-i( 2) 

:  :  :  :  :  (?) 

.  m(n)  p2yz{n)  •••  pn~lyn{n)  yfin) 

where  [yfi  1),  •  •  •  ,  yi{n)]T  is  a  point  in  complex  lattice  Tn(Gn>p)  with  a  generating  matrix  Gn,$  over 
Z[<n]  for  l  =  ,n, 


'  1 

1 

1 

'  1 

0  ••• 

0 

Cn 

c2 

•••  C 

0 

P  0 

0 

/•n— 1  >2(n— 1) 

L  S  n  S  n 

>n(n— 1) 

_  0 

0  ••• 

ff1-1  _ 

(8) 


Therefore,  (y1(l),y1(2),-  --  ,Vl(n),  m(l),py2(2),  -  •  •  ,  m(n),  ■■  ■  , p^ynil), Pn~lyn(2),  •  •  •  ,pn-lyn{n)\ 
can  be  considered  as  a  point  of  larger  complex  lattice  Tnn{Gn>p,  pGnfi,  •  ■  ■  ,  pn~lGn^),  with  a  gen- 
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erating  matrix  Q(Gn,p,  pGn> p,  ■■■  ,pn  1Gritp)  over  Z[(n]  x  Z[(ra]  x  •  •  •  Z[(n] 


Q{Gn^pGn,p,---  ,pn~lGn>p) 


Gn,p  0 

0  pGn,p 


0 

0 


(9) 


0  0  •••  pn-lGntf}  J 

The  absolute  value  of  the  determinant  |  det {Q{Gn^,  pGn,p ,  •  •  •  ,  Pn~lGn>p))\ of  the  generating  matrix 
G(Gn,p, pGn,fii  ■  •  ■  ip  Gn^fi)  is 


I  det {Q{Gn^pGn^  -  -  ,pn~lGn>p))\  =  |  det {Gn)\np^-^pr2^2. 


(10) 


As  a  result,  for  a  given  integer  n,  we  can  compare  two  n-dimensional  cyclic  algebraic  space-time 
codes  X(Ki/Fi,  Pi,  &i,p\)  and  X(K2/F2,  P2,  o 2,  pi)  with  the  following  criterion  by  using  the  results 
in  [21]. 

Lemma  1  For  any  two  n-dimensional  cyclic  algebraic  space-time  codes  X(Ki/Fi,  P\,oi,  pi)  and 
X(K2/¥2,  P‘2,  <72,  pi),  the  former  is  better  than  the  latter  ifdmin(K 2/F2,  P2, 02,  pi)  =  dmm(lKi/Fi,  pu  <j\,  pi) 
and  \Pipi\  <  \p2P2\,  where 

dminCK/F,  p,  a,  p)  =  min|det(X)|  (11) 

X^0,XeX(K/lF  ,/3,a,p) 

is  called  the  minimal  determinant  of  code  X(K/F,  P,  a,  p). 

In  a  special  case  where  F  =  Q(£n)  and  a  =  (n,  the  cyclotomic  space-time  code  is  a  cyclic 
algebraic  space-time  code.  In  all  the  cyclotomic  space-time  codes  with  non-vanishing  determinants, 

F  =  Q«s)  =  Q(<«),  orF  =  Q(£4).  In  this  paper,  we  only  consider  the  special  case  where  F  =  Q((s)  = 

Q(Ce),  or  F  =  Q(£4).  In  the  following,  we  will  design  optimal  cyclic  algebraic  space-time  codes  for 
three,  four,  and  six  transmit  antennas  over  F  =  Q(£3)  =  Q(Ce)>  or  F  =  Q(£4)  under  the  diversity 
product  criterion. 


Definition  5  [18]  Let  K  be  a  field  and  F  a  cyclic  extension  of  degree  d  of  ¥,  a  be  the  generator 
of  Galois  cyclic  Galois  group  Gal(K/ F).  Take  7  €  F*.  The  algebra  A  —  (K/F,  <7, 7)  generated  by  K 
and  an  element  e  is  called  a  cyclic  algebra,  where  ed  =  7,  and  e.  x  =  x .  a(e),  for  any  x  €  K.  This 
algebra  A  —  (K/F,  <7, 7)  can  be  written  as 

A  =  (K/F, <7,7)  =  K®e.I©"-0 ed~x . K.  (12) 
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This  algebra  can  be  constructed  as  a  subalgebra  of  M<j(K),  the  d-dimensional  matrix  algebra  over 
K,  by  setting 


0  1  ••  0 


0  0  ■*.  1 

7  0  •  •  •  0 


x  =  diag(al(x)),  1=0,  ■■■,d  —  1, 


From  the  definitions  of  cyclic  algebra  and  cyclic  algebraic  space-time  code,  we  know  that  a  n- 
dimensional  cyclic  algebraic  space-time  code  X(K/¥,  [5,  o,  p)  is  a  degree  n  cyclic  algebra  (K/F,  a,  pn). 
In  the  design  and  proof  of  the  optimality  of  a  cyclic  algebraic  space-time  code,  we  need  to  introduce 
the  following  lemma. 


Lemma  2  [18], [26]  A  cyclic  algebra  in  the  definition  5  is  a  division  algebra  if  and  only  if  7, 72, • -  •  ,  7^  1 
are  not  algebraic  norms  of  elements  in  K(7,72,  •  •  •  ,7d_1). 


3.1  Optimal  Cyclic  Algebraic  Space-time  Code  for  Three  Transmit  Antennas 

In  this  subsection,  we  assume  that  F  =  Q((3)  and  3K  =  ¥(p)  is  a  cyclic  field  extension  of  degree  3 
for  some  p  with  P3  =  a  €  Z[£3]. 

Now  let  us  consider  field  extension  towers  Q  C  Q(£3)  C  Q((3,P)-  Then,  any  2  e  Z[£3,  0\  can 
be  written  as  z  ;=  Z\  +  Z2P  +  Z3P 2  €  Q(£3 ,/?)  =  F(/3),  where  Zk  e  Z[£3],  fc  =  1,2,3.  Notice  that 
nQ(6^)/Q(6)(z)  “  IlLi  ak(z)>  where  a k,k  =  1,2,3  are  the  three  embedding  of  Q(£3 ,P)  to  C  that 
is  fixed  over  Q(£3)  with  <Tk{P)  =  (Jf~xp.  Hence,  from  the  definition  of  the  relative  algebraic  norm 
we  have 

—  z\  +  ^  z2  +  f^zz  +  3/93(C3  +  d)ziZ2Z3 

=  zl+p3z%+p6zl-3P3z1z2z3.  (14) 

If  we  let  p  =  21/3,  then,  we  know  from  the  algebraic  number  theory  that  Q(£3, 21/3)/Q((3)  is  a 
cyclic  field  extension  of  degree  3.  In  the  following  we  will  prove  that  (3  and  (e  are  not  the  algebraic 
norm  of  some  integer  of  Q(£3, 21/3)  over  Q((3). 


Theorem  1  Foranyx,y  €  Z[<3 ,2x/3\,  i/NQ(0,f2i/3)/Q(&)(ar)  =  C6%C3)2i yij/Q^Cy),  Nq^i/s )/q(<,)(®) 
^^QfCs^/si/QfCs)^)'  or  NQ(C3,21/3)/Q(C3)(a;)  =  ^3NQ(Cs ,2V3)/q(c3)  (y)  tflen>  x  =  V  =  °- 
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The  proof  of  Theorem  1  is  given  in  Appendix. 

Similarly,  we  can  prove  that  2  is  not  the  algebraic  norm  of  some  integer  of  Q(C3,Ce)  =  Q(Cis) 
over  Q(C3). 

Theorem  2  Foranyx,y  €  Z[Cis]  (x,y  €  Z[<g]),  */Nq(Ci8)/q(6)(*)  =  2N q(c18)/q(c6)(v)  (nQ(6)/QKs)(®)  = 
2NQ(C9)/Q(C3)(y)),  then,  x  =  y  =  0. 

For  any  cyclic  algebraic  space-time  codeword  X  €  X(Q((3,  P)/Q((s),  (3,  <J,  p)  with  the  following 
form 

x  py  p2z 

X  =  pV(z)  a(x)  pa{y)  ,  (15) 

_  pa\y)  p2cr2(z )  cr2(a:)  _ 

where  a:  =  V  =  Z)f=i  V*/5*-1.  z  =  Z)f=i  ^l-1>  with  %uVuzi  e  ^tCa] ,  /33  €  Z[C3],p3  G 

Z[£3],  we  can  prove  that  its  determinant  belongs  to  Z[£3];  i.e., 

Theorem  3  T/ie  determinant  det(X)  of  matrix  X  with  the  form  of  (15)  belongs  to  Z[£3],  i.e., 

det(X)  €  Z[<3].  (16) 

The  proof  of  Theorem  3  is  provided  in  Appendix. 

The  following  theorem  guarantees  that  cyclic  algebraic  space-time  codes  X(<Q>(£3, 2  |/3)/Q((3),  21/3,  a ,  (g), 
^(Q(C3)21/3)/Q(C3),21/3,o-, <«),  A(Q(C9)/Q(C3)!<9>cr,21/3))  and  X(Q(Ci8)/Q(C«), <18, <r, 21/3)  for 
three  transmit  antennas  provide  full  rate  full  diversity  with  non-vanishing  determinants. 

Theorem  4  The  cyclic  algebraic  space-time  code  X(Q((z,  21/3)/Q(^3),  21/3,  o,  (9), 

^(Q«3,  21/8)/Q(Cs),  21/3,  a,  (is),  X(Q(C9)/Q(C3),C9,a,21/3))  and  X(Q{Sis)/Q{C*),<i&,<r,  21/3),  are 
full  diversity  space-time  code  with  the  minimum  determinant  1. 

Proof  of  Theorem  4:  From  Theorem  3  we  know  that  the  minimal  determinants  of  these  cyclic 
algebraic  space-time  codes  are  either  1  or  0.  Therefore,  we  only  need  to  prove  that  these  cyclic 
algebraic  space-time  codes  enable  full  diversity,  i.e.,  if  det(X)  —  0,  then,  X  =  0,  since  they  are 
linear.  In  addition,  from  Theorem  1,  Theorem  2  and  Lemma  2,  we  know  that  any  non-zero  codeword 
X  in  the  algebraic  space-time  codes  of  Theorem  4  does  not  have  a  zero-divisor,  i.e.,  if  det(X)  =  0, 
then,  X  =  0. 
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Q.E.D. 

Furthermore,  we  can  prove  that  the  resulting  cyclic  algebraic  space-time  codes  for  three  transmit 
anetnnas  in  Theorem  4  is  optimal  based  on  the  criterion  given  in  Lemma  1. 

Theorem  5  Under  the  diversity  product  criterion,  X(Q((s,  21/3) /Q(£3),  21/3,  o,  (g),  X(Q(£3,  21/3) /Q((3), 
21/3,<r,<is),  X(Q(C9)/Q(C3),C9,o',21/3)  andX(Q«i8)/Q(^6),^18,a>21/3)  are  the  optimal  cyclic  space- 
time  codes  for  three  transmit  antennas  with  the  minimal  determinant  1. 

The  proof  of  Theorem  5  is  given  in  Appendix.  From  Lemma  1  we  know  that  diversity  product 
of  a  cyclic  algebraic  space-time  code  is  determined  by  \ftp\,  the  absolute  value  of  product  ft  and  p. 
Therefore,  the  main  idea  to  prove  the  optimality  of  these  codes  is  to  prove  ft3  =  2  and  p3  =  £3  are 
integers  in  Z[(3j  with  the  minimum  product  \ftp\  such  that  p,  p2  are  not  relative  algebraic  norm  of 
any  elements  of  Q((z,ft)  over  Q((3). 

Remarks:  The  cyclic  algebraic  code  ^(Q(C9)/Q((3),  Co,  <?,  21/3)  is  a  cyclotomie  space-time 
code,  but  AT(Q(^3,21/3)/Q((3),21/3,<t,  £9)  is  not.  Although  they  have  the  same  diversity  product, 
i.e.,  the  same  performance  under  the  diversity  product  criterion,  the  former  with  p  =  21/3  >  1  has 
different  average  power  at  each  layer  and  the  latter  with  |p|  =  |£9|  =  1  has  the  same  average  power 
at  each  layer.  Therefore,  the  latter  has  a  lower  peak  to  average  power  ratio. 

3.2  Cyclic  Algebraic  Space-Time  Code  for  Four  and  Six  Transmit  Antennas 

In  this  subsection,  we  will  generate  cyclic  algebraic  codes  for  four  and  six  transmit  antennas. 

Let  a  —  ft4  €  Z[£4],  ft,  ft2,  ft3  £  Q«4),  Q(C4,/?)/Q(C4)  be  a  4-dimensional  cyclic  field  extension; 
he.,  [Q(C4,/3)  :  Q(C4)]  =  4  and  l,  ft,  ft2,  ft3  is  a  base  of  Q(C4,  ft)  over  Q(C4)-  Let  a  be  the  generator 
of  the  Galois  group  of  a  field  extension  Q(Q,  ft) /Q(U)  with  ak(ft)  =  <4_1/?,fc  =  1,2, 3, 4.  For  any 
element  x  =  X\  +  Xgft  +  xgft'2  4-  x^ft3  with  Xk  €  Z[£4],  its  relative  algebraic  norm  NQ(c4,p)/Q(c4)(a:) 
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can  be  expressed  as 


where 


4 

N®(C4,p)m<4)(x)  = 

k—  1 

=  —x\  +  /34(x|  +  2*1*3  +  4xjx2X4  -  4*i*2*3) 

+  /?8(4x2*3*4  -  4x1*3*4  -  2x\x\  -  a:4)  +  (312x\ 


—  —X4  +  a(x4  +  2*f*3  +  4x(*2*4  -  AxiX^Xz) 

+  a2(4x2x2x4  -  43:1*3X4  -  2*2*4  ~  *3)  +  a3*4 
=  -*1  +  a/(*lt*2,*3,*4,Q!)  e  Z[<4], 


(17) 


/(*1,*2,*3,*4,«)  =  (*2  +  2*1*3  +  4*j*2*4  “  4*i*2*s) 

+  ft(4*2*3*4  -  4*i*3*4  -  2*2*4  -  *3)  +  a2*4  £  Z[<4].  (18) 

Similar  to  Theorem  4  for  the  three  transmit  antennas  case,  we  have  the  following  non- vanishing 
determinant  theorem  for  four  transmit  antennas. 


Theorem  6  X(Q((4, P)/Q(Ci),P,  cr,  C16)  where  fi  =  (2  -±  j)1/4  is  a  full  rate,  full  diversity  non¬ 
vanishing  determinant  code  with  the  minimum  determinant  1  and  its  average  power  at  each  layer 
is  the  same. 

The  proof  of  Theorem  6  consists  of  two  parts.  In  the  first  part,  we  first  prove  that  j,  -j,  and  -1 
can  not  be  a  algebraic  norm  of  any  element  of  Q(( 4,  /?)  over  Q«4).  Then,  by  Lemma  2,  we  know  that 
^(<□(<4,  (x,  Cie)  is  a  full  diversity  space-time  code.  This  part  constitutes  Lemma  3  in  the 

Appendix.  In  the  second  part,  we  prove  that,  for  any  codeword  X  £  -XXQ(<4,/?)/Q(<4),/?,a,<16), 
its  determinant  det(X)  €  Z[£4],  which  forms  Lemma  4  in  Appendix. 

Finally,  combining  Lemma  3  and  Lemma  4,  we  prove  that  X(Q((4,  P)/Q((a),  A  <x,  (i6)  is  a  full 
diversity  space-time  code  with  the  minimal  determinant  1.  Since  |p|  =  |(16|  =  1,  the  average  power 
of  this  code  at  every  layer  is  the  same.  The  detail  of  the  proof  is  given  in  Appendix. 

Analogously,  we  obtain  the  following  theorem. 

Theorem  7  X(Q(C4,0)/Q(Q),  with  (3  =  (2  +  j)1/4  is  an  optimal  cyclic  algebraic  space- 

time  code  with  the  same  average  power  at  each  layer. 
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The  proof  of  Theorem  7  is  similar  to  that  of  Theorem  5.  Here  we  only  give  the  outline  of  its 
proof.  In  order  to  prove  this  theorem,  it  suffices  to  verify  that  /34  =  2  +j  and  p4  =  j  are  the  integers 
in  Z[Ci]  with  the  minimal  value  \(3p\  such  that  p,/?2,p3  are  not  the  relative  norm  of  any  elements 
of  Q(Q,0)  over  Q«4). 

Similarly,  we  have  the  following  two  theorems  for  six  transmit  antennas. 

Theorem  8  The  cyclic  algebraic  code  X(Q(Q,/3)/Q((6),  (2  +  (e)1/6,  cr,  (36)  is  a  full  diversity,  full 
rate,  determinant  non-vanishing  code  with  the  minimal  determinant  1  and  the  average  power  at 
each  layer  is  the  same. 

Theorem  9  X (Q(£6,  /?) /Q((e) ,  /3,  <x,  (m)  with  (3  =  (2  +  (g)1/6  is  an  optimal  cyclic  algebraic  space- 
time  code  with  the  same  average  power  at  each  layer.  . 

4  Simulation  Results 

In  this  section  we  consider  a  MIMO  system  with  three  transmit  antennas  and  two  receive  antennas. 
We  compare  the  error  performance  of  the  space-time  code  A'(Q(C9)/Q(C3)>  (9, 0-,  (3+exp(j7r/,3))1/3) 
proposed  in  [18]  with  our  proposed  code  X(Q((z, 21/3)/<Q>((3),  21/3, a, ($).  While  both  have  non¬ 
vanishing  determinants,  the  latter  has  a  lower  peak  to  average  power  ratio  than  the  former,  since 
|p|  =  |(3  +  exp(j7r/3))1/,3|  >  1  in  the  former  whereas  |p|  =  [Co]  =  1  in  the  latter  code.  Moreover, 
according  to  the  diversity  product  criterion,  the  proposed  code  also  has  better  codeword  error 
performance.  Our  simulation  results  depicted  in  Fig.  1  show  that  our  new  code  achieveds  about 
2dB  gains  over  the  one  proposed  in  [18]. 

5  Conclusion 

A  systemic  non-vanishing  determinant  space-time  code  design  has  been  proposed  using  some  cyclic 
field  extensions.  Based  on  the  diversity  product  criterion,  the  optimal  space-time  codes  for  three, 
four  and  six  transmit  antennas  have  been  obtained.  We  have  proved  that  these  optimal  codes  have 
the  same  average  power  at  each  layer  and,  subsequently  a  low  peak  to  average  power  ratio. 
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Appendix 

Proof  of  Theorem  1 

Here  we  only  prove  the  first  part  of  Theorem  1.  The  second  part  can  be  proved  in  the  similar  way. 
Let  x  =  xi  +  /3x2  +/32x3,  y  =  yx  +/%2  +P2yz,  with  p  =  21/3,  xk,  yk  €  Z[C3],  k  =  1,2, 3.  Since  2Z([<3] 
is  a  ideal  of  ring  Z[(3],  for  any  given  x  and  y,i.e.,  xk,yk  G  Z[£3],  x  and  y  can  be  written  as  follows 
with  some  integer  Iq 

lg  Iq 

xk  =  ^2  2  l~lxkti,  yk  =  5^2  l~lyk,h  (19) 

i=i  i=i 
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where  Xk,i,yk,i  €  {0,exp(j'2p7r/6),  \/3exp(j7r/6)  exp(j2px/6),p  =  l,  --  ,6}  for  k  =  1,2,3  and  l  = 
1,  •  •  •  ,  lo-  If  NQ(Cc3,a)/Q(C3)(x)  ■=  C6NQ(Cc3]Q)/Q(C3)(y),  then,  from  (14)  we  have  the  following  equation 

y\  -  (exf  =  2  {  (x\  +  2x\  -  3xix2x3)  -  Ce  (yl  +  2 y|  -  Zymyz)  }  ,  (20) 

that  is 

Via  ~  Ce^li  =  2  { (xf  +  2xl  ~  3xix2x3)  -  Ce  (j/|  +  2vl  ~  Zymyz) } 

+  2  {  (3®?,!*!,!  +  6x1,1®?  !  +  4x?;1)  -  Ce  (3y?,i?/i,i  +  63/1,1^1  +  4yf;X) }  ,  (21) 

where  xi,i  =  Sp=2  2p-2xj,p,  1/1,1  =  Sp=22p_22/i,p.  The  term  on  the  right  hand  side  of  (21)  belongs 
to  2Z[£3],  so  does  the  term  on  the  left  side  hand  side  of  (21);  i.e., 

x\,i  ~  Cel/u  e  2Z[£3].  (22) 

Checking  (22)  with 

*i,i.Vi,i  €  { 0 , exp(j 2jm / 6) ,  \/3 exp(  j?r/6) exp (J 2pn /6) , p  =  1,...,6},  (23) 

we  can  find  that  it  is  true  only  when  Xi,i  =  j/1,1  =  0.  In  this  case,  equation  (21)  becomes 

y\  -  Ce^2  =  2{(x|  +  2x?il-3xi,ix2x3) -C6(yf +  2y?,i-3yi,iy2y3)},  (24) 

that  is 

yi,i-C«*a,i  =  2  {  (x|  +  2x?  i  —  3xi,ix2x3)  —  (y3  +  2y|,i  —  SyiAV^Vz)  } 

+  2  {  (3x|ix2,i  +  6x2,ix?  i  +  4x|x)  -  <6  (3y|, 1^2,1  +  63/2,11/2,1  +  4y|i) }  ,  (25) 

where  x2,i  =  Y?p=- 2  2p~2x2,p,  y2,i  =  Sp=2  2p_22/2,p-  The  term  on  the  right  hand  side  of  (25)  belongs 
to  2Z[£3],  so  does  the  term  on  the  left  hand  side  of  (25),  i.e., 

i4,i-&4i€2Z[<3].  (26) 

Similar  to  the  proof  for  X14  =0  and  3/1,1  =  0,  we  can  obtain  that 

£2,1  =  2/2,1  =  0.  (27) 

and  furthermore,  ®i,2  =  yi,2  =  x3,i  =  y3,i  =  0.  Finally,  we  have  x  =  y  =  0. 

Q.E.D. 
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Proof  of  Theorem  3 


It  can  be  directly  verified  by  computation  that 

det(X)  =  a(x)(r2(x)a3(x)  +  p3a(y)a2{y)a3(y)  +  p6a(z)a2(z)a3(z) 

~  PS(^(x)(r2{y)(T3{z)+a(y)a2(z)a3(x)  +  (x(z)a2{x)(T3(y)).  (28) 

Since  a(x)cr2(x)cr3(x)  =  NQ(C3i/3)/Q(0l)(a;)  €  Z[<3],  o{y)a2{y)a3{y)  =  Nq(C3)/3)/q (&)(y)  e  Z[Cs], 
o-(2:)(T2(z)o-3(2:)  =  e  Z[£3]  and  p3,p6  6  Z[£3],  the  first  three  terms  of  (28)  be¬ 

long  to  Z[C3J.  Now  we  only  need  to  prove  that  the  coefficient  of  p3  in  (28)  also  belongs  Z[(3], 
i.e.,o-(a;)<72(?/)0-3(z)  +  a(y)a2(z)cr3(x)  +  cr(z)a2(x)a3(y)  €  Z[<3] .  Notice  that  a(x)a2(y)a3(z)  + 
a(y)a2(z)a3(x)  +  a{z)a2{x)a3(y)  in  (28)  can  be  rewritten  as 

a{x)a2{y)a3{z)  +  a(y)a2(z)a3{x)  +  a{z)a2(x)az(y) 

E  Ph+h+h-\iyhzh  (tl>+2h  +  C33+2Zl  +  C31+2'2)  •  (29) 

h,h,h€{  1>2)3} 

If  h  +  l2  +  k  is  divided  by  3,  then,  Ph+h+h~3xhyhzh  (<f?+2*8  +  <£+afl  +  C3l+2/2)  €  Z[(3].  If 
h  +  h  +  h  is  not  divided  by  3,  then,  li  -f  2/3,  ls  +  2 ly  and  l\  +  2l2  are  not  congruent  to  each  other 
under  modulo  3.  Therefore,  they  run  through  a  complete  residue  system  of  3.  As  a  result,  C^+2h  + 
C33+2Zl  +  C31+2Za  =  1  +  (3  +  C3  =  0-  Therefore,  a(x)a2(y)a3(z)  +  cr(y)cr2(z)cr3(a;)  +  a(z)a2(x)a3(y) 
and  hence,  det(X)  belongs  to  Z[(3], 

Q.E.D. 

Proof  of  Theorem  5 

First  we  prove  that  ^(Q^1/3)/^)^1/3,^),  A(Q«3, 21/3)/Q«3), 21/3,  <7,  Cis)  and 
^f(Q(Ci8)/Q(C6),Ci8)Cr,21/3)  have  the  same  average  powers.  Since 

px(Q(C3,2i/3)/Q(c3),2i/3)<7iC9)  =  -px(Q(C3,2V3)/Q(C3)i2i/3,<r)Cl8)  =  12  det(G3,3)|3  =  23j  det(G3>3)j3,  (30) 

P^(Q(Ci8)/Q(C6),Ci8,<T,2V3)  =  2  x  22) det(G3)3)|3  =  23|  det(G3>3)|3,  (31) 

and  the  minimum  determinants  of  cyclic  algebraic  space-time  codes  X(Q((3, 21/3)  /Q(C3),  21/3,  a,  (9), 
•^(Q(C3)21/,3)/Q(C3),21'/3,cr,  Cis)  and  A"(Q(C18)/Q(^6),  Ci8)  <r,  21/3)  are  1,  therefore,  they  have  the 
same  performance  based  on  the  diversity  product  criterion. 


80 


In  the  following,  we  will  prove  that  any  other  cyclic  algebraic  space-time  code  over  Q(£3)  defined 
in  this  paper  has  a  less  diversity  product  than  X(Q(£3, 21//3)/Q(C3),  21/3,  a,  (g)  and  ^(Q(Cis)/Q(C6), 
(is,  ex,  21/3)  have.  From  the  definition  of  cyclic  algebraic  space-time  codes,  we  know  that  ft3  G 
Z[(e]  =  Z[(3].  Now  we  consider  the  following  four  cases  according  to  the  different  values  of  ft: 

1-  ft2  €  {  exp(;2p7r/6), p  =  1,  •  •  •  , 6}. 

2.  ft 3  €  { \/3Ci2 exp(j 2p7r /6) ,  p  =  l,-  -  ,6}. 

3.  ft 3  G  {2exp(j'2p7r/6),p  =  l,--  -  ,6}. 

4.  |/?3|  >  2. 

Case  1:  /33  G  {  exp(y 2p7r/6) ,  p  =  1,  *  -  •  ,6}.  If  ft3  =  1,  or  ft 3  =  -1,  then,  ft  G  Z«6)]  C  QCCs)- 
So  <Q>(/3, Ce)  =  Q(Ce),-  In  this  case,  ft  can  not  be  used  for  the  cyclic  space-time  code  design.  If 
fto  =  <6,<6><6  or  Cl>  then.  for  any  |p3 1  <  2  with  p3  G  Z[C3],  X(Q((3,fto)/Q((3 ),ft,<r,p)  can  not 
generate  any  full  rate  full  diversity  space-time  code.  It  is  not  difficulty  to  find  that  every  element  in 
the  set  S  =  { exp(j2pir/6),  \/3exp(j7r/6) exp(j2p7r/6), p  =  I,---  ,6}  is  an  algebraic  norm  of  some 
elements  in  Qfe,fto)  over  Q((e)-  From  Lemma  2,  we  know  that  X (Q((3,  fto / CMC3) >  (ft  p)  cannot 
generate  any  full  rate  space-time  code  for  p3  G  S.  If  choose  p  such  that  |p|3  >  2,  the  average  power 
of  X(Q«3,  /?o)/<Q«3),  ft,  <r,  p)  satisfies 

P*(Q(c3,Pb)/Q(C3 ),/W)  =  I/5!9!  det(G3j3)|3  >  23|  det(C3,3)]3  =  <?,&)■  (32) 

Case  2:  ft3  g  {  V3Ci2 exp(j2p7r/6), p  =  1, •••  ,6}.  In  this  case,  it  is  not  hard  to  find  that 
exp(j2p7r/6),p  =  1,  •••  ,6  are  algebraic  norms  of  some  elements  in  Q(&,ft)  over  Q(Ce)-  To  gen¬ 
erate  a  full  rate  cyclic  algebraic  space-time  code  X(Q((3,  ft)/Q(&),  ft,  cr,  p),  p  can  not  be  cho¬ 
sen  from  set  exp(j2p7r/6),p  =  1,  •••,6.  If  choose  p  with  |p|3  >  -y/3,  the  average  power  of 
*(Qtts,  ft)/Q(Ca),  ft,  <rt p)  satisfies 

-pa:(q(C3,P)/Q(C3),P,M  =  |p|9j/?|9|det(G3>3)|3  >  ^1  det(G3;3)|3 

>  23|det(G3;3)|3  =  -pX(q(C9)/q(C3))C9j<7i2i/3).  (33) 

So  it  can  not  achieve  optimal  cyclic  algebraic  space-time  code  in  this  case. 
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Case  3:  /?3  €  {2exp(j2p7r/6),p  =  1,  •  •  *  ,6}.  In  this  case,  by  Theorem  1,  £6,  (3  and  (f  are  not 
the  algebraic  norms  of  any  element  of  Q((3, /3)  over  Q((3) .  Therefore,  X(Q«3, 21/3) /Q«3),21/3,  a,  p ) 
is  a  full  diversity  space-time  code  with  diversity  product  1  with  p3  =  (3  or  p3  =  (6.  Therefore,  these 
pare  the  numbers  with  smallest  absolute  value  that  make  <Q>(£3,  j3)  over  Q((3).  X(Q((3, 21/3) /Q((3),  21/3,  a,  p) 
has  minimum  determinate  1. 

Case  4:  |/?3|  >  2.  In  this  case,  for  any  p  >  0  with  p3  e  Z[£3],  |p|  >  1. 

-Px'(q(C3,/3)/<Q(C3),^,ct,p)  —  l/^l  IpI  | det(C3i3)|  >2  |det(G3)3)|  =  -fx(Q(C3,21/3)/Q(C3),21/3i0->^9)-  (34) 

So  the  optimal  cyclic  algebraic  space-time  code  can  not  be  obtained  in  this  case.  This  completes 
the  proof  of  Theorem. 

Q.E.D. 


Proof  of  Theorem  6 

In  order  to  prove  Theorem  6,  we  first  establish  the  following  Lemma  3  and  Lemma  4. 

Lemma  3  For  any  two  algebraic  integers  x  andy  o/Q( 6,0),  */NQ(C4>j9)/Q{C4)(ar)  =  C^q 

for  k  =  1, 2  or  3,  then  x  =  y  =  0,  i.e.,  j,  — j,  and  —1  can  not  be  a  algebraic  norm  of  any  element  of 

Q(C4 ,P)  over  <Q>(£4),  where  (3  =  (2  ±  j)1^. 

Proof.  Here  we  only  give  the  proof  in  the  case  where  (3  =  ax/A  =  (2  +  j)1/4  and  k  =  2.  The 
other  cases  can  be  proved  in  a  similar  way. 

For  any  given  x,  y  with  a:  =  xx  +  x20  +  x3(32  +  x4(33,  y  =  yi  +  y2(3  +  y3(32  +  y4(33,  xm,  yrn  € 
^[(4], m  =  1, 2, 3, 4,  because  (2  +  i) Z{£4]  is  a  ideal  of  ring  Z[£4],  there  is  an  integer  Iq  such  that 

*  =  ^x^a1-1  +  pJ2x2,ial~1  +  02J2x3jial~1  +  f 

i=l  z=l  1=1  1=1 

Jo  Iq  Iq  Iq 

y  =  +p2Y^y^al~x +/33J2v4,ia1-1,  (35) 

1= 1  i=i  1= 1  j=i 

where  xm>i, ym>i  €  S  =  {0, 1,-1,  j, - j ,  1+j,1-j,  -1+j,  -l-j,2,  -2,2j,-2j}.  = 

-NQ(C4,/3)/<Q(C4 )(y)>  by  (17)  and  (18),  we  have 

®1+»1  =  -a(f(xuX2,x3,x4,a)  +  f(yi,y2,y3,y4,a)).  (36) 
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Substituting  (35)  into  (36)  yields 

«f,i  +  2/i,i  =  -at(f(zi,X2,X3,x4,a)  +  /(yi,2/2,2/3,2/4,a))  +  a(g(xltziti)  +  g(yi,yi,i)),  (37) 
where 


y(*ii*i,i)  =  a  -xftl)  e  Z[<4],  (38) 

5(2/1, 2/1, 1)  =  a~\yi  -vii)  €  z[(4]-  (39) 

Then  the  term  on  the  right  hand  side  of  equation  (37)  belongs  to  aZ[£4],  so  does  the  term  on  the 
left  hand  side  of  equation  (37),  i.e., 


®i,i  +2/1,1  €  IC43  - 


(40) 


Checking  (40)  for  xi,i,  2/1,1  e  S,  we  can  find  that  (40)  is  true  only  when  xiti  —  1/1,1  —  0.  In  this 
case, 

*  =  p  (y!X2-io/1~1 + p!l,x  +02^x4,ial~1 +0Z^xi,ioi1~1'\  =  0x, 

\l=l  1=1  1=1  1=2  J 

y  =  0  f^2/2liai_1  +  P^2,yz,ial~l  +  2/4,1^“ 1  +  03  ^yi,ial~l\  =  0y,  (41) 

\i=i  1=1  1=1  1=2  / 

and 


NQ«4,/J)/Q(C4)(a:)  =  /^NQ(C4,0)/Q(C4)(aO  =  aNQ(C4,/3)/Q(C4)(®)  =  a(4  +  /(®2,  £3,  #4,  aq),  a), 
NQ(C4,/3)/Q(C4)(?/)  =  /^NQ(C4,/?)/Q(C4)(2/)  =  aNQ(C4,/?)/Q(C4)(^)  =  a(V2  +  f(y 2. 2/3, 2/4,  2/l),  a), 

(42) 

where 

xi  =  Y1  xi,iPl~2’  vi=Yl 

1=2  1=2 

x  =  J2xva1-1  +  0'}Tx3tial~1  +/?2^*4,iaz_1  +/?3  J^aq^a*"1, 
i=i  1=1  i=i  1=2 

/o  Jo  Jo  Jo 

y  =  V2,ial~1  +&Y1  ?/3,iaZ“ 1  +  (32  Vu®l~l  +  P3J2  V^~X  -  (43) 

1=1  1=1  1=1  1=2 

With  the  above  notation,  equation  (37)  becomes 

4,1  +  2/2,1  =  -a(/(«2,®3,*4,ai,a)'+/(lte,2te,y4,yi,a))  +  a0/(*2,*2,i)  +  5(512,2/2,1))-  (44) 
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The  term  on  the  right  hand  side  of  the  equation  (44)  belongs  to  a:Z[<(4],  so  does  the  term  on  the 
left  hand  side  of  equation  (44),  i.e., 

*ti+l4i€aZ[Ct].  (45) 

Following  a  similar  step  to  prove  aq,i  =  yi,i  =  0  in  equation  (40),  we  can  get  X2,i  =  j/2,1  =  0,  and 
furthermore,  0:3,1  =  y3j  1  =  0,  a:4ii  =  y4)i  =  0,  aqi2  =  S/i,2  =  0,  •  •  •  ,  xm,i0  =  ym,i0  =  0-  Hence,  we  get 
a:  =  y  =  0. 

Q.E.D. 


Lemma  4  For  any  codeword  X  e  X(Q(C4,/3)/<Q>(<4),/?,<7,  Cis),  det(X)  €  Z[<4], 


Proof.  For  4x4  matrix 


Xx 


its  determinant  det(Xi)  is 


ai  b\  C\  di 
d2  a2  62  c2 
c3  d3  a3  63 
64  c4  d4  a4 


(46) 


det(Xi)  =  aia2a3a4  —  61 6263^4  +  cic2c3c4  —  did2d3d4 

—  (aia263d4  +  ai&2d3a4  +  &id2a3a4  +  dia2a364) 

+  (01^2^3^  +  6ic2a364  +  cia26364  +  f>i62C3a4) 

+  (cid2d3a4  +  did2a3c4  +  dia2c3d4  +  aic2d3«;4) 

—  (cid263c4  +  cic2d364  -f  &ic2c3d4  +  di&2c3c4) 

+  (6id263d4  +  di62d364)  —  (aic2a3c4  +  cia2c3a4). 

So,  for  a  codeword  X  €  X(Q(£4,/?)/Q(<(),/?,  c,  p)  with  the  form  of 

a:  py  p2z  p3w 

X  _  pM™)  <r(s)  po-(y)  p2cr(z) 

p2<r2(z)  p3a2(w )  a2  (a;)  pcr2{y) 

pv3(y)  p2cr3(z)  p3a3(w)  a3(w) 


(47) 


(48) 
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we  have 


det(X)  =  xa(x)a2(x)cr3(x)  -  p4y<x(y)a2(y)a3(y)  +  p8za(z)a2(z)a3(z)  -  p12wa(w)a2(w)a3(w) 

-  p4(xa(x)a2(y)a3(w)  +  xa(y)a2(w)a3(x)  +  ya(w)a2(x)a:i(x)  +  wa(x)a2(x)a3(y)) 

+  P4{%<r(y)<r2(y)<73(z)  +  ya(y)a2(z)<j3(x)  +  ya(z)a2(x)a3(y)  +  zcr(x)a2(y)a3(y)) 

+  p8(za(w)a2(w)cr3(x)  +  wa{w)a2{x)az(z)  +  wcr(x)a2(z)a3(w)  +  xa(z)a2(w)a3(w)) 

-  p8{za(w)a2(y)a3(z)  +  wa(y)a2(z)a3(z)  +  ya(z)a2(z)a3(w)) 

+  p8(ya(w)a2(y)a3(w)  +  wa(y)a2(w)a3(y)) 

-  p4(xa(z)a2(x)a3(z) +  za(x)a2(z)a3(x)).  (49) 


Since  the  symmetry  of  the  elements  in  (49),  and  xa(x)a2(x)a3(x)  =  ^q^){x)  G  Z[(4],  pA  e  Z[<4], 
we  only  need  to  prove  that 


xa(x)a2(y)a3(w)  +  xa(y)a2(w)a3{x)  +  ya(w)a2(x)a3(x)  +wa(x)a2(x)a3(y)  G  Z[(4] 


and 

xa(z)a2(x)a3(z)  +  za(x)a2(z)a3(x)  G  Z[(4). 

Let  X  =  EL  1  Xmr~\  y  =  ELl  Vmf?n-\  Z  =  ELl  *™r-\  w  =  ELl  Wmr-1. 

xa(x)a2(y)a3(w)  +  xa(y)a2(w)cr3{x)  +  ya(w)a2(x)a3(x) +wa(x)a2(x)a3(y) 

Y  xm1xrn2ym3zrni{i6mia(pTn2)a2(l3m3)a3(Pmi)  +  ft"12  a(pm3)cr2  (pmi)a3  ([3mi) 
+0m3a(pmA)a2(0mi)a3(/3m2)  +  pm4a{pmi)a2(pm2)a3(pm3)} 

Y  xm1Xm2ym3Zmi  {j3mi+^+m3+mi-4 ^m2+2m3+3mi  +  J-m3+2m4+3mi 


mi  ,77^2,7713,7^4 


mi  ,m2  >7773,7724 


_)_j>n4+2roi+3m2  _j_  jtoi+2to2+3to3^  j 


(50) 


where,  1  <  mi,  m2,  m3,  m4  <  4.  It  is  easy  to  verify  that  jm2 +2rn3  +3m4 _|_ jm3 +2m4 +3mi  _^j7i4+2mi+377i2_|_ 
J-mi+2m2+3m3  =  0  if  mi  +  m2  +  m3  +  m4  is  not  divided  by  4.  Therefore,  a(x)a2(y)cr3(w)  + 
xa(y)a2(w)a3(x)  +  ya(w)a2(x)a3(x)  +  wa(x)a2(x)cr3(y)  g  Z[<4]  and  as  a  result,  det(X)  G  Z[(4], 

Q.E.D. 
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Proof  of  Theorem  8 


The  proof  of  Theorem  8  can  be  done  by  establishing  the  following  Lemma  5  and  Lemma  6  and 
using  Lemma  2  with  |£36|  =  1. 

Lemma  5  For  any  two  algebraic  integers x  andy  o/Q«6,/?),  */%(<«, /3)/Q(C6)(>)  =  P)M&)(y)’ 

fork  =  1,2, 3, 4,  or  5,  thenx  =  y  =  0,  i.e.,  can  not  be  a  algebraic  norm  of  any  element  of  Q((e,(3) 
over  <Q(Ce),  where  @  =  (2  +  C6)1//6- 

Proof.  Similar  to  the  proof  of  lemma  3  for  four  transmitter  antennas  and  using  the  fact  that 
x  =  y  =  0,  if  x6  -  < 6V  €  (2  +  <6)Z(<6),  for  x,y  e  {0, <?,  s/Z^l 2$ ,p  =  1,  •  •  •  ,  6}. 

Q.E.D. 

Lemma  6  For  any  codeword  X  e  X(Q((e,  f3)/Q((e),  (2  +  C6)1/6,er,<36),  then  det(X)  e  Z[C6], 

J 

Proof.  Similarly,  the  proof  of  this  lemma  can  be  done  by  using  the  fact  that  fi+h=0  if 
m\  +  m2  +  m3  +  TO4  +  ms  +  me  is  not  divided  by  6,  where 

f  1  _  £Tn2+2m3+3m4+4m5+5m6  _j_  ^13 +2m4  +3m5  +4m6 +5mi  _j_  +2m5  +3m{  +4mi  +5rri2 
f2  =  ^l5+2?ri6+3mi4-47n2+5m3  _|_  ^m6+2mi+3m2+4m3+5rri4  _j_  ^miH-2m2+3m3+4m4+5m5 

Q.E.D. 
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Codeword  Error  Rate 


CER  of  Full  Rate  Cyclotomic  Space  Time  Code  for  Lt=3,  Lr=2,  Rate=6  bits 


Figure  1:  Codeword  Error  Probabilities  of  Cyclic  Algebraic  Space-Time  Code  and  Cyclotomic 
Space-Time  Code  with  Three  Transmit  and  Two  Receive  Antennas. 
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Networks* 

Gjergji  Kurti1,  Yimin  Zhang1-,  Morgan  Watson1,  and  Moeness  Amin1 

1  Center  for  Advanced  Communications 
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Abstract 

Unmanned  surface  vehicles  (USVs)  are  envisioned  to  be  a  multi-mission,  rapidly  config¬ 
urable,  long  endurance  unmanned  platform  providing  offboard  sensor  and  weapon  capability 
to  perform  critical  missions  like  intelligence,  surveillance,  reconnaissance,  anti-submarine 
warfare,  mine  countermeasures,  communications,  and  navigation.  In  this  paper,  we  consider 
the  propagation  environment  and  the  related  technology  of  multi-sensor  communications  in 
unmanned  surface  vehicle  networks  located  in  a  rough  sea  environment.  We  first  investigate 
the  statistical  model  of  the  propagation  environment,  where  we  provide  the  analytical  expres¬ 
sions  of  the  wave  height,  the  probability  of  having  line-of-sight  between  two  USVs  in  terms 
of  the  weather  condition  and  antenna  height,  and  the  spatial  correlation  of  the  propagation 
channels  related  to  different  antennas.  We  then  investigate  the  fading  reduction  capability 
of  a  multiple-input-multiple-output  (MIMO)  system  and  the  performance  in  general  Ricean 
channels  with  different  channel  correlations  is  examined.  The  results  are  important  to  the 
understanding  of  the  communication  environment  where  USV  communications  take  place. 


‘This  work  was  supported  in  part  by  the  ONR  under  Grant  No.  N00014-04-1-0617. 
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1  INTRODUCTION 


Unmanned  surface  vehicles  (USVs)  are  envisioned  to  be  a  multi-mission,  rapidly  configurable, 
long  endurance  unmanned  platform  providing  offboard  sensor  and  weapon  capability  to  perform 
critical  missions  like  intelligence,  surveillance,  reconnaissance,  anti-submarine  Warfare,  mine 
countermeasures,  communications,  and  navigation.  The  main  technical  challenges  during  these 
USV  missions  are  to  provide  optimal  performance  while  being  subjected  to  adverse  motions, 
multi-vehicle  control  issues,  data  fusion  challenges  and  obstacle  avoidance  maneuvering.  Most 
importantly,  in  an  autonomous  environment,  vehicle  communications  are  likely  subjected  to  fad¬ 
ing,  jamming  and  interference,  thus  subjecting  the  vehicle  to  robust  and  reliable  communication 
challenges. 

In  this  paper,  we  consider  the  propagation  environment  and  technology  of  multi-sensor  com¬ 
munications  in  unmanned  surface  vehicle  networks  located  in  a  rough  sea  environment.  We  first 
investigate  the  statistical  model  of  the  propagation  environment,  where  we  provide  the  analyt¬ 
ical  expressions  of  the  wave  height,  the  probability  of  having  line-of-sight  (LOS)  between  two 
USVs  in  terms  of  the  weather  condition  and  antenna  height,  and  the  spatial  correlation  of  the 
propagation  channels  related  to  different  antennas.  We  then  investigate  the  fading  reduction 
capability  of  a  multiple-input-multiple-output  (MIMO)  system  and  the  performance  in  general 
Ricean  channels  with  different  channel  correlations  is  examined.  The  results  are  important  to 
the  understanding  of  the  communication  environment  where  USV  communications  take  place. 

The  offerings  and  design  architectures  of  multi-sensor  communication  methodologies  are  then 
discussed.  In  combating  channel  impairment  problem  caused  from  multipath  fading,  diversity 
and  array  processing  methods  have  been  effectively  used  for  many  years.  These  methods  improve 
the  communication  quality  and  enhance  the  communication  quality  without  increasing  the  total 
transmit  power  or  consuming  additional  bandwidth  resources.  While  traditional  array  process¬ 
ing  and  diversity  techniques  use  multiple  antennas  at  the  receiver  and  a  single  antenna  at  the 
transmitter  (referred  to  as  single-input-multiple-output,  or  SIMO,  system),  the  development  of 
space-time  coding  techniques  also  permits  the  use  of  multiple  transmit  antennas  as  well,  result¬ 
ing  in  more  effective  MIMO  systems.  We  pointed  out  that,  depending  on  the  communication 
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characteristics,  particularly  the  statistical  distribution  and  spatial  correlation  characteristics, 
different  multi-sensor  techniques  may  be  preferred  in  different  fading  environments  to  achieve 
robust  and  reliable  USV  networking  connectivity.  Analytical  as  well  as  numerical  results  will  be 
provided  to  demonstrate  the  effectiveness  of  the  multi-sensor  techniques  for  various  communi¬ 
cation  scenarios. 

This  paper  is  organized  as  follows.  Section  2  discusses  the  wave  height  and  period  dis¬ 
tributions,  and  their  relationship  to  the  weather  status  is  examined.  Section  3  considers  the 
probability  of  sea  status.  The  probability  of  having  LOS  is  investigated  in  Section  4.  Section  5 
examines  the  space-time  correlation  coefficients  between  propagation  channels  related  to  differ¬ 
ent  antennas.  Section  6  considers  the  structure  and  advantages  of  MIMO  systems  in  combating 
multipath  fading  and  the  bit  error  rate  (BER)  performance  is  examined  numerically. 

2  WAVE  HEIGHT  AND  PERIOD  DISTRIBUTIONS 

In  this  section,  we  review  the  wave  height  and  period  distributions  in  a  sea  surface.  Wave 
height  and  wave  period  are  highly  correlated  random  variables  [Och98].  Therefore,  in  order  to 
study  the  wave  distribution  we  need  to  consider  the  joint  probability  of  the  wave  height  and 
the  wave  period.  The  simplest  case  is  when  there  is  a  fully  developed  sea,  which  is  described 
by  the  Pierson-Moskowitz  (P-M)  spectrum  [PM64],  In  this  case,  there  is  only  one  crest  and  one 
trough  for  each  half-cycle  of  the  wave.  Therefore,  the  wave  spectrum  is  generally  considered  to 
be  narrowband.  In  order  to  simplify  the  problem,  we  introduce  the  following  assumptions: 

a.  Waves  are  considered  to  be  a  steady-state  ergodic  random  process; 

b.  Waves  are  a  Gaussian  random  process  with  zero  mean  and  a  variance  representing  the  sea 
severity; 

c.  The  wave  spectral  density  function  is  narrowband; 

d.  Wave  peaks  and  troughs  are  statistically  independent. 
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2.1  Wave  Spectrum 


The  P-M  spectrum  is  defined  as  [PM64] 

where  w  is  the  wave  frequency,  g  =  9.81  m/s 2  is  the  gravitational  acceleration,  U  is  the  wind 
speed  measured  at  19.5  m  above  the  sea  surface,  and  a  =  8.10  x  10~3  and  b  =  0.74  are  constants. 
The  modal  frequency  w  and,  consequently,  the  wave  period  T  and  the  wave  length  L,  are  given 
by  a  function  of  wind  speed  as 


w 

=  0.87  (g/U) 

(2) 

2itU 

(3) 

T 

w  0.87 g 

L 

27 xg 
w 2  ’ 

(4) 

It  is  evident  from  Eqs.  (1)— (4)  that  the  P-M  spectrum  depends  only  on  the  wind  speed  U. 


Usually,  the  wave  spectrum  and  the  sea  state  are  described  in  terms  of  significant  wave  height 
H. s,  instead  of  the  wind  speed  U.  The  significant  wave  height  Hs  is  defined  as  the  average  height 
of  1/3  of  the  highest  waves  and  is  computed  through  the  following  procedure:  Measure  the 
wave  height  for  a  few  minutes  and  select  90  wave  crests  and  record  their  heights;  Then  choose 
30  highest  waves  and  calculate  the  average  height.  A  more  useful  definition  of  II s ,  however,  is 
given  by  the  integration  of  the  spectral  density  function  given  in  Eq.  (1),  yielding, 

f°°  a  r/4 

l  s^dw  =  T^-  <*> 

Moreover,  by  assuming  that  the  wave  spectrum  is  narrowband,  the  area  under  the  spectral 
density  function  is  equal  to  [Och98] 


Combining  Eqs.  (5)  and  (6),  the  following  relationship  can  be  derived  between  the  wind  speed 
U  and  the  significant  wave  height  Hs, 


Hs  =  2^/aJb(U2jg)  =  0.21  U2/g. 


(7) 
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Using  Eqs.  (1)  and  (7),  the  P-M  spectrum  can  be  expressed  in  terms  of  the  significant  height 
as 

SW  -  8.10  x  l<r3|lexp  |  -  O-rcytfj.  (8) 

Fig.  1  depicts  the  significant  wave  height  Hs  and  the  average  period  T  versus  the  wind  speed 
U  for  the  P-M  spectrum  as  given  by  Eqs.  (3)  and  (7). 

2.2  Wave  Height  and  Wave  Length 


While  the  wave  spectrum  gives  a  general  result  of  the  sea  state,  it  does  not  fully  describe  it.  As 
we  mentioned  earlier,  the  wave  dispersion,  period  and  height  are  not  statistically  independent. 
For  example,  the  higher  the  waves  the  longer  the  period.  Therefore,  we  need  to  consider  the 
joint  probability  density  function  between  these  variables. 


Consider  fully  developed  seas,  which  are  assumed  by  the  P-M  spectrum,  and  we  assume  that 
waves  are  narrowband  Gaussian  random  processes  whose  spectrum  is  highly  concentrated  in  the 
vicinity  of  the  modal  frequency  w.  To  examine  the  bandwidth  of  the  wave  spectrum,  we  define 
a  spectrum  width  parameter  v  as  [LH83] 


/mom2  jN1/2  [&}_ 

V  ml  )  \J  no  w’ 


where 


poo 

rrij  =  /  up  S{w)dw 
Jo 

is  the  jth  moment  of  the  spectrum,  and  Ho  and  /x2  are  expressed  as 


(9) 


(10) 


respectively. 


poo 

Ho  =  I  S(w)dw  —  mo, 
Jo 


poo 

H  2  =  /  (t»  —  w)2S(w)dw  =  (rooTO2  —  m\)/mo, 

Jo 


(11) 

(12) 


For  a  general  wave  spectrum  defined  as 

S(v>)  =  ^exp(-B/«;4), 


(13) 
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the  moments  can  be  expressed  in  terms  of  A  and  B  as  [Och98] 


mo  = 
mi  = 

m2  = 
m  4  = 


A_ 

45’ 


0.306 

A 


A 

£3/4’ 

A 


4  yfB' 


00. 


(14) 

(15) 

(16) 
(17) 


Note  that  the  P-M  spectrum  is  a  special  case  of  the  model  depicted  in  Eq.  (13)  where  A  =  ag 2 
and  B  =  b(g/U)A.  By  substituting  Eqs.  (14)-(16)  into  Eq.  (9),  the  width  of  the  wave  spectrum 
defined  by  Eq.  (13)  is  found  to  be 

A  A 


„=  (M±M  X»  (  X* 


0.43. 


(18) 


Eq.  (18)  shows  that  the  parameter  v  is  a  constant  in  this  case,  which  is  in  agreement  with  the 
normalized  P-M  spectrum,  where  the  frequency  axis  is  normalized  by  its  modal  frequency  w 
[Och98,  p.  34].  Moreover,  since  v  <  0.6  the  narrowband  assumption  remains  valid  [LM83]. 

2.3  Joint  Probability  Function 


Define 


as  the  dimensionless  wave  height,  and 


H 


^  ~  T 


(19) 


(20) 


as  the  dimensionless  period,  and  let  T  represent  the  mean  (average)  wave  period.  The  joint 
probability  density  function  (pdf)  of  wave  height  and  period  in  dimensionless  form  is  given  by 
[Och98] 

/(C’,)  =  ^(1  +  T)?eXP{-y[1  +  (1-^)  ^]}'  0«<oo.0<fl<oo.  (21) 

Since  the  P-M  spectrum  has  a  constant  wave  bandwidth  as  given  by  v  in  Eq.  (18),  f((,  rj)  is 
independent  of  the  wind  speed.  Fig.  2  depicts  the  value  of  /(£,  rj)  versus  the  dimensionless 
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wave  height  and  dimensionless  wave  period.  Moreover,  the  marginal  pdf  of  wave  height  can  be 
obtained  from  Eq.  (21)  as 


/(0=J f  /(C^  =  Cexp(-C2/2)^l  +  ^$(C/v),  0 


<  £  <  oo, 


(22) 


where 


$(«)  =  “/==  /  exp(— tx2/2)dn.  (23) 

\  Z7T  J— oo 

Note  that,  for  narrowband  wave  spectrum,  i.e.,  v  <  0.6,  the  value  of  (1  +  (v2/4))$(C/v)  in 
Eq.  (22)  is  very  close  to  unity. 

The  most  probable  wave  height  and  wave  period  are  found,  respectively,  from  df((,  rj)/dC,  =0 
and  9/(C,  v)/^  =  0,  and  are  expressed  as 


and 


C  = 


V  = 


1  +  vi 


1 


1  +  V2' 

The  wave  period  T  can  be  expressed  in  terms  of  the  wave  frequency  by 

T  —  2n/w. 


(24) 


(25) 


(26) 


Therefore,  the  dimensionless  wave  period  can  be  expressed  in  terms  of  the  wave  frequency  as 

(27) 


T  2-n/w  w 
V  =  7¥,~ 


T  2  Txjw  w 

Moreover,  using  the  direct  relationship  between  the  wave  frequency  and  the  associated  wave 
length  expressed  in  Eq.  (4),  the  dimensionless  wave  period  can  be  expressed  as 


w  =  V^ng/L  =  fl 
w  y/27 vq/L  V  L 


V2ng/L 

3  PROBABILITY  OF  SEA  CONDITIONS 


(28) 


In  the  previous  section,  we  have  discussed  the  joint  pdf  of  the  wave  height  and  wave  length 
which  provides  a  general  view  of  the  probability  that  a  certain  scale  of  wave  occurs  in  a  given 
sea  state.  In  order  to  fully  describe  the  situation  in  a  probabilistic  sense,  we  introduce  in  this 
section  the  probability  of  the  sea  states  that  a  USV  is  subjected  to  [Och98]. 
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In  general,  the  sea  severity  is  a  function  of  the  frequency  of  occurrence  of  storms,  water  depth, 
wind  direction,  etc.  Thus,  there  is  no  scientific  basis  for  selecting  a  specific  probability  distri¬ 
bution  function  to  represent  the  statistical  properties  of  the  sea  state  (or  the  significant  wave 
heights).  Nevertheless,  data  collected  over  the  years  has  yielded  many  different  distributions. 
It  is  generally  accepted  that  a  great  part  of  the  significant  wave  height  data  is  well  represented 
by  the  log-normal  probability  distribution  [Och78],  However,  the  Weibull  probability  distrib¬ 
ution  is  accepted  when  waves  are  high  [BS86,  MB90].  The  generalized  gamma  distribution  is 
introduced  by  [Och98]  which  has  the  advantage  of  representing  both  the  log-normal  and  Weibull 
distributions.  The  pdf  f{x )  and  the  cumulative  distribution  function  F(x)  of  the  generalized 
gamma  distribution  are  given  by 


where  0  <  x  <  oo, 


/(*)  =  \C"V’”-1  rapMAan 

(29) 

F(x)  =r{m,(Az)c}/r(m), 

(30) 

roc 

r  (m)  =  /  xm~le~xdx, 

Jo 

(31) 

T(m,x)  =  [  xm~le~xdx 

Jo 

(32) 

are  gamma  functions,  and  m,  c,  and  A  are  constants  dependent  to  the  actual  sea  environment. 


Table  1  shows  the  data  collected  in  the  North  Sea  over  a  period  of  three  years  [Bou78]. 
For  this  set  of  data,  the  parameters  for  the  generalized  gamma  distribution  are  estimated  as 
m  =  1.60,  c  =  0.98  and  A  =  1.37.  Both  f(x)  and  F(x)  are  plotted  in  Fig.  3.  In  this  case,  there 
is  a  5%  probability  that  the  significant  wave  height  is  larger  than  3  m  (Hs  >  3  m). 

4  PROBABILITY  OF  HAVING  LOS 

In  this  section,  we  consider  the  problem  whether  two  USVs  that  are  in  a  rough  sea  environment 
will  have  an  LOS  communication.  The  existence  or  absence  of  the  LOS  is  a  key  factor  that 
affects  the  propagation  channel  characteristics. 

Depending  on  the  location  of  the  two  USVs,  we  can  classify  the  communication  scenario  into 
two  categories:  long-distance  and  short-distance  communications.  Refer  to  Fig.  4,  the  former 
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Table  1:  Significant  wave  heights  [Bou78]. 


Significant  wave 
height  (m) 

No.  of  observations 

0-0.5 

1280 

0.5-1.0 

1549 

1.0-1. 5 

1088 

1. 5-2.0 

628 

2.0-2.5 

402 

2. 5-3.0 

192 

3.0-3. 5 

115 

3.5-4.0 

63 

4.0-4.5 

38 

4. 5-5.0 

18 

5.0-5. 5 

21 

5.5-6.0 

7 

6.0-6.5 

8 

6.5-7.0 

2 

7.0-7.5 

1 

Total  5412  in  3  years 

occurs  when  the  distance  between  the  two  USVs,  La,  is  much  larger  than  the  most  probable 
wave  length,  L.  The  latter  occurs  when  the  USVs  are  closer  and  their  distance  is  comparable 
to  the  most  probable  wave  length.  In  this  paper,  we  only  consider  long-distance  scenarios. 

4.1  Single- Wave  Scenario 

We  first  consider  a  scenario  where  only  one  significant  wave  is  present  between  the  two  USVs. 
Denote  He  as  the  height  of  the  wave  where  a  USV  is  located,  and  Hw  as  the  height  of  the 
significant  wave.  The  pdf  of  both  He  and  Hw  can  be  expressed  by  the  following  formula 

/(H)  =  ^exp|-2^y|,/f>0,  (33) 

where  Hs  is  the  significant  wave  height.  Assume  that  the  heights  He  and  Hw  are  independent, 
then  their  joint  pdf  is  expressed  by 

/(ft.ff„)  =  ^exp{-2^2|^exp|-2^y|,  He>0,Hw>0.  (34) 

Denote  Ha  as  the  antenna  height  relative  to  the  wave  level,  then  the  height  of  the  antenna  is 
Ha  +  He.  For  simplicity  of  analysis,  we  consider  that  the  probability  of  having  LOS  between  the 
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two  USVs  equals  to  the  probability  that  the  antennas  of  both  USVs  are  higher  than  the  peak 
of  the  significant  wave.  The  probability  that  an  antenna  is  higher  than  the  significant  wave  is 
given  by 


<■00  rHa+He 


P(Ha+He>  Hw)  =  f°°  [  f(He>Hw)dHwdHe 

Jo  Jo 

foo  rHa+He  2 t  /  it  \  2 


ni 


Hs 


■exp 


d/ffiud/ffe  ♦ 
(35) 


Therefore,  the  probability  of  having  LOS  between  the  two  USVs  is  a  function  of  the  antenna 
height  Ha  and  the  significant  wave  height  Hs,  and  is  expressed  as 

P[LOS\  =  {P(Ha  +  He>  Hw)}2 ,  (36) 

where  the  square  operation  emphasizes  the  fact  that  the  antennas  of  both  USVs  have  to  be 
higher  than  the  significant  wave.  The  results  for  two  different  sea  states  represented  by  Hs  =  2 
m  and  H s  =  3  m  are  depicted  in  Fig.  5.  It  is  evident  that  the  probability  of  having  LOS  is  very 
sensitive  to  the  antenna  height  Ha  as  well  as  the  significant  wave  height  Hs. 

4.2  Multi- Wave  Scenario 

Now  we  consider  the  situation  where  multiple  significant  waves  are  present  between  the  two 
USVs.  Since  the  significant  waves  constitute  of  only  one-third  (1/3)  of  the  waves,  the  number 
of  significant  waves  between  the  two  USVs  is  determined  from  the  ratio  between  the  distance 
between  the  boats,  La,  and  the  most  probable  wave  length,  L,  as 

w-se-  (37) 

where  the  most  probable  wave  length  L  is  given  by  Eq.  (24).  By  using  v  —  0.43  as  derived  in 
(18),  we  have 

”=  \/f=TTibp=0'844’ 

f  _  2irg  __  27r Hs 
0A2g/Hs  ~  0.42  ' 
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From  these  equation,  we  can  express  L  in  terms  of  Hs  as 


L  =  0.8442L  = 


0.844227rifs 

0.42 


Assume  that  the  N  significant  waves  are  statistically  independent,  we  can  express  the  prob¬ 
ability  of  having  LOS  as 

P\LOS\  =  {P(Ha  +  He>  Hw)}2N  ,  (39) 


Fig.  6  shows  the  P[LOS ]  for  different  distances  between  the  two  boats  defined  by  the  pa¬ 
rameter  N  and  different  sea  states  defined  by  Hs  —  2  m  and  Hs  =  3  m,  respectively.  The 
probability  of  having  LOS  decreases  as  the  number  of  significant  waves  increases. 

5  SPATIAL  CORRELATION 


In  this  section,  we  consider  the  spatial  correlation  between  the  propagation  channels  related 
to  two  spatially  separated  antennas.  We  first  review  the  general  concept  of  spatial  correlation 
developed  in  a  wideband  Nakagami-Rice  fading  environment.  Denote  Ax  as  the  distance  between 
two  receive  antennas,  denoted  as  Rxl  and  Rx2  in  Fig.  7,  and  B  as  the  bandwidth  of  the  signal. 
Then  the  wideband  spatial  correlation  function  is  given  by  [KI00] 

/—  r--f 

PSp{Ax)  =  — \ - ,  (40) 

f*  ['  {K/x(A/,0;/)|2}dA/d/ 

J  2  J  2  7 

where  pa  jx(A /,  Ax;  /)  is  the  space-frequency  correlation  coefficient  between  the  two  receive 
antennas  with  a  frequency  difference  A/.  When  the  time  delay  and  the  angle-of-arrival  (AO A) 
of  a  particular  wave  are  statistically  independent,  pa>fx(Af,Ax-,f)  can  be  decoupled  into  the 
frequency  correlation  coefficient  paj{Af)  and  the  spatial  correlation  coefficient  pajX(Ax;  /),  i.e., 

Pa,fx(Af ,  Aa3  /)  »  Pa,f(Af)pa,x(Ax'’  /)•  (41) 

Substituting  Eq.  (41)  into  Eq.  (40)  simplifies  the  double  integral  into  two  single  integrals. 
Denote  ps(r)  as  the  normalized  power  delay  profile.  Then,  paj(Af)  is  given  by  [Cla68] 

/OO 

Ps(t)  exp( — j  2-7T  A/r)  dr.  (42) 

■OO 
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Similarly,  the  spatial  correlation  coefficient  pa,x(Ax;  f)  can  be  obtained  as  [KI00] 

ft,„(As;/)  -  £ nWexp  (43) 

where  £1(8)  is  the  normalized  angular  power  profile,  and  C  is  the  velocity  of  light. 

^  - 

5.1  Proposed  Model 

We  now  propose  a  model  to  evaluate  the  value  of  pax(Ax\  /),  whose  general  expression  is  given 
in  (43).  In  a  rough  sea,  the  sea  surface  can  be  described  as  a  combination  of  a  large  number 
of  small  waves  and  a  small  number  of  high  waves.  We  consider  the  case  that  a  high  wave  with 
height  h  obstructs  the  LOS  between  the  two  communicating  USVs,  whereas  a  number  of  random 
smaller  waves  surrounding  both  USVs  are  modeled  using  the  ring  scattering  model.  Therefore, 
the  propagation  from  a  transmit  antenna  to  a  receive  antenna  experiences  scattering  as  well  as 
diffraction. 


In  considering  the  diffraction  effect,  the  diffration  edge  in  the  top  of  a  high  wave  is  assumed 
to  have  an  angular  power  profile  £1(9)  which  is  composed  of  two  components:  (i)  the  diffraction 
component  which  is  assumed  to  be  uniformly  distributed  about  a  mean  angle  9b  with  a  diffraction 
width  /?,  and  (ii)  the  scattering  component  which  is  assumed  to  be  uniformly  distributed  over 
[0,2tt].  We  define  S  as  the  strength  between  the  scattered  and  diffracted  components  of  the 
received  signals.  The  composite  angular  profile  for  multipath  waves  is  then  given  by 


f  1  S 

(0(1  +  5)  +  2tt(1  +  5)’ 


m  =  < 


if  \9\  < 


2’ 


(44) 


_ S— 

12tt(1  +  5)’ 


otherwise. 


Therefore,  the  new  spatial  correlation  coefficient  is  obtained  by  combining  Eqs.  (43)  and  (44), 


i.e., 


pa,x{ Arc;  /) 


r2-K 

Jo 


+ 


2jt(1  +  S) 

Ob+P/2  1 


m- 

Jeh- 


exp  [jkAxcos(9  —  9b)\  dB 
exp  [jkAxcos(9))  d9, 


(45) 


I 6 b-0/2  P0-  +  S) 

where  k  =  2nf/C.  The  first  integral  in  Eq.  (45)  is  independent  of  9f,  and  is  the  Oth  order 
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In  this  section,  we  consider  MIMO  communications  in  rough  sea  environment.  The  Ricean 
channel  model  is  first  described,  and  the  MIMO  performance  is  then  investigated.  Cooperative 
diversity  is  also  an  alternative  technology  in  rough  sea  environment  [ZWW04]. 

6.1  Ricean  Channel  Model 


It  is  well  known  that  the  general  model  of  the  envelope  (magnitude)  of  a  communication  channel 
is  well  represented  by  the  Ricean  distribution  [UJOO,  Rap02],  The  Ricean  distribution  of  an 
envelope  is  the  result  of  a  complex  Gaussian  channel  with  non-zero  mean.  The  mean  value 
usually  corresponds  to  the  determinant  component  due  to  direct  path  propagation  in  the  presence 
of  LOS,  whereas  the  random  component  is  due  to  multipath  scattering. 


The  Ricean  distribution  of  the  envelope  x  is  given  by 

(  x  f  x2  +  A2\  ( xA\ 

p(x)  =  <  (T%  exp  \  2og  J  °W  / 

u 


x  >  0, 
x  <  0, 


where  A  is  the  peak  amplitude  of  the  determinant  component  and  Jo(-)  is  the  modified  Bessel 
function  of  the  first  kind  and  zero-order.  The  Ricean  distribution  is  often  described  in  terms  of 
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a  parameter  K  which  is  defined  as  the  ratio  between  the  determinant  component  power  and  the 
variance  of  the  random  components,  i.e., 


a*_ 

2a f 


Note  that,  the  total  power  of  a  Ricean  distributed  envelope  is 


P  =  A2+2al 


(50) 


(51) 


In  particular,  when  A  =  0,  i.e.,  there  is  no  determinant  component,  the  Ricean  distribution 
reduces  to  Rayleigh  fading. 

6.2  MIMO  Performance 


MIMO  technology  has  ignited  a  new  revolution  in  information  communications.  By  using  mul¬ 
tiple  antennas  at  the  two  ends  of  a  communication  link,  a  high  diversity  gain  is  achieved,  and 
the  channel  capacity  between  a  transmitter  and  a  receiver  in  a  fading,  scattering  environment 
grows  linearly  with  the  number  of  antennas  [Fos96,  Tel99,  FG98,  TSC98]. 

In  this  paper,  we  consider  a  simple  setting  where  each  USV  is  equipped  with  two  antennas. 
The  channels  between  the  two  transmit  and  two  receive  antennas  can  be  denoted  using  the 
following  matrix 

where  hltk  is  the  complex  coefficient  of  the  channel  connecting  the  ith  transmit  antenna  and  the 
kth  receive  antenna,  i,k  =  1,2. 

As  we  discussed  earlier,  the  channels  are  in  general  Ricean  distributed  and  have  correlation 
between  them.  In  a  Ricean  channel,  as  the  value  of  K  increases,  the  random  component  becomes 
less  significant  and,  therefore,  the  correlation  between  the  channels  becomes  less  significant.  Note 
that  the  determinant  components  are  time-invariant  and,  therefore,  are  coherent. 

The  performance  of  the  MIMO  system  is  evaluated  for  different  sea  scenarios.  The  Alamouti 
space-time  code  [Ala98]  with  quadrature  phase  shift  keying  (QPSK)  modulation  is  used  in 
transmitting  the  MIMO  signals.  Fig.  9  shows  the  bit  error  rate  (BER)  performance  for  different 
values  of  K  =  0, 1,  and  5,  whereas  for  each  plot,  the  spatial  correlation  coefficient  p  between  the 


101 


two  antennas  takes  values  of  0,  0.5,  0.7,  and  1.  When  p  =  0  the  channel  are  uncorrelated  and 
when  p  =  1  the  channels  are  fully  correlated  (coherent).  The  results  clearly  show  the  advantage 
of  the  MIMO  system  compared  to  the  corresponding  single-antenna  system,  particularly  in 
low  spatial  correlation  scenarios.  When  the  channels  are  coherent,  the  MIMO  only  achieves 
3  dB  array  gain  compared  to  the  single  antenna  case  (Note  that  Alamouti  code  assumes  no 
channel  state  information  at  the  transmitter  and,  therefore,  no  beamforming  gain  is  achieved  at 
the  transmitter).  Therefore,  it  becomes  clear  that  higher  diversity  gain  is  achieved  with  lower 
spatial  correlation,  which  underscores  the  importance  of  having  low  channel  correlation  in  a 
fading  environment. 

7  CONCLUSIONS 

In  this  paper,  we  have  considered  the  application  of  multiple-input-multiple-output  (MIMO) 
technology  in  networked  communications  between  unmanned  surface  vehicles  (USVs)  in  a  rough 
sea  environment.  Specifically,  we  have  taken  a  detailed  investigation  in  the  channel  modeling, 
which  includes  the  wave  height,  period,  and  length.  The  probability  of  having  line-of-sight 
between  two  USVs  is  derived,  and  its  relationship  to  the  antenna  height  and  significant  wave 
height  is  revealed.  We  then  developed  the  spatial  correlation  between  the  propagation  channels 
associated  to  two  antennas  in  a  rough  sea  environment.  Finally,  the  MIMO  technology  is  exam¬ 
ined  in  such  an  environment,  and  the  advantages  of  using  MIMO  technology  is  demonstrated. 
The  results  presented  in  this  paper  are  important  to  the  understanding  of  the  communication 
environment  where  USV  communications  take  place. 

REFERENCES 

[Ala98]  Alamouti,  S.  M.,  “A  simple  transmit  diversity  technique  for  wireless  communications,” 
IEEE  J,  Sel.  Areas  in  Comm.,  vol.  16,  no.  8,  pp.  1451-1458,  Oct.  1998. 

[Bou78]  Bouws,  E.,  “Wind  and  wave  climate  in  the  Netherlands  sectors  of  the  North  Sea  between 
53°  and  54°  north  latitude,”  Scientific  Report  W.  R.  Koninklijk  Netherlands  Meterologisch  Inst. 
78-9,  1978. 


102 


[BS86]  Burrows,  R.  and  Salih,  B.  A.,  “Statistical  modelling  of  long-term  wave  climates,”  Proc. 
20th  Int.  Conf.  Coastal  Engng.,  Taipei,  vol.  1,  pp.  42-56,  1986. 

[Cla68]  Clarke,  R.  H.,  “A  statistical  theory  of  mobile-radio  reception,”  Bell  System  Tech.  J.,  pp. 
957-1000,  July/Aug.  1968. 

[FG98]  Foschini,  G.  J.  and  Gans,  M.  J.,  “On  limits  of  wireless  communications  in  a  fading 
environment  when  using  multiple  antennas,”  Wireless  Personal  Communications,  vol.  6,  no.  3, 
pp.  311-335,  March  1998. 

[Fos96]  Foschini,  G.  J.,  “Layered  space-time  architecture  for  wireless  communications  in  a  fading 
environment  when  using  multi-element  antennas,”  Bell  Labs  Tech.  J.,  pp.  41-59,  1996. 

[KC97]  Kalkan,  M.,  and  Clarke,  R.  H.,  “Prediction  of  the  space-frequency  correlation  function 
for  base  station  diversity  reception,”  IEEE  Trans.  Vehi.  Tech.,  vol.  46,  no.  1,  pp.  176-184,  1997. 

[KI00]  Karasawa  Y.  and  Iwai  H.,  “Formulation  of  spatial  correlation  statistics  in  Nakagami-Rice 
fading  environments,”  IEEE  Trans.  Antennas  Propagat.,  vol.  48,  no.  1,  pp.  12-18,  2000. 

[LM83]  Longuet-Higgins,  M.  S.,  “On  the  joint  distribution  of  wave  periods  and  amplitudes  in  a 
random  wave  field,”  Proc.  Roy.  Soc.,  London,  Ser.  A,  289,  pp.  241-58,  1983. 

[MB90]  Mathiesen,  J.  and  Bitner-Gregersen,  E.,  “Joint  distribution  for  significant  wave  height 
and  zero-crossing  period,”  Appl.  Ocean  Res.,  vol.  12,  no.  2,  pp.  93-103,  1990. 

[Och78]  Ochi,  M.  K.,  “On  long-term  statistics  for  ocean  and  coastal  waves,”  Proc.  16th  conf. 
Coastal  Engng.,  Hamburg,  vol.  1,  pp.  59-75,  1978. 

[Och92]  Ochi,  M.  K.,  “New  approach  for  estimating  the  severest  sea  state  from  statistical  data,” 
Proc.  23rd  Int.  Conf.  Coastal  Engng.,  Venice,  vol.  1,  pp,  512-25,  1992. 

[Och98]  Ochi,  M.  K.,  Ocean  Waves:  The  Stochastic  Approach,  Cambrige,  1998. 

[PM64]  Pierson,  W.  J.  and  Moskowitz,  L.,  “A  proposed  spectral  form  for  fully  developed  wind 


103 


seas  based  on  the  similarity  theory  of  S.  A.  Kitaigorodskii,”  J.  Geophys.  Res.,  vol.  69,  no.  24, 
pp.  5181-80,  1964. 

[Rap02]  Rappaport,  T.  S.,  Wireless  Communications:  Principles  and  Practice,  Second  Edition, 
Prentice-Hall,  2002. 

[Tel99]  Telatar,  I.  E.,  “Capacity  of  multi-antenna  Gaussian  channels,”  European  Transactions 
on  Telecommunications,  vol.  10,  no.  6,  pp.  585-595,  Nov.-Dee.  1999. 

[TSC98]  Tarokh,  V.,  Seshadri,  N.,  and  Calderbank,  A.  R.,  “Space-time  codes  for  high  data  rate 
wireless  communication:  Performance  criterion  and  code  construction,”  IEEE  Trans.  Inform. 
Theory,  vol.  44,  pp.  744-765,  1998. 

[UJ00]  Ungan  B.  U.  and  Johnson  J.  T.  “Time  statistics  of  propagation  over  the  ocean  surface: 
A  numerical  study,”  IEEE  Trans.  Geoscience  and  Remote  Sensing,  vol.  38,  no.  4,  pp.  1626-1634, 
2000. 

[ZWW04]  Zhang,  Y.,  Wang,  G.,  Watson,  M.,  and  Amin,  M.,  “Cooperative  wireless  network  for 
unmanned  surface  vehicles,”  AUVSI  Unmanned  Systems,  Anaheim,  CA,  Aug.  2004. 


104 


Figure  1:  Significant  wave  height  versus  wind  speed. 


Figure  2:  Significant  wave  height  versus  dimensionless  wave  height  and  dimensionless  wave 
period. 
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P[LOS]  for  a  single  significant 


Figure  6:  Probability  of  having  LOS  P[LOS\  for  multiple  significant  waves. 


Figure  7:  Illustration  of  the  communication  environment. 
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(a)  K  =  0 


(b)  K  =  1 


(c)  =  5 


Figure  9:  BER  performance  for  a  MIMO  system  (2  transmit  antennas,  2  receive  antennas). 
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